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If light nuclei are approximated by the alpha-model, with the particles localized in clusters 
of four (or less), the particle with uncompensated spin is accelerated by the relatively weak 
attraction of alphas, or by the difference of the attractions of two clusters. Then the Thomas 
relativistic term is not in general larger than the Larmor magnetic term in the spin-orbit coup- 
ling in contrast to the result of the central model, wherein the acceleration may be attributed 
to unsaturated forces toward one center. The relative magnitudes of these terms are estimated 
for the alpha-model of Li’ and C", two of the few nuclei to which it is here considered that the 
model might apply. The magnetic moment consistent with the Larmor-Thomas coupling in the 
alpha-model of B" is found to be considerably larger than the experimental value. The quad- 
rupole moment of the deuteron implies a spin-orbit coupling arising directly from angle- 
dependent nuclear forces. The appropriate ‘‘spin-orbit-spin’”’ angle dependence suggested by 
the meson theory causes no first-order coupling in the alpha-model of Li’, however, and a rough 
estimate of the second-order doublet splitting shows that it may be considerably smaller than 


that due to the Larmor and Thomas terms. 





OR nuclei which may be approximated by a 
central-field model, it has been pointed out! 
that the Thomas relativistic term may be ex- 
pected to be considerably larger than the Larmor 
magnetic term, and to determine the sign and 
order of magnitude of the spin-orbit coupling 
due to spherically symmetric forces between the 
particles in light nuclei. This may not be said in 
general for nuclei which may be approximated 
by the alpha-model, in which the protons and 
neutrons are treated?: * as being clustered into 
alphas and possibly an extra particle or a triton 
or He*. The difference arises largely from the 
fact that the Thomas term contains the acceler- 
ation of the particle with uncompensated spin 
toward the center of mass of the system, and in 
1D. R. Inglis, Phys. Rev. 50, 783 (1936). 


?L. R. Hafstad and E. Teller, Phys. Rev. 54, 681 (1938). 
3H. Bethe, Phys. Rev. 53, 842 (1938). 


the alpha-model this acceleration is due to 
inter-alpha-forces which are supposedly con- 
siderably weaker than intra-alpha-forces, in 
keeping with the saturation requirement. 

The validity of explaining nuclear doublet 
splittings by means of the Larmor-Thomas pre- 
cession depends on whether or not its effect is 
masked by the effect of the angle-dependent 
nuclear interaction responsible for the deuteron 
quadrupole moment. Such an interaction, of 
“spin-orbit-spin” angle dependence (s;-fr)(S2-T), 
is suggested by the meson theory, but specifi- 
cation of its radial dependence requires an 
arbitrary cut-off. The contribution of this 
interaction to the doublet splitting in the alpha- 
model, of Li’, for example, is a second-order 
effect and is expected to be quite sensitive to 
details of the model and of the radial dependence, 


1175 





1176 D. R. 


so it does not at present seem wise to try to 
compute the splitting. Since a very rough 
estimate below indicates that this contribution 
and the Larmor-Thomas splitting have about 
the same order of magnitude, an estimate of the 
Larmor and Thomas terms.remains physically 
interesting, especially in view of such minor 
successes as its explanation of the sign of the 
spin-orbit coupling in the relatively simple 
nuclei Li’ and O”. 

Among the (4n—1)-particle stable nuclei, Li’ 
is probably the only promising application of the 
alpha-model.* This is indicated first by the fact 
that the magnetic moment of Li’ calculated by 
the alpha-model* disagrees> with experiment 
considerably less severely than that calculated 
by a perturbation theory starting from the 
central model;* second, by the relatively high 
ratio of “internal” to ‘external’ binding of the 
triton in Li’; and third, by the failure of the 
alpha-model to furnish a correlation between the 
electrostatic energies of the various (4n—1)- 
particle nuclei.’ For the (4%+1)-particle nuclei, 
such as C!’, which do not contain an easily dis- 
torted triton, the use of the alpha-model seems 
more plausible. 

It is here shown that the magnetic term is at 
least of the same order of magnitude as the 
Thomas term in the alpha-model of the ground 
state of Li’ and of C", and in Li’ probably even 
larger. Only in the odd-neutron case C' (and 
Be’) can a preponderant magnetic term change 
the sign of the spin-orbit coupling, since the 


4 Compare, however, B.O. Gyhblom and R. E. Marshak, 
Phys. Rev. 55, 229 (1939). 

5D. R. Inglis, Phys. Rev. $5, 329 (1939). 

®D. R. Inglis, Phys. Rev /53, 882 (1938). 

7H. Brown and D. R. Inglis, Phys. Rev. 55, 1182 (1939). 
In that paper it appears ¥hat the electrostatic energy of a 
light nucleus is as much & about half as great as the binding 
energy between alphasf This fact furnishes an easy rough 
indication of the pr#ponderance of the magnetic term in 
the alpha-model if he spin stays with a triton; it favors the 
Thomas term byAnly a factor 2 in a comparison of nuclei 
with the well-kn6wn situation in atoms (wherein the electric 
field is the accelerating field, and the magnetic term is twice 
as large as the Thomas term). The magnetic term in nuclei 
is equally favored by the anomalously large gyromagnetic 
ratios of the protori and neutron, leaving the magnetic term 
preponderant as in atoms. Values of nuclear masses newer 
than used in that paper (compare W. H. Barkas, Phys. 
Rev. 55, 691 (1939)) alter the comparison with experiment 
but do not decrease the discrepancy. The following may be 
substituted in Table 1:  - 








H3—He3 Li7—Be? Bu-Cu Nb-OH 
a-model 0.86* 2.00 3.14* 4.28 
obs. 0.86 1.87 3.14 3.66 
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negative gyromagnetic ratio of the neutron here 
makes the magnetic and Thomas terms have 
opposite signs. It is further shown that the 
Larmor-Thomas splitting in B" is too small to 
lead to agreement of the result of the alpha- 
model with the observed magnetic moment. 


THE LARMOR-THOMAS PRECESSION IN Li’ 


For the sake of numerical simplicity in the 
dynamical problem we approximate the alpha- 
model of Li’ by two clusters of 3} particles each, 
having their center of mass at the geometrical 
center. Each cluster has a displacement r and a 
velocity v, relative to the center of mass, and a 
mass (7/2)M. The magnetic field at the triton 
due to the motion of the alpha relative to the 
triton is then 


H = (2e/c)(2r) X (2v)/(2r)?= (eL/7 Mcr'*), 


where L=7MrXv is the orbital angular mo- 
mentum. 
The Larmor energy is 


E.=—(ge/2Mc)S-H, 


where S is the spin angular momentum, (ge/2Mc)S 
the spin magnetic moment, and g the spin 
gyromagnetic ratio (5.6 for the proton in Li’ and 
about —4 for the neutron in Be’). S-L con- 
tributes a factor 3h?/2 to the doublet splitting, 
AE=E(?P1;2) —E(?P 3/2), so we have the Larmor 
splitting of the Li? doublet 


AE, = 3ge*h?/(28.M?cr*) 
= (3g/28)(137/1830)?(e?/mec*r) mc? 
= 0.003;(e?/mce*r)®mc?. 


The Thomas relativistic term is 
Er=S-aXv/2¢e 


the mean value of which is due mainly to the 
radial component of the acceleration a and the 
tangential component of the velocity v. We may 
consider the radial acceleration in two parts. 
First there is a centripetal acceleration a, 
accompanying the exchange of particles between 
the two clusters, and, second, the familiar 
centripetal acceleration a,=v*/r of a particle at 
rest in a rotating coordinate system. 

In treating the exchange, each particle may 
be described by a wave function having two 
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extremes, one at each of the clusters.? Each of 
these ‘molecular orbital’’ wave functions may be 
considered in the distant-cluster approximation 
to be a symmetric or antisymmetric combina- 
tion of the s functions of the two clusters, as in 
the perturbation treatment of the H2* ion. If we 
call the energy difference between the antisym- 
metric and symmetric states iw, the “frequency 
of exchange”’ between the clusters in the single- 
particle problem is w/27. The order of magnitude 
of the corresponding acceleration we take from 
simple harmonic motion, a,.~ —w*r. Using this 
and L=7MrxXv, we have 


AEr = (3/28M)(hw/c)?. 


The magnitude of fw is given in Fig. 4 of refer- 
ence 2 as iw=Q=mc’, but this depends upon the 
application of the alpha-model to B" and N'™® 
where it is probably not valid. A more generous 
estimate of iw may be made as follows: The 
binding energy of the alphas in Be® would be 
about 5mc* if there were no Coulomb energy,’ 
and the binding energy of the triton to the alpha 
in Li? would be about 7mc*. Of this, roughly 
(3/4)5mc?~4mc? might be due to the same sort 
of interaction as binds Be’, there being ¢ as 
many intercluster heavy-particle interactions in 
Li’, so we might attribute the remaining 3c? 
to ‘‘exchange’’—that is, to the fact that Li’ has 
one more particle in the low energy symmetric 
state than in the antisymmetric state, whilst in 
Be® both states are filled. Using hw=3mc*, we 
obtain 
AEr=(m/M)mce?=5 X 10-4 mc’. 


The “rotational” centripetal acceleration is 
a,= —(v/r)*r= —(L/7Mr’)*r 
and the consequent doublet splitting is 


AE?’ = (3/4) L?/ (7? M3cr*) = (3/2)h*/ (72 Mecr*) 
= 3X 10-*(e?/mc*r)*mc?. 


For the value of r suggested below, this is con- 
siderably larger than the previous term AEr. 
The contribution of the tangential component 
of a and the radial component of v to |aXv|« is 
((#/r)(d/dt)r*¢)», which does not in general van- 
ish in the alpha-model as it does in the central 
model. In the case of a simple harmonic motion 
along the figure axis (if one considers ¢ constant 





and so retains only the Coriolis acceleration), 
this becomes w’*r?¢, which is just |a.r¢|. The 
effect of the Coriolis term is thus just equal to 
the rather small effect AEr of the radial exchange 
acceleration, in this approximation. In the alpha- 
model of more complex nuclei, in which the 
particle is not so nearly constrained to move on 
a rotating line, the tangential acceleration would 
be expected to be even less important than here. 

The value of r to be used in the comparison of 
the various terms should be indicated by the 
mass difference of Li? and Be’, which is attributed 
to Coulomb energy according to the usual 
assumption that the Hamiltonian is otherwise 
symmetrical between protons and neutrons. The 
atomic mass difference Be’—Li’ is roughly® 1.0 
Mev=2.0mc, and that of He*—H®* is about 
zero. The Coulomb energy difference of Be’ — Li’ 
minus that of He*—H® is then about 2mc*® and 
the corresponding radius r= }ftq=}-2e*/(2mc*) 
=}e?/mc?. Using this value we obtain AE; 
=T7AE7'’~50AEr and a total splitting AE 
=0.03mc?=15 kev. Although these estimates 
are rough, they indicate that the magnetic term 
is almost certainly larger than the Thomas term 
in the alpha-model of Li’. The same is true of 
the alpha-model of Be’, the sign of its Larmor- 
Thomas coupling being thus altered and its 
ground state *P4;2. 

The fact that the ‘‘rotational’’ centripetal 
acceleration appears to be greater than the 
acceleration associated with exchange is an 
indication of the manner in which the alpha- 
model should be treated in rotational problems. 
Since the triton and alpha change places by 
rotation more frequently than they exchange 
identity by migration of a particle, it is more 
consistent to consider a triton and an alpha 
rotating about a center of gravity somewhat 
nearer the alpha than it is to picture two clusters 
of 3} particles each rotating about the geometrical 
center. This makes a slight difference*:* in the 
expected magnetic moment of Li’, favoring the 


8’ L. H. Rumbaugh, R. B. Roberts and L. R. Hafstad, 
Phys. Rev. 54, 657 (1939), especially p. 680. 

* R. G. Sachs, Phys. Rev. 55, 825 (1939). I wish to thank 
Dr. Sachs for discussion of his results, especially on C™, 
in connection with the present paper, and be pointing out 
that his reason for neglecting the Larmor term was a desire 
to make a formal comparison between the results of the 
two models, introducing the same assumptions, as nearly 
as possible, in both. 
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value 0.405 uy+u, of reference 3 over the value 
0.433 uv+u, of reference 9. 


THE LARMOR-THOMAS PRECESSION IN C!3 


C'’ has an extra neutron which, in the alpha- 
model, occupies free space on both sides of the 
plane containing the three alphas. It is attracted 
by the alphas and may be pictured as being 
accelerated continually toward the plane in 
which they lie. In the ground state, the nucleus 
rotates about an axis in the plane of the alphas 
with one unit of angular momentum.’ About this 
axis the moment of inertia is almost the same as 
that of the Li’ nucleus, if the alpha-alpha 
distance raa is nearly equal to the triton-alpha 
distance ria. The angular velocities are therefore 
about the same. The C' calculation differs from 
that for Li’ principally in the fact that the extra 
particle has a greater centripetal acceleration 
because it is not bound into a cluster with other 
particles. 

We can get a rough idea of the average distance 
of the neutron from the plane of the alphas in C™ 
by assuming that it is the same as (or slightly 
less than) that of the proton in N!*, and by 
using the observed Coulomb energy of the latter, 
3.21 mMU=6mc?. This is 6e?/rra, where fra 
is the average proton-alpha distance in N'*. We 
thus have for the average neutron-alpha distance 
in C!3, 7,.~e?/mc?. Taking raa~e?/mc*, as sug- 
gested by the Li’— Be’? Coulomb energy, we have 
the average position of the neutron on each side 
of the plane forming a regular tetrahedron with 
the alphas. The average radial component of the 
electric field at the neutron is only slightly 
reduced by the vectorial addition of the fields 
of the three alphas. This electric field is thus 


F = 6€/? xa? = 6m?c*/e8 


or somewhat less.. With the magnetic field 
H=vxXF/c and with vXr=L/2M, the Larmor 
splitting is roughly 


AE, = (9/4)g(137/1830)?mc? = —0.05mce?. 


The binding energy of the last neutron in C'® 
is about 10mc*. Although it moves in a non- 
Coulomb field so that we cannot apply the virial 
theorem exactly, we may assume as a rough 
approximation that the average kinetic energy 
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has the same magnitude. Of this, part corre- 
sponds to motion in the two dimensions parallel 
to the plane of the alphas, but probably as much 
or more corresponds to motion along the z axis. 
To calculate the partition would require quite 
detailed treatment of the forces, so we merely 
assume that about 5mc? corresponds to z motion. 
Approximating the z motion by a simple har- 
monic motion of frequency w/2m7 and amplitude 
r, and considering Er to be mainly due to this 
motion, we have 7= —w*z = —4(5mc?/Mr*)z and, 
due to the rotation, Utangentiat=(L/2MTaa’)z. 
These lead to a Thomas splitting 


AEr = (15/4) (137/1830)?mc? = 0.02mc? 


(or somewhat more with the Coriolis term). 
This is thus weaker than the magnetic effect, 
but so slightly that it remains for a more accurate 
estimate to decide definitely whether the spin- 
orbit coupling arising from the Larmor-Thomas 
precession is positive or negative in the alpha- 
model of C'*. This rough estimate serves at 
least to show that the Larmor term may not be 
merely neglected, and suggests that the coupling 
is more apt to be negative, in agreement with the 
central model!® and with experimental indica- 
tions." 


LARMOR-THOMAS COUPLING IN B" AND THE 
CONSEQUENT MAGNETIC MOMENT 


A study of Coulomb energies does not exclude 
the possibility that the alpha-model applies to 
B", but makes this seem unlikely by showing 
that the model does not apply both to B"™ and 
to N'®. Recent papers’: * suggest that the alpha- 
model of B", with the Larmor-Thomas coupling, 
stands in agreement with the observed magnetic 
moment, but this is due to apparently excessive 
caution in avoiding specification of parameters 
in the model: Even a relatively safe order-of- 
magnitude estimate of the spin-orbit coupling 
limits the plausible » to a narrow range which 
does not include the experimental value,” 
2.68 uv. The treatment of Li’ suffices to estimate 
the coupling due to body rotation and exchange 
acceleration (corresponding to 8 of reference 9) 


10 M.E. Rose and H. A. Bethe, Phys. Rev. 51, 205 (1937). 

1 C, H. Townes, Phys. Rev. 56, 850 (1939). 

22S. Millman, P. Kusch and I. I. Rabi, Phys. Rev. 56, 
165 (1939). 
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as 10-*mc?. An upper limit of the coupling due to 
rotation of the “hole’’ around the triangle (the 
parameter a of reference 9) is had by considering 
the entire axis-projection A= of the angular 
momentum to be due to rotation of one particle 
about the circumscribed circle. This gives an 
upper limit because a hole has less mobility 
(greater effective mass in the theory of metals) 
than a particle. Taking a radius 3~'e?/mc*, one 
finds the Thomas term AEr = }mc (or a= mc?/6) 
and the magnetic splitting much smaller. The 
ground state might differ from ?P3;2 by admixture 
of *D3;2. The rotational kinetic energy difference 
of the states is e = 3h?/2I ~ 7mc*, so the admixture 
of D wave function is of the order of magnitude 
imc?/7mc? =0.03. The operator g,L.+gs5, being 
diagonal, the magnetic moment yu of the ground 
state would thus be expected to differ from 
u(?? P32) =ue+5/11=3.23un by only about one- 
tenth of one percent, (0.03)*. (Specializing the 
formulas of reference 9 gives b2=3'a/16e and 
u=(1—10-)u(?P32). The lower limit 1.5 uy 
suggested in reference 9 and quoted in reference 
12 corresponds to about equal mixture of P and 
D states, which is quite implausible.) 


ANGLE-DEPENDENT NUCLEAR FORCES IN Li’ 


The above remarks have been directed toward 
an investigation of the possibility that nuclear 
properties may be derived from rather simple 
spherically symmetric exchange forces between 
the particles. The quadrupole moment of the 
deuteron indicates that such forces (if retained 
at all) must be supplemented by a torque between 
spin and position vectors, so as to mix virtual 
states with the low S state of the deuteron. The 
appropriate angle-dependent interaction, or 
“spin-orbit-spin coupling,” suggested by the 
meson theory, has the form 


H'= {8,°S2—3(s;-u)(Se-u) } f(r), 


where u=r/r, and r=r.—m, the difference of the 
position vectors of the particles. Because of the 
present difficulties in the meson theory," it still 
seems desirable to investigate the possibility of 
supplementing the old central exchange forces 
by a term H’ only strong enough to account for 
the deuteron quadrupole moment. We shall here 


13H. A. Bethe, Phys. Rev. 55, 1261 (1939). 
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estimate the magnitude of the spin-orbit coupling 
thereby introduced in the alpha-model of the 
relatively simple nucleus Li’. 

Considering H’ as a perturbation, the first- 
order spin-orbit coupling would be due to a 
difference of the ?P1/2 and *P 3/2 diagonal elements 
of H’. According to the sum rule, this may be 
calculated as the difference between the sum of 
the diagonal elements for the two states with 
(M1, Ms)=1, —}3, or 0, } and that for the state 
with (M,, Ms)=1, 4. (Here M, and Ms are the 
usual orbit and spin projection quantum num- 
bers.) The wave function of the ground state in 
the alpha-model, as in a molecule, is a product of 
a function of the internal coordinates, with Ms 
as a subscript, and a function of the orientation 
of the line between the clusters, with My, as a 
subscript. The function of internal coordinates 
may be written as an antisymmetric sum of 
products of the form 


= = _ i 
gs’ ts” gttst a’ta’ a*® : 


where s and a are the symmetric and antisym- 


metric ‘“‘molecular orbital’’ wave functions, 
orthogonal to one another, the superscripts v 
and 7 are the isotopic spin quantum numbers" 
referring to neutron and proton, respectively, 
and the superscripts + and — refer to single- 
particle spin projection m,, or to the spin func- 
tions a and 8, respectively. For calculating the 
matrix elements, the essential “‘spin-orbit-spin”’ 
coupling factor of H’ may be written 


(Si- r)(Se- r) - iL (sz+isy)1(S2—isy)2 
+ (s2—iSy)1(Sz+isy)2 ](x*?+y") 
+515 ,22?+terms which alter Ms, 


and the spin operators may be used in the 
form 


(9)-(0). ees(3)-() 
a/2 


a 


By) \-8/2)° 


The diagonal element of this coupling factor 
consists of direct terms, between the product 
written above and itself, and of exchange terms, 
between the product as written and a similar 
product with two of its factors interchanged, with 


Sz 


14 B, Cassen and E, U. Condon, Phys. Rev. 50, 846 (1936). 
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a negative sign because of the Pauli antisym- 
metry. One thus finds six direct terms in —2*/4 
from like-particle, unlike-spin pairs (two s—s, 
three a—s, and one a—a). There are also three 
terms in 2°/4 from like-particle, like-spin, a—s 
pairs, making the direct terms from a—s pairs 
vanish. Unlike-particle pairs contribute as many 
positive as negative terms. (x*?+y’) is absent from 
the direct terms because of orthogonality intro- 
duced by the accompanying spin operations. 
Among the exchange terms, we first have three 
terms —}(x*+y*) from the like-particle, s—s 
or a—a pairs (two s—s and one a—a). We may 
combine the terms involving s—s pairs into a 
space integral of —}(x?+y*+2*) f(r), independent 
of the sign of Ms, +. Due to its symmetry it is 
also independent of M/_;, and gives no splitting. 
Similarly for the a—a pairs. Further, we have 
the exchange terms due to like-particle, a—s 
pairs which amount to —3(x*+~y*) from the 
unlike-spin pairs and —#?z*? from the like-spin 
pairs, all involved in the same “‘internal’’ space 


integral of an exchange nature, again centrally 


symmetric, so not contributing to the splitting. 
Unlike-particle pairs have no exchange terms for 
an Hf’ independent of isotopic spin. There is thus 
no first-order doublet splitting in the alpha- 
model of Li’. 

The second-order splitting arises from admix- 
ture of excited P and F states to the two levels 
of the low ?P. (Excited S and D states are not 
admixed because the space factor of //’ trans- 
forms as the product of two vectors.) The mag- 
nitude of the matrix elements responsible for the 
admixture would depend in great detail upon the 
form of f(r) and of the alpha-model. We shall 
therefore make only a very rough order-of- 
magnitude estimate of the splitting by com- 
parison of this problem with that of the deuteron 
quadrupole moment. 

By numerical integration of an f(r) based on 
cut-off of the neutral meson theory, Bethe has 
shown that the 7-percent admixture of D wave 
function which he obtains in the ground state of 
the deuteron is quite adequate to explain the 
quadrupole moment.'* While dependent on a 
special form of the interaction, his result indi- 
cates roughly the amount of D admixture 
required by the quadrupole moment with any 
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interaction giving about the same ‘‘size of the 
deuteron,”’ except for the possibility'® of a con- 
siderably larger admixture which requires the 
small quadrupole moment to be a difference of the 
effects of the S—D and D—D quadrupole matrix 
elements. For interactions of which H’ is a small 
perturbation this latter possibility seems to be 
ruled out by comparison of the stability of the 
deuteron and alpha-particle,® together with new 
measurements of the range of the proton-neutron 
interaction.!® So it seems that the matrix element 
responsible for the D admixture is about 
(S| H’|D) ~0.07(Ep—Es) in order of magnitude, 
and this might be about 1mc*. This gives us an 
idea of the magnitude of the corresponding 
matrix elements in the alpha-model of Li’, 
between the excited 2} 4P and ‘4F states and the 
low *P states. The ‘internal’ space factors of 
these excited states have an extra node making 
them orthogonal to the ground state, so that the 
matrix elements in question would vanish for 
very short or very long range of H’, and for 
intermediate range the integrations must contain 
strong cancellation of positive and negative 
parts. This leads one to expect the matrix ele- 
ments to be much smaller than in the deuteron, 
of the order of magnitude 10~'mc?. The energy 
involved in the excitation is intermediate be- 
tween that for a triton and that for an alpha, 
perhaps about 30mc*. For the depression of one 
of the low *P levels by one excited F state we 
may thus estimate (F{|H’|P)?/(Er—Ep)=3 
< 10-%mc*. The excited P states depress both the 
levels *Pi;2 and ?P3/2, but the *F depresses only 
the *P32, suggesting that this type of doublet 
splitting would also be apt to make the J,;'=3, 
as observed. Since it depends on the detailed 
difference of the depressions of the two low ?P 
levels due to several excited states, the splitting 
might be either larger or smaller than the de- 
pression (~ 3X 10-*mc*) due to one excited state. 
We may at least make the very rough estimate 
that the second-order splitting due to H’ is of 
about the same order of magnitude as, and 
might be considerably smaller than, the Larmor- 
Thomas splitting (~3X10-*mc*) in the alpha- 
model of Li’. 


1% D. R. Inglis, Phys. Rev. 55, 988 (1939). 
16 F. O. Salant, R. B. Roberts and P. Wang, Phys. Rev. 
55, 984 (1939). 
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The Scattering of Alpha-Particles by Nitrogen 


GorRDON BRUBAKER* 
Sloane Physics Laboratory, Yale University, New Haven, Connecticut 


(Received October 13, 1939) 


A study has been made of the scattering of alpha-particles by nitrogen into four angular 
ranges of mean angles 53°, 66°, 88°, and 104°. Anomalies in the scattering are observed for 
energies as low as 3.5 Mev. Curves showing the ratio of observed to Coulombian scattering as 
a function of energy show several maxima and minima. The form of the curves varies greatly 


with the scattering angle. 





EXT to those elements which may be 

obtained in the form of thin foils, elements 
which exist in gaseous form are the most con- 
venient for use in investigations of scattering. 
Nitrogen is such an element, and one of the 
lightest. Although, presumably, the scattering of 
alpha-particles by nitrogen nuclei would show 
marked anomalies, it was not investigated during 
the early work on alpha-particle scattering since 
the scintillation screen then used as a detector of 
the scattered particles would not distinguish 
sharply between scattered alpha-particles and 
the protons produced by the transmutation of 
nitrogen. The development in recent years of 
electrical detectors which will distinguish between 
the two kinds of particles has made it possible 
to study the scattering of alpha-particles by 
nitrogen nuclei. Such studies have been carried 
out by Devons,! who has investigated the scat- 
tering at 90°, and by the author. 

In the present work the scattering at four 
different mean angles has been studied. The 
apparatus used has been previously described.’ 
The detectors used were toroidal proportional 
counters which detected particles scattered from 
a small volume on the axis of the toroid. The 
counters were filled with the gas whose scattering 
was being studied. Proper adjustment of the 
voltage applied to the counters and the ampli- 
fication of the impulses produced by them made 
it possible to detect alpha-particles but not 
protons. 

The use of two of the toroidal counters made 
it possible simultaneously to observe particles 
scattered into two different angular ranges. By 

* Sterling Fellow. 


1S. Devons, Proc. Roy. Soc. A172, 127 (1939). 
? Gordon Brubaker, Phys. Rev. 54, 1011 (1938). 


changing the positions of the edges defining 
incident and scattered beams of particles, scat- 
tering into two other angular ranges could be 
studied. The scattering of alpha-particles having 
energies between 3.5 Mev and 7.0 Mev into 
four somewhat overlapping angular ranges whose 
mean angles are 53°, 66°, 88°, and 104° was 
studied. 

Since its original description, the apparatus 
has been altered to improve its resolving power. 
The size of the angular ranges has been reduced 
by use of an alpha-particle source 3 mm instead 
of 9 mm in diameter, and by slightly reducing 
the angular width of the beam of scattered par- 
ticles. These changes have also reduced somewhat 
the spread in energy of incident particles in the 
scattering volume. The total spread in scattering 
angle for the 53°, 66°, 88°, and 104° ranges, 
respectively, are 35°, 36°, 41°, and 37°. A com- 
putation of the number of recorded scattered 
particles as a function of scattering angle, on the 
assumption of Coulomb forces, shows that, 
unless the scattering differs a great deal from 
Coulomb scattering, very few particles are 
recorded which have been scattered through the 
very smallest or very largest angles in each 
angular range. The effective angles of scattering, 
may, then, be given as 53°+11°, 66°+12°, 
88°+14°, and 104°+12°. 

The’spread in energy of the incident particles 
is due to the thickness of the scattering volume 
and to the straggling of the alpha-particles. The 
thickness of the scattering volume in each case 
was about 2 mm air equivalent. For an alpha- 
particle energy of 5 Mev this amounts to 0.2 Mev 
spread in energy. At the same energy the 
straggling in energy of the alpha-particles 
amounts to about 0.09 Mev. The total spread 
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Fic. 1. Experimental results for the scattering of alpha- 
particles by nitrogen. 


in energy of incident particles in the scattering 
volume was, therefore, about 0.38 Mev. 

For calibration of the apparatus, the scattering 
of alpha-particles by argon was used. The num- 
ber of particles scattered by argon at energies 
near 4 Mev (for which the scattering is known 
to be classical) was determined and, by use of 
the Rutherford scattering formula, the number 
which would be scattered by nitrogen on the 
assumption of Coulomb forces was calculated. 

To make sure that no protons were being 
detected, a foil thick enough to stop the scat- 
tered alpha-particles but not the protons was 
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placed between the scattering volume and 
detectors and alpha-particles of various energies 
were allowed to pass through the scattering 
volume. It was found that no protons were 
counted. 

RESULTS 


The results which were obtained are shown in 
the curves of Fig. 1, in which the ratio of ob- 
served to Coulomb scattering is plotted against 
the energy of the incident alpha-particles. From 
200 to 1500 particles were counted for each of 
the experimental points. The curves show that 
anomalies in the scattering exist for energies as 
low as 3.5 Mev, as is to be expected since the 
penetrability of the nuclear potential barrier for 
this energy is at least 10 percent, and that the 
course of the anomalies with increasing energy 
of the particles is very irregular. Particularly 
noteworthy is the high and narrow peak at 5.3 
Mev in the curve for 88° mean scattering angle, 
the appearance of which strongly suggests that 
it is due to resonance. The whole curve is similar 
to that obtained by Devons! at 90° mean angle 
of scattering. The second smaller resonance peak 
at about 4.6 Mev shown in Devons’ curve does 
not appear here, probably because of lack of 
experimental points in that region. The difference 
in height of the higher peak as shown in Devons’ 
curve and in Fig. 1 may be caused by differences 
in resolving power of the two sets of apparatus. 

A somewhat broader peak occurs in the curve 
for mean scattering angle 104° with its maximum 
at about 5.0 Mev. It is possible that this is com- 
posed of two peaks and that further and more 
accurate data would separate the two. The 
subsequent rise in the curve at higher energies 
is similar to that found in the scattering by 
oxygen at these same angles,” and to those found 
by Riezler*® in the scattering at high angles by 
beryllium, boron, and carbon. 

The curves for the scattering into the two 
lower angular ranges resemble each other, but 
bear little resemblance to the curves for the two 
higher angular ranges. The first anomaly evident 
from the data is a depression of the scattering 
ratio below unity having its minimum at about 
5.2 Mev. Following this there are two maxima 
in each of the curves. More data are needed to 


3W. Riezler, Ann. d. Physik 23, 198 (1935). 
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determine the exact positions and heights of 
these maxima. It seems probable that the depres- 
sion at 5.2 Mev is produced by the same 
resonance level as is responsible for the high 
peak in the curve for angle 88°. The difference 
(0.1 Mev) in the positions of the minimum and 
the maximum is felt to be real and not the result 
of experimental error. Such an error could arise 
only through errors in the measurement of the 
dimensions of the apparatus. This measurement 
was carefully repeated and found to be correct. 


DISCUSSION 


Bethe‘ has given a formula based on the many- 
body theory of scattering which is valid in the 
neighborhood of resonance for any spin of scat- 
tering nucleus and incident particle but for only 
zero orbital angular momentum of the particle. 
For the particular case of nitrogen this formula is 


a(@) p?>+2psin ¢+2px cos ¢ 
eae 
oo 1+x? 





’ 


where o(@)/ao is the ratio of actual to Coulomb 
scattering, {=a log sin? 6/2, a=e’2zZ/hv, @ is the 
angle of scattering in the center of gravity 
system, 


p=(2/a)(''p,/T,) sin? 0/2, x=2(E-—E,)/T,, 


and I'’p,/T, is the ratio of the partial width of 
the resonance level corresponding to scattering 
to the total level width. 

Consideration of this formula shows that, a 
given resonance energy having been assumed, 
the behavior of o(@)/ao as the energy is varied 
slightly from the resonance value is as follows: 


4H. A. Bethe, Rev. Mod. Phys. 9, 176 (1937). 





For angles @ greater than a certain high value 
6,, the numerator of the fraction is positive and 
a(@)/oo9 decreases as E either increases or de- 
creases: that is, the evidence of resonance is a 
peak in the curve of o(@)/oo against E. 

For another range of angles #.<@<6, the 
numerator is negative and the evidence of 
resonance is a depression in the curve. 

This alteration in the character of the curve 
continues as the angle is further decreased. 

In the curves of Fig. 1 there is apparent 
evidence of resonance in the high peak of the 
curve for angle 88°, and further evidence of the 
same resonance level in the depression at approxi- 
mately the same energy in the curves for the two 
lower angles. This change from a peak to a 
depression agrees qualitatively with the behavior 
predicted by the theory. Quantitatively, how- 
ever, the formula cited above predicts only 
depressions at an assumed resonance energy of 
5.3 Mev for all the angles concerned. This indi- 
cates that incident particles of zero angular 
momentum are not involved in the resonance 
process. This is not in disagreement with the 
conclusion of Devons that particles of angular 
momentum two are involved. 

The irregularities at the higher energies in the 
curves for the two lower angular ranges are not 
well defined by the present data. For this reason 
and because of the variability in the character 
of evidence for resonance discussed above, no 
assignment of the energies of the resonance 
level or levels presumably responsible for the 
irregularities has been made. 

The author wishes to express his gratitude to 
Professor A. F. Kovarik and Professor E. Pollard 
for helpful discussions during the course of this 
work. 
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The phenomenon of double @-disintegration is one for which there is a marked difference 
between the results of Majorana’s symmetrical theory of the neutrino and those of the original 
Dirac-Fermi theory. In the older theory double 8-disintegration involves the emission of four 
particles, two electrons (or positrons) and two antineutrinos (or neutrinos), and the prob- 
ability of disintegration is extremely small. In the Majorana theory only two particles—the 
electrons or positrons—have to be emitted, and the transition probability is much larger. 
Approximate values of this probability are calculated on the Majorana theory for the various 
Fermi and Konopinski-Uhlenbeck expressions for the interaction energy. The selection rules 
are derived, and are found in all cases to allow transitions with Ai = +1,0. The results obtained 
with the Majorana theory indicate that it is not at all certain that double 8-disintegration can 
never be observed. Indeed, if in this theory the interaction expression were of Konopinski- 
Uhlenbeck type this process would be quite likely to have a bearing on the abundances of 
isotopes and on the occurrence of observed long-lived radioactivities. If it is of Fermi type 


this could be so only if the mass difference were fairly large (e=20, AM7=0.01 unit). 





I. INTRODUCTION 


HE probability of double §8-disintegration 
was calculated some years ago by Goeppert- 
Mayer' on the basis of the Fermi theory.?: * The 
result obtained was extremely small, correspond- 
ing to a lifetime of the order of 105 years in the 
case of two isobars whose masses differ by 0.002 
mass unit and whose atomic numbers differ by 
two units. Thus one can account for the large 
number of isobaric pairs with AZ=2, as com- 
pared to the scarcity of isobars with AZ=1. 
Although not strictly stable, the heavier isobar 
of a pair with AZ=2 may be supposed to be 
metastable, having a lifetime large compared to 
geologic time. 

An inspection of the calculations shows that 
the results would not be changed by any sig- 
nificant factor by the use of an expression for the 
interaction energy involving derivatives of the 
neutrino wave function, as suggested by Kono- 
pinski and Uhlenbeck.‘ The same is true as 
regards the generalizations in structure of this 
expression,® which make it possible to obtain 


1M. Goeppert-Mayer, Phys. Rev. 48, 512 (1935). 

2 E. Fermi, Zeits. f. Physik 88, 161 (1934). 

5 For a review of the theory and its various modifications 
and applications up to the beginning of 1936, see pages 186- 
206 of the article by Bethe and Bacher, Rev. Mod. Phys. 
8 (1936). 

*E. J. Konopinski and G. E. Uhlenbeck, Phys. Rev. 48, 
7 (1935). 

5 Cf. reference 3, pp. 190-192. 


selection rules* for ordinary 8-decay decidedly 
different from those originally given by Fermi. 
The original Fermi picture of the fundamental 
interaction processes concerned in 8-decay has 
now been generally supplanted by a picture in 
which mesons play the role of mediaries between 
the heavy particles and the electrons and 
neutrinos. In this case also, as is evident from 
a consideration of the arguments of Yukawa,’ 
the results of Goeppert-Mayer will remain un- 
changed. 

The situation is, however, decidedly altered if 
one admits a change in the theory of the neutrino 
itself as an elementary particle. Such a change 
was suggested by Majorana® in a paper on the 
symmetry properties of the Dirac theory. 
Majorana’s suggestions have been more generally 
developed in the case of the positron theory by 
Kramers,® and for the neutral particle by the 
writer.’° Racah" has also discussed the applica- 
tion to the neutrino theory of 8-decay. 

The essential difference between the Majorana 
theory of the neutrino and the usual form of the 
Dirac theory is that in the former there are only 


6G. Gamow and E. Teller, Phys. Rev. 49, 895 (1936). 

7H. Yukawa, Proc. Phys.-Math. Soc. Japan, 17, 55-56 
(1935); H. Yukawa, S. Sakata, M. Kobayasa and M. 
Taketani, Proc. Phys.-Math. Soc. Japan 20, 731-734 (1938). 

8 E. Majorana, Nuovo Cimento 14, 171 (1937). 

9H. A. Kramers, Proc. Amsterdam Akad. 40, 814 (1937). 

10 W. H. Furry, Phys. Rev. 54, 56 (1938), referred toas N. 

1G, Racah, Nuovo Cimento 14, 322 (1937). 
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TRANSITION 


two states for a given momentum, corresponding 
to the two possibilities for the spin; there are no 
negative energy neutrinos, and no “holes,” or 
antineutrinos. If in applying the usual Dirac- 
Fermi theory to 6-decay one assumes that the 
emission of a neutrino accompanies that of a 
positron, then the emission of an electron must 
be accompanied by the emission of an anti- 
neutrino, or else the absorption of a neutrino. 
In the Majorana theory the emission of either 
an electron or a positron is accompanied by either 
the emission or the absorption of a neutrino. 

It can be shown that the use of the Majorana 
form of neutrino theory instead of the usual 
theory makes no difference in the case of ordi- 
nary 8-decay.” For the double 8-disintegration, 
however, there is a marked qualitative difference 
between the results of the two theories. In the 
ordinary form of the theory four particles must 
be emitted in such a process: two neutrinos (or 
antineutrinos) must accompany the emission of 
two positrons (or electrons). In the Majorana 
theory there can occur not only these four- 
particle disintegrations, but also disintegrations 


in which only the two charged particles—elec- 


trons or positrons—are emitted, unaccompanied 
by neutrinos. In these two-particle disintegra- 
tions the neutrino plays only a transitory or 
virtual part, such as is played by electron- 
positron pairs in certain hypothetical radiative 
processes.'* Subject to the usual limitations on 
the meaning of such language, one can say that a 
(virtual) neutrino is emitted together with one 
of the electrons (or positrons), and reabsorbed 
when the other electron (or positron) is emitted. 

In calculating the probability amplitude for 
such a process, one must integrate over the 
momentum space of the neutrino. This intro- 
duces a quantity which is a sort of ‘‘form factor” 
of the nucleus. In most cases the integral con- 
verges, giving a finite ‘‘form factor’ without any 
arbitrary “cutting off” procedure; one has to 
resort to “cutting off’’ only for certain types of 
Konopinski-Uhlenbeck expressions for the inter- 
action energy. The results obtained for the dis- 
integration probability are, nevertheless, greater 


12 Reference 10, pp. 66-67; but see limitation stated in 
footnote 22 of reference 10. 

13 Cf. for example M. Delbriick, Zeits. f. Physik 84, 144 
(1933); N. Kemmer, Helv. Phys. Acta, 10, 112 (1937); 
H. Euler, Ann. d. Physik 26, 398 (1936). 
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than those for the four-particle process by a 
factor which ranges from 10° to 10° or more, 
depending on the particular interaction ex- 
pression used. The fact that the neutrino 
momentum is not limited by the amount of 
energy available in the transition also makes the 
probability less strongly dependent on the 
energy. 

In the following section we carry through the 
calculation for the simplest case, using the 
“scalar-type’’ interaction energy expression" and 
assuming that the nuclear spin does not change 
and that the transition is of even-even or odd-odd 
type. In Section III we consider the effects of 
using various other expressions for the inter- 
action energy, and determine the selection rules 
governing changes in the nuclear spin 7 and the 
even-odd or even-even, odd-odd character of 
the transitions. 

Before we continue beyond the introduction it 
should be-remarked that the Majorana form of 
the theory is not the only one which permits 
this new type of disintegration. The same sort 
of result can be obtained with the more usual 
theory if one introduces an interaction expression 
involving a linear combination of the neutrino 
wave function and its conjugate. Such an intro- 
duction has never been suggested, however, in 
the neutrino-antineutrino type of theory. The 
Majorana theory is the natural one to use for 
our purpose, both because, in the elimination of 
the antineutrino, it provides an independent 
motive for the use of such expressions, and 
because it provides, so to speak, a canonical 
form for them. 


II. ScALAR-TyPE INTERACTION WITHOUT 
DERIVATIVES; At=0, EVEN-EVEN, 
Opp-Opp TRANSITIONS 


We consider the transition of two neutrons in 
the nucleus into two protons, with the emission 
of electrons into states having wave functions 
y, and y¥; and energies H, and H;. The formula 
for the probability amplitude must contain two 
terms, since either electron may be emitted first. 
These terms are of opposite signs because of the 
antisymmetry of the electronic wave functions; 
whether or not they actually cancel strongly 
depends on the interaction expression used. 


4 Cf. Eq. (208a) of reference 3. 
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The neutrino wave functions are taken to be 
expanded in terms of plane waves subject to a 
periodic boundary condition in a box of volume 
V, as expressed in Eqs. (30) and (31) of N. Only 
the neutrino amplitudes are regarded as quan- 
tized, and all other wave functions are treated 
as c-numbers. 

Apart from a resonance factor, the square of 
whose absolute value may be set equal to 
(2nt/h)6(Ev+H,+H,—Ewm), the probability am- 
plitude is 


ane =e DUE f vs" {SB.0*(ri|k)} Vidr 


f 1! OB,0,t(4/' | —k)} Va'd" 
2 





{ (y.'Ba(k))(y."8a(—k)) 
Eu—E.—H,—Ky. 
—ee] 1) 


Ey—E,—T1,-—K;, 





Here g is Fermi’s constant; K, is the neutrino 
energy in the virtual intermediate state; and 
vu, wv, and Wy, are the wave functions of the 
nucleus in the initial, final, and intermediate 
states. They are functions not only of the space 
coordinates of all the heavy particles, but also 
of the four-valued Dirac spin coordinates on 
which matrices such as 8 operate, and of a two- 
valued fifth coordinate for each particle which 
indicates whether it is a neutron or a proton. 
The indicated integrations must be understood 
as including summation over the fourth and 
fifth coordinates. It is on the fifth coordinates 
that the Fermi operators Q;* act in such a way 
as to change a neutron into a proton; if the ith 
particle is already a proton, the operator Q;' 
gives the result zero. The operators 6,Q;' are 
symmetrized by summing over all the particles 
in the nucleus, and the factor (r;|k) appears 
because the neutrino wave function must be 
given its value at the position of the heavy 
particle. The electron wave functions y, and ¥; 
are to be evaluated at the “outer edge of the 
nucleus,”’ as is usual in calculations on 6-decay.?:? 

According to the Majorana theory, the four- 
component neutrino amplitudes a(k) are given 
by Eqs. (36), (50), and (53) of N. We have to 
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substitute from these equations into (1), and 
then put 
B,B;=0, B*B;=0 (2) 
and 
B(k)B;*(k)=6;;,, B*B*=0. (3) 


Eq. (2) corresponds to the initial condition that 
no neutrinos are present, and (3) indicates our 
intention to calculate the process in which no 
neutrinos are emitted and not that in which two 
neutrinos are emitted. As indicated in N, the 
calculation can be simplified by using Eq. (51) 
of N instead of (50), and then discarding all 
terms in b3, 54, b3*, and }4*. Our conditions (2) 
and (3) can be satisfied by using instead of 
Eq. (51) of N just 


b(k)—b(k) (4) 
for the left-hand factor and 
b(k)—>T™* (k)b*(—k) (4’) 


for the right-hand factor. Then, remembering 
that in the Majorana form of 8-decay theory it 
is consistent to regard the factors (2)~! given by 
Eq. (53) of N as absorbed in the disposable 
constant g, we obtain 


(v.'Ba(k))(y."Ba(—k)) 
= L (WBS(k))(Y.18S(—k)T*(—k)); 


i=1, 2 


= Wit BS(k) -2(1+ ps): T1(—k)S'(—k)p'y.*. (5) 


Here p3 is a matrix originally defined by Dirac,'® 
and quoted in N, Eq. (39). Then we can use 
the value of 7(k) as given in N, Eq. (43), and 
obtain 


(v.'Ba(k)) (y."Ba(—k)) = ¥.'8S(k) 
2(1+ps)St(kK)AtBY *=yptMy.*. (6) 


It is readily verified that (6) has the same 
invariance properties as Eq. (56) of N with 
respect to changes in the form of the matrices 
a;, 8, and S(k), and that if the neutrino mass is 
equal to zero M is a symmetrical matrix, inde- 
pendently of such changes. We shall use the 
representation of a;, 8 which was originally given 
by Dirac, and has usually been employed in 
explicit calculations; it is quoted in N, Eggs. 
(59), (62). We can accordingly substitute Eqs. 


1 PA. M. Dirac, Proc. Roy. Soc. A117, 614 (1928). 
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(60)-(62) of N into (6); we shall, however, 
replace both ¢ and e+1 by k, which corresponds 
to setting the neutrino mass equal to zero. 
The result is 


M=—}k"{k+(a-k)}-ia,B=M’. (7) 


It is certainly correct to set the neutrino mass 
equal to zero, since important contributions to 
the result will come only from large values of 
the neutrino momentum. This step leads to a 
formal difficulty, however, because in N the 
vector k represents the propagation vector in 
units mc/h, where m is the neutrino mass. We 
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must, accordingly, change the units of k; this 
makes no difficulty, because (7) is homogeneous 
of degree zero in k. We shall from now on take k 
to be expressed in c.g.s. units, and replace Eq. 
(31) of N by 


(r|k) = V~- exp (ck-r). (8) 


When (7) is substituted into (1), we see that 
on account of the symmetry of the matrix WM the 
two numerators are equal. The minus sign due 
to the antisymmetry of the electronic wave 
function thus causes a strong cancellation, and 
we get 


Qvem=2?°VD (k)D f Vy { OB,0¥t exp (ik-1;)} Widr- f W1't{ 58,0; exp (—ik-1;')} y'dr 
L i 7 
-(H,—H,)(Em—E.—H,—K,)-“(Eu—E1.—Hi—Ki)-Wi My,*. (9) 


The main difficulty in the calculation occurs in the evaluation of the factor 


raed f vw Z6.00 exp (ik-r,)} vidr- fw" ZB,0h exp (—ik-r/)} Wa/dr' -f(Ex) (10) 
L ; j 


where f(Ex) = (Ey —E,—H,—K,)-“(Eu —E,.—H.-—K,). 


For values of the neutrino energy small enough 
to allow us to put exp (¢k-r;)=1, the first 
factor of the summand is 


J vat 8:0.) Wade. (11) 


This is just the sort of factor which occurs in 
calculations on single §8-disintegration. It is 
usually assumed to have the value 1 if the angular 
momentum selection rule is obeyed, and, of 
course, the value 0 otherwise. Actually this is a 
decidedly crude assumption, as has been pointed 
out by Nordheim and Yost;!® for moderately 
heavy nuclei it seems likely that the expression 
(11) cannot have a value larger than perhaps 0.1. 
Since, however, the value of g given by Fermi 
was assigned to agree with experiment on the 
basis of setting an expression such as (11) equal 
to 1, we shall not make any very large net error 
if we use this value of g and assume that there is 
at least one value of L for which both of the first 
two factors of the summand in (10) can be set 
equal to 1. There still remains, however, the 

16. W. Nordheim and F. L. Yost, Phys. Rev. 51, 942 
(1937). Very recently (Phys. Rev. 56, 519 (1939)) the 


matrix elements for a number of light radioactive elements 
have been calculated explicitly by Wigner. 





question as to how many terms in the summation 
will yield similarly appreciable contributions to 
the value of I. 

At first one might suppose that there would be 
a great many such values of L, because of the 
large number of energy levels possesed by any 
moderately heavy nucleus and because k is not 
restricted to small values, so that limitations on 
the angular momentum of the state Z do not 
appear in the usual way. Such a conclusion, 
however, would in all probability be erroneous. 
The nature of nuclear states is such that any 
appreciable excitation energy is practically cer- 
tain to be distributed effectively among a con- 
siderable number of particles, so that a factor 
exp (tk-r;), which involves the coordinates of one 
particle only, is unlikely to be very effective in 
removing the orthogonality of two wave func- 
tions. The assumption that there is only one 
intermediate state L which makes an important 
contribution to the sum is accordingly not an 
unreasonable one; we shall proceed in this way, 
with the reservation that our result is perhaps 
best regarded as a lower limit, from which the 
correct result should not differ by any very 
large factor. 
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In accordance with the above considerations, we shall use in our further calculations the simple 


l= g?f(E1)(4rp*, 3*f f fexp (ie-n)dr- ff f exp (—ik-r’)dr’. 
r<p r’<p 


Here the integral is over the ‘interior of a three-dimensional sphere of radius p, the nuclear radius. 
‘ The factor (42p*/3)-* comes from normalization. The formula has been written on the assumption 
that the angular momentum of the state L is the same as that of M and N, but this assumption 
has actually rather little to do with our final result. If, for example, i,—iy=i,—iv=+1, we 
should have to use some such expression as 


I’ =g°f(E1)(4rp' 3-3-f ff cos 6-exp (ik-nydr- ff f cos 6’ exp (—ik-r’)dr’, =(12’) 
r<p r’<p 


and we shall verify later that this makes no decided difference. 

Further simplification in Eq. (9) can easily be made. Since & is in c.g.s. units, the energy of the 
neutrino is K,=/ck. The main contributions to (9) come from values of k so large that the other 
energies can be neglected compared to K,; the occurrence of a vanishing denominator is of 
course ruled out by the postulated impossibility of single 6-disintegration for the nucleus in 


question. Accordingly we set 


formula 
(12) 


f(Ex) = Ki? = (he) *k. (13) 


The sum over the allowed values of k in the box of volume V can be replaced by an integral in the 


usual way: 


E(k): = (nav f mak [ f ao--. (14) 
0 
Substituting (12), (13) and (14) in (9), we obtain 
dyeu = — (2r)-9(g/hc)?C??- (H,— Hi) Wi - dia, By,*, (15) 
where 
| 2 
; (16) 











cum fae f fan 
0 


the superscripts 0,0 refer to the values of 7,—i 
and i,—tiy. The second term in (7) is omitted 
from (15), since it makes no contribution for 
even-even or odd-odd transitions. 


(3 eff for (ik-r)dr "=(9 avon fal f f foo 


approximation, valid because | S| <1: 
R, 2(162/|P(3-+25S) | 2)(me/h)*(2mep/h)2s 
-Qryh./p.. (17') 





Following Goeppert-Mayer, we express the 
energies of the electrons and the total available 
energy in terms of dimensionless quantities: 
H,=hme, H,=himec, Eu—Ey=emc, where m 
is the mass of an electron. Let us also put 
pe=(h.2—1)', p:=(h?—1)', and make the abbre- 
viation 
R= (16m/ | T'(3 +25) |*)(mc/h)*(2mcpp./h)’s 

-exp (wyh,/p.) |T(1+S+ivh./p.) |, (17) 


where y=Z/137, S=(1—vy*)!—1. We follow 
Goeppert-Mayer in replacing R, by the following 


A similar definition and approximation hold 
for R:. 

When the resonance factor is included, the 
squared amplitude for the transition to two 
particular states s and ¢ is (2at/h)|awew|? 
5(emc?—h.mc?—hmc?). The total probability 
per unit time for transitions in which one electron 
gets an energy between h,mc*? and (h,+dh,)mc 
and the other receives the rest of the energy 
emc’ is then 
P(h,)dh,=(22/hmc?)(dh,) > > | avem|?. 


dhe dhs 





(18) 








_- -_, wv \v 
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Here the symbol >> means the sum over all states 
dht 


having energy between hymc? and (hi+dh,)me’; 
the result is of course proportional to dh,. Using 
the Dirac wave functions for an electron in the 
Coulomb field, evaluated at r=p, we obtain, to 
the lowest order in y’: 


pm » | vit , >1a,By.* 5 
dhs dht 
=1R.Rippihshi—1)dhdh. (19) 
Then 
P(h,)dh,=2-*x-*g*h-5mc—* | C®.°| ? 
XR: Rip pilhshi- 1)(h,—h,)*dh,. (20) 





The energy spectrum of the electrons is, of 
course, symmetrical around the value }emce’. 
The factor (4,—h;)? in (20), which makes the 
intensity vanish for the energy }emc*, comes 
directly from the cancellation due to the anti- 
symmetry of the electronic wave function. 

To get the total probability of disintegration 
per unit time, we must integrate (20) from 
h,=1 to h,=e—1. According to the approxi- 
mation (17’) the quantity R,p,/h, is essentially 
independent of h,. This enables us to write for 
the total probability of disintegration per unit 
time: 


e~} e—1 
Pi=f Ph )dhy=2-*e-gth-Imc*| C%8|(Rypa/he)® f h.(e—h,) 
1 1 


Then 


X {hs(e—h,) —1}(€—2h,)*dh,. (21) 


P,=2-"x—5g*h-mc—| C%|*(Rip./hs)?- o(e—2), (22) 


where 


v(x) = (x*/3) {1+ (11x%/10) + (x?/5)+(x*/70)}. (23) 


We shall now evaluate P; numerically for Z=31, the value used as an example by Goeppert- 
Mayer. Taking p=8X10-" cm, we get from (17’) the value R,p,/h,=1.52X10* cm-*. The 
integral (16) can be evaluated by a rather tedious elementary calculation, or more elegantly 


by means of Fourier’s integral theorem: 


foal ff femas 





=f ae: f [edsaxdy- f f fewvaran'ay'=nf ff f faxdyas'ay’as, 
— r<p r’<p 


where in the last expression the limits for x’ and y’ depend on the value of z in the same way as 


those for x and y. Then 
fra 
0 


C= 1242/5p. (24) 


A similar calculation shows that if (12’) were 
used instead of (12) we should have simply 
to replace C®.° by a quantity C':!=92?/5p 
= (3/4)C® °. Using g=4X10-* erg cm’, we ob- 
tain finally 
P,=1.07 X10-**g(e—2) sec. 
=3.4X10-*9(e—2) year. (25) 





Thus 


Some values of ¢ are as follows: 


e= 3 4 6 8 10 12 20 


g(e—2) =0.77 22 8.110? 7.7 X10® 4.1X10# 1.54105 5.9 x 108 (26) 





2 p 
SS fever =-rf w?(p? —2*)*dz= 167% p*/15. 
r<p ~ 


For small values of « the P found here is about 
10° times that given by Goeppert-Mayer. For 
larger values of «€ the ratio is not so large. This 
is natural, since for large « there is plenty of 
energy available for the emission of neutrinos 
along with the electrons. 


III. OrHER INTERACTION EXPRESSIONS 


In the preceding section we calculated the 
value of P by using the scalar-type expression 
for the interaction energy. We now want to 
consider what would be the effect of using any 
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of the various other energy expressions which 
have been proposed. 

Since we are interested only in the order of 
magnitude of P, we shall neglect the effect of 
using a different matrix instead of }ia,6 in calcu- 
lating the expression corresponding to (19). This 
would at most replace the factor (h,4,—1) by 
some factor such as hh; or (h,t,+1), and for 
¢— 21 such a change amounts to only a factor 
of perhaps two in the answer. Hence we shall 
simply keep the factor (h,4,—1) throughout. 

To begin with we shall consider the two main 
differences in the size of P which can arise from 
changes in the interaction energy expression. 
We shall then list the various expressions and 
indicate the main types of terms which they give 
in the integrand of the equation corresponding 
to Eq. (1) in the general case. This classification 
makes it possible to state the selection rules for 
the different types of interaction expression. 


Differences in order of magnitude of P 


The first important difference which may come 
from the use of a different interaction expression 
is the absence of the strong cancellation due to 
the antisymmetry of the electronic wave function 
which occurred in Section II. If the matrix!’ M is 
not purely symmetric, the main contribution to 
ay+m will come from the antisymmetric part of 
M, Ma=}(M—M’). Eq. (15) is then replaced by 


dvem = (2r)~*(g?/hc)D° "Yt May.* (27) 
and Eq. (20) by 
P(h,)dh,=2-*2-5g*h-3m-'c-4| D992 
XR:Ripsbil(hshi—1)dh,, (28) 


where 


prv= f rade f fav 
x\(3 aro’) { { [exp (ike-n)dr| . (29) 


17 Tn general it is actually a matrix function which has to 
be used, as is indicated below in Eq. (45). In the case 
treated in Section II this function could conveniently be 
split into factors, only one of which is a matrix (Cf. Eq. 
(46)). In any case it is only the symmetry or antisymmetry 
under the interchange of two four-valued indices which 
comes into question in determining the presence or absence 
of the cancellation. 
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The spectrum now has a maximum at h,=}e, 
instead of vanishing there. 

The total probability of disintegration per 
unit time is now 


P_=2-'2-5g!h-§m-'c-4| D9.) 2 
X(Rips/hs)?-x(e—2), (30) 


where 


e—1 
x(e-2)= f h.(e—h,) {hs(e—hs) —1}dh, 
1 


or 


x(x) =x?{1+(7x/6)+(x?/3) +(x /30)}. (31) 


Integration of (29) by elementary methods gives 


D®.%= 187/ p?. (32) 
Then 
P.=1.43 X10~*4x(€—2) sec.~! 
(33) 
=4.6X10~-'"x(e—2) year—!. 
«=3 4 6 s 10 12 Dt 


x(e—2) =2.5 19.7 2.110? 9.8 X10? 3.1108 7.9108 1.05 x 105 


Another important effect on the disintegration 
probability is obtained if we use a Konopinski- 
Uhlenbeck form of interaction expression. The 
use of such an expression is subject to certain 
theoretical difficulties,!* and it has been sug- 
gested” that the experimental results can perhaps 
be fully explained without using an interaction 
expression involving derivatives; nevertheless it 
is not without interest to see what such ex- 
pressions give for the double 6-disintegration. 

The use of a Konopinski-Uhlenbeck expression 
means essentially the insertion of a factor k? into 
the integrand in (16) or (29), and the replacement 
of the constant g by a constant which differs 
from it by the dimensions of a length. The length 
by which g is multiplied to give the new constant 
should be (27) times a typical de Broglie 
wave-length of the neutrino in an ordinary 
8-disintegration, in order that the probabilities 
for single 8-disintegration may be kept the same. 
Hence we shall replace g by (4/2mc)g. The total 
disintegration probability per unit time is then 


P3=2-"a-5g4h-lm-3¢-8| F0.0|2 
X(R.ps/hs)*- e(e—2), 
18 Cf. M. Fierz, Helv. Phys. Acta 10, 123, 284 (1937). 


19H. A. Bethe, F. Hoyle and R. Peierls, Nature 143, 200 
(1939). 


(35) 











TRANSITION 


where 


Poon f eak { [ao 
0 
x|(s trot) { [ fexp (tk-r)dr 
r<p 


if M is a pure symmetric matrix. If M is asym- 
metric we get 


(36) 





Py=2-x-Sy4hm-Fc-8 | G0.0| 2 


X(Rsps/hs)?x(e—2), (37) 


where 


Grom f 2erdk f fan 
0 
x|(3 tro’) f [ fexp (ik-r)dr 


The integral (36) can be evaluated elegantly 
by using Fourier’s integral theorem: 


F.0= (3 trp)? [ ak. f ak, | dk, 
ff fesp (tk: r)dxdydz 
r<p 
ff fexe (—1tk-r’)dx'dy'dz’ 
r’<p 


= (3/4mp*)?-(2x)?- f f f dxdydz = 6r?/ p?, 
r<p 


(39) 


(38) 





Numerical computation then gives: 


P3=2.28X10-*y(e—2) sec.— 
=7.17X10-y(e—2) year. (40) 


The integral (38) diverges logarithmically, so 
that in order to obtain a finite value it is neces- 
sary to “cut off” at some value of k. It is generally 
assumed that the.value at which the theory 
becomes unreliable and the cutting off should be 
done is about mc*/e®. In evaluating G®° it is 
particularly convenient to cut off at a value 
which makes kp=7. One then readily finds by 





PROBABILITIES 1191 


elementary methods that 


G°%= (727 09) f sin? y-dy/y—} ° 
0 


=162/p*. (41) 


The cut-off occurs rather beyond the first 
maximum of the integrand. The range 0<k < 2/p 
includes the essential contributions to the 
integrals CC®, D®°, and F®*, so that their 
values are genuinely independent of whether one 
cuts off or not. 

Using the value (41), one obtains 


P,=3.99X10-*x(e—2) sec. 
= 1.26 10-x(e—2) year-!. 


(42) 


We have now listed four typical formulas for 
the disintegration constant P. Next we must 
classify the various possible expressions for the 
interaction energy in such a way as to show 
what value of P each would give. 


Classification of interaction energy expressions 


By methods analogous to those used in de- 
riving Eq. (7) one can readily calculate the value 
of adyey corresponding to any particular ex- 
pression for the interaction energy. The formulas 
obtained are, however, for the most part very 
long and cumbersome. Accordingly we shall 
indicate only the typical characteristics of the 
dominant terms, which suffice to determine the 
order of magnitude of the disintegration constant 
and to determine the selection rules. 

Usually in calculations using the neutrino 
theory of 8-decay one neglects terms involving 
the small components of the wave functions of 
the heavy particles compared to those which 
involve only the large components. This is done 
because the small components are of order i/c 
compared to the large, 6 being a typical velocity 
of a heavy particle, and i/c is only about 0.2. 
For our present purposes we can neglect such 
terms usually, but not always; if the terms in- 
volving only the large components give purely 
symmetric matrices M and other terms give 
asymmetric matrices, then the latter will be more 
important because of the absence of any strong 
cancellation. This situation arises with Fermi’s 
original ‘‘vector-type’”’ interaction expression. 
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Usually only the product of the time components 
is used in calculations, because the product of 
the space components brings in small com- 
ponents of the heavy particle wave functions. 
In our present case, however, such a procedure 
gives a symmetric matrix M, and the space 
components appear in the actual dominant 
terms. 

To avoid useless multiplication of the number 
of types listed, we shall write the typical ex- 
pressions in terms of just the matrices 1 and ¢ for 
the heavy particles, whose wave functions may 
then be thought of as two-component Pauli wave 
functions Any occurrence of the small com- 
ponents in the dominant terms will then be 
indicated by a suitable power of (#/c). Explicitly, 
so far as the heavy particles are concerned we set 


B—1, y¥s=p1—>(8/c) 
Boo, a-—(i/c)e, Ba—-(v/c)e. (43) 
We use the abbreviations: 
Wit {LON | —k)}Vyw=a,, 
Wit{ dS e.0i(r;| —k)} Vu =bi, 
, (44) 


Wyt{ Ot (ri |k)} Vi =az, 


Vy} o0;" (1; | k)V;=be. 


The quantity ay+, can always be expressed in 
the form 


aveu=eE(k) f dr fadr’S Lvu* 


im 





Aim’(r, ¥! 3 kk) 
Am“(r, ¥’ 3k) 
e-ink, 








Hen (45) 


In the case treated in Section II the matrix 
function A was simply 


A(r, r’; k) =ae(r; k)M(k)ai(r’;k), (46) 


with M given by Eq. (7). Eq. (9) gives the result 
of substituting (46) and (44) in (45), provided 
(43) is taken into account. 

We now list the types of expression which can 
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appear as the dominant terms of the matrix 
function A: 


Symmetric Terms: 


(47a) 
(47b) 


S1=d2"ta,Bay, 


S2=d2(a-n)ia,Ba, (n=k/k). 


Many other types of symmetric term occur, but 
only in combination with antisymmetric terms 
which give the dominant contributions to the 
result. 


Antisymmetric Terms: 


A,=i([bi Xbe]-e)ia,£, (48a) 
A,=i([bi Xbe]-n)piia, 8, (48b) 
A;=i([b:Xb2]-n)(n-o)ia,B, (48c) 
A,=6((b:1Xn]-[beXn))ia,B, (48d) 
A;=i([(aibe—asb;) Xn]-e)ia,B. (48e) 


We give the different expressions for the inter- 
action energy the same designations as used by 
Bethe and Bacher.* The types of the dominant 
terms in A and the order of magnitude of the 


disintegration constant P are given in the 
following list : 
Fermi Expressions: 
Scalar : A~Si, Se; P=P,, (49a) 
Vector: A~(b/c)As; P~(8/c)*Pe, (49b) 
Tensor : A~A,,A2; P=P2o, (49c) 
Fseudovector: A~A;,A2; P=P2, (49d) 
Pseudoscalar : 
A~(i/c)*Si, (6/c)*?S2; P~(i/c)*Pi. (49e) 
Konopinski-Uhlenbeck Expressions: 
Vector: A~qSi, gS2; P=P3, (50a) 
Tensor : A~qAo,qA3; P=P,, (50b) 
Pseudovector : 
A~qA1, gA2, gA3,qAs; PP; (50c) 
Here , 
q=k?(h/2mc)* (51) 


is the essential factor which distinguishes the 
Konopinski-Uhlenbeck from the Fermi formulas, 
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and which was taken into account in calculating 
P;and P,. 


Selection rules 


The expressions 
DLviiSw.* and DytAw.* 
L L 


are all scalars with respect to the three-dimen- 
sional rotation-reflection group: so far as trans- 
formation properties are concerned we can think 
of ia,By,* as replaced by a function y,.” To find 
the selection rules corresponding to the ex- 
pressions (47) and (48) we may inspect the 
transformation properties of either the factors 
involving @,, b,, and n or the matrix factors, 
apart from the factors ia,8. We obtain: 


S; and Ay: Scalars: Ai=0; odd-odd and even- 


even transitions. 


As: Scalar, but with reflection character —1: 
Ai=0; odd-even and even-odd transitions. 


Se: Polar vector: At=0 or +1, but not 0-0; 
odd-even and even-odd transitions. 


A, A3, and A;: Axial vectors: Ai=0 or +1, but 
not 0-0; odd-odd and even-even transitions. 


The selection rules for each interaction ex- 
pression can then be read off from (49) or (50). 
The selection rules in general bear more re- 
semblance to those of Gamow and Teller than to 
those of Fermi, since some Ai= +1 transitions 
are always allowed. Transitions forbidden by 
these rules can occur only if at least one electron 
is emitted in a state with 7+}3, and their proba- 
bilities will be decidedly smaller than those 
given by (49) and (50). For e—2~2 the proba- 
bilities will be decreased by a factor of roughly 
{(€—2)pmc/2h}* for a change in parity type of 
the transition or for each extra unit of Az. For 
smaller values of «—2 the decrease will be more 
drastic. 


20 Cf. discussion in N, p. 66; also W. H. Furry, Phys. 
Rev. 51, 125 (1937), Eq. (8) with A =ita,@. The free electron 
functions form a complete set. 





Dependence of transition probability on atomic 
number 


The numerical results which we have given 
were all computed for Z=31, p=8X10-" cm. 
The quantity (R,p./h,) is roughly proportional 
to Z, and p may be taken to be proportional 
to Z'. Approximate relations for changing the 
results to different values of Z or a different 
estimate of the nuclear radius are, accordingly : 


Py x Z%p-2 « Z4/3, P. x Z%p-4 « Z3, 
P;<Z*p-*xcons., Pyx«xZ*%p-§«Z-!, (52) 


The results as calculated are for emission of 
two negative electrons. In the case of positron 
emission the probabilities will be smaller, because 
R, contains an additional factor exp (—2yh,/p.), 
and R, a similar factor. These factors decidedly 
complicate the evaluation of the integrals from 
which the functions ¢ and x have to be deter- 
mined. It is evident, however, that for «—2~2 
the decrease will be by a factor not much smaller 
than e~**7=¢-°-°%Z_ As is the case for the for- 
bidden transitions, the decrease is more decided 
for smaller values of «—2. 


IV. CONCLUSION 


We have seen that the phenomenon of double 
beta-disintegration is one for which there is a 
decided difference between the results of the 
Majorana theory and those of the older theory 
of the neutrino. According to the older theory it 
seemed certain that double beta-disintegration 
could never be capable of observation because 
of its extremely minute probability, but the 
Majorana theory indicates that this is by no 
means necessarily the case. Indeed, if the inter- 
action expression were of Konopinski-Uhlenbeck 
type this process would be quite likely to have 
a bearing on the abundances of isotopes and on 
the occurrence of observed long-lived radio- 
activities. If it is of Fermi type this could be 
so only if the mass difference were fairly large 
(eX20, AM0.01 unit). 
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Phase Shift Calculations for Proton-Proton Scattering at High Energies 
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The theoretical s wave shift for the square well and the Gauss error well which is fitted to 
experiment from 0.7 to 2.6 Mev is calculated up to energies of 10 Mev. Coulomb functions 
needed for this purpose are tabulated and the phase shift, as well as the ratio of theoretical 
scattering to that expected according to Mott’s formula, is represented graphically as a func- 
tion of the energy. The Gauss error well and the square well are found to give very similar 
extrapolations of the experimental phase shifts to high energies. 


NALYSIS':? of the experimental data* on 

proton-proton scattering has shown that 
the p state scattering is small in the energy 
interval 0.9 to 2.4 Mev. The analysis also shows 
that on the assumption of pure s state scattering 
in this energy region the experimental data are 
closely fitted by the scattering expected from a 
square potential well of depth 10.50 Mev and 
range e*/mc* without Coulomb potential inside, 
or from a Gauss error potential well Ae-*”’ 
where A =51.44 mc? and a=21.59 Mmc?/h?, cut 
off at r=3e?/mc? and with Coulomb potential 
inside. Comparisons with experimental data were 
also made? with the use of square wells of depth 
6.3452 Mev and range 1.25 e?/mc*, and of depth 
19.690 Mev and range 0.75 e/mc’. 

The probability that proton-proton scattering 
experiments will soon be done at energies from 
3 to 10 Mev has made it desirable to have the 
phase shifts and the expected scattering for the 
above wells calculated for these higher energies. 
Graphs up to 9 Mev have already been pub- 
lished? for the square wells of depth 10.50 Mev 
and radius e*/mc?, and of depth 46.78 Mev and 
radius 0.5 e?/mc?. In this paper the high energy 
calculations for the others of the above potential 
wells are reported. 

The notation and the values of the physical 
constants are the same as those used in BTE. 


* Now at the Agricultural and Technical College of 
North Carolina. 

1G. Breit, E. U. Condon and R. D. Present, Phys. Rev. 
50, 825 (1936). This paper will be referred to as BCP. 

2G. Breit, H. M. Thaxton and L. Eisenbud, Phys. Rev. 
55, 1018 (1939). This paper will be referred to as BTE. 

*R. G. Herb, D. W. Kerst, D. B. Parkinson and G. J. 
Plain, Phys. Rev. 55, 603(A) and 998 (1939); N. P. Hey- 
denberg, L. R. Hafstad and M. A. Tuve, Phys. Rev. 55, 
en (1939). For analysis of the data, see BTE pp. 

—6. 


The quantities o*/Po, oO and Cy?pPo” used in 
the calculation of the phase shift Ko (BCP Eq. 
(7.8)) are given in Table I. In the other tables 
are listed the values of the quantities necessary 
to find the ratio ® of the expected scattering to 
Mott scattering (BTE Eqs. (2.1) and (2.2), and 
see BCP Eqs. (6) to (6.7)). Tables II, [V and VI 
give the values of the coefficient of —sin K, cos K, 
in the formula for ® for L=0, 1, 2. Tables III, 
V and VII give the values of the coefficient of 
sin? K;, for L=0, 1, 2. Tables II to VII are ex- 
tensions of the BCP tables of the same numbers, 
and Tables II and III are extensions of BTE 
Tables I and II. 

The calculations with the Gauss error potential 
were made at energy intervals of 1 Mev. The 
joining to the Coulomb functions was made at 
r= 3e?/mc*. Interpolations were necessary to ob- 
tain the value of xj}’/§ at the value of x (here 
x=a'r) corresponding to r= 3e?/mc*; these inter- 
polations were made parabolically. At energies 
of 6 and 7 Mev the value of xj}’/§ changes so 
rapidly with x that interpolations of § and ’ 
separately were necessary. 

Figures 1 to 4 show the phase shifts Ko for 
the different potential wells as a function of the 
energy E of the incident protons. Figs. 1, 2 and 3 
illustrate the effects of varying the depths of the 
square wells of radius 0.75, 1.00-and 1.25 e?/mc* 
without Coulomb potential inside. Increasing the 
depth of the well mainly raises the Ko, E curve 
as a whole, and does not change its shape or 
slope much. Fig. 4 shows how the shape of the 
Ko, E curve for the square well without Coulomb 
potential inside is affected by changing the 
radius of the well. Decreasing the radius in- 
creases the slope and slightly reduces the curva- 
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TABLE I. Coulomb functions. This table is probably accurate to +0.02 percent. It is an extension of the table in the appendix 
of BTE and gives the quantities needed for Eq. (7.8) of BCP. 



































r=0.5 e2/mce r =0.75 e2?/mce r =e2?/mct 
E (MEv) n Po*/Po D080 Co2pho? Po*/Po PoGo Co%pPo? Do*/ Po Poo Co? pPo? 
3.0 | 0.091279) 1.0005 0.9228 0.2034 0.9829 0.8789 0.3029 0.9530 0.8215 0.3972 
3.5 .084508| .9965 .9189 .2239 .9739 .8678 .3319 .9367 .8003 4324 
4.0 .079050;} .9925 .9146 .2429 .9649 .8561 .3583 .9204 .7784 4635 
4.5 .074529| .9886 .9098 .2606 .9559 .8439 .3824 .9040 .7562 4912 
5.0 .070704| .9845 .9048 .2772 .9470 .8314 4047 .8875 .7338 5161 
5.5 .067414| .9806 .8996 .2929 .9377 .8186 4254 .8708 .7114 5385 
6.0 .064544| .9765 .8942 .3075 9285 8057 4446 .8540 .6890 .5587 
6.5 .062012| .9725 .8886 .3215 .9193 .1927 4625 .8372 .6667 5770 
7.0 .059756| .9685 .8829 .3348 .9101 .7795 4792 .8202 .6446 .5935 
7.5 .057730| .9645 8771 .3476 .9009 .7663 4949 .8031 .6225 .6084 
8.0 .055897| .9604 .8712 .3597 8916 .7531 .5096 .7858 .6009 .6219 
8.5 .054228| .9564 8653 3713 .8823 -7400 5234 .7685 .5793 .6340 
9.0 .052700; .9524 8593 3825 .8729 .7268 5364 .7510 5583 .6449 
9.5 .051294)| .9483 8532 .3932 .8635 .7136 .5486 .7334 .5372 .6547 
10.0 .049996 8541 .7006 .5602 
r =1.25 e2/mc r=2 e2/me r=3 e2/mc 
E (Mev) ” Do*/ Po $080 Cope? Do*/Do P80 Co®pPo? Po*/ Ho PoGo Co2pHo? 
3.0 0.9108 0.7516 0.4818 0.7004 0.4781 0.6544 0.1728 0.08249 0.6451 
35 .8851 .7183 .5194 .6282 4063 .6764 — .03225 .6027 
4.0 .8591 .6847 5515 .5538 3387 .6880 — .2555 — .07356 .5510 
4.5 .8329 .6512 5791 4771 .2755 .6914 — 5005 4942 
5.0 .8063 .6179 .6028 3980 .2169 .6883 —.7715 — .16425 4359 
Be .7795 5851 .6231 3163 .1623 .6796 — 1.0742 3780 
6.0 .7524 .5529 .6404 .2318 .1133 .6666 — 1.416 — .2045 3223 
6.5 .7250 5213 .6550 .1444 .06805 .6500 — 1.806 .2700 
7.0 .6972 4904 .6672 .05367 .02701 .6306 — 2.259 — .2080 .2219 
7.5 .6692 4601 .6773 — .04038 — .01005 .6089 — 2.795 .1784 
8.0 .6409 4306 -6855 — .1382 — .04325 5854 — 3.440 — .1867 .1399 
8.5 .6122 4020 .6919 — .2398 — .07283 .5606 —4.239 .1064 
9.0 .5831 .3740 .6967 — 3460 — .09894 5349 — 5.262 — .1498 .07794 
9.5 | 5537 .3467 .7001 | —.4568 — .12187 .5086 — 6.627 .05437 
10.0 | .5241 .3205 .7022 | — 8.479 — .10579 .03601 





| 
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ture. Fig. 4 also shows the close agreement of 
the Ko’s as calculated for the square well of 
radius e?/mc*? and depth 10.5 Mev and for the 
Gauss error well Ae~*” with A=51.44 mc and 


TABLE II. Values of the coefficient (i.e., 2X/nM) of 
—sin Ky cos Ko for P/Py. This table ts an extension of Table 


a= 21.59 Mmc*/h*, the latter potential well being 
supposed to be superposed on the Coulomb 
potential. A slight change in the depth and 
radius of the square well should bring the two 
curves almost into coincidence. This indicates 


I of BTE. For expansions of X and in terms of E see BTE 


Tables III and V. 


TABLE IV. Values of the coefficient of —sin K, cos Ky 


about +0.1 percent. 











for P/Py. Interpolations linear in E may be correct only to 


















































oy ; =15° ° ° ° ° ° ° 
E (MEv) n @=15° 20 25 30 35 40 45 E (Mev) . @=15° 20° 25° 30° 35° 45° 
5 ‘ J J . . S$ 9 a a a a a aah a aa, ai | “i J 
60 0.068584 2348 4592 8069 13.23 2020 2751 3096 J) Qooi9 1134 1731 2137 2116 1813 5820 0000 
85 0.054228 2 806 5.482 0621 15.77 24.06 32.76 36.85 40 0.079050 13.10 19.99 2467 24.43 17.46 6.122 0.000 
10.0 0.049996 3.050 5.952 10.44 17 "10 26.11 35.54 39.98 50 0.070704 14.66 22.36 27.58 27.30 19.51 6.841 0.000 
, . si =? . , ; = : 70 0.059756 17.35 2645 32.67 32.30 23.07 8.090 0.000 
10.0 0.049996 20.75 31.62 38.99 38.59 27.56 9.664 0.000 

TABLE III. Values of the coefficient (i.e. 4/n?IN+2Y/nN : 
of sin® Ko for P/Pu. This table ts an extension of Table II of ante V. Values of the coefficient of sin? K, for P/Pwy. 

BTE. For expansions of 2¥/n and JM in terms of E see BTE Interpolations linear in. E are very good, 

Tables IV and V. —— 
E(Mev) 9 @=15° 20° 25° 30° 35° 40° 45° 
E (Mev) n @=15° 20° 25° 30° 3s° 40° 45° — - —— 
3.0 0.091279 45.69 1164 2053 2600 218.6 84.38 0.000 
6.0 0.064544 4.225 14.26 36.11 76.17. 137.3. 205.7 238.7 40 0.079050 61.27 155.8 273.7 3468 291.5 112.5 0.000 
7.0 0.059756 4.994 16.73 42.26 89.02 160.4 240.2 278.7 50 0.070704 76.85 1950 3423 433.6 360.1 140.6 0.000 
8.5 0.054228 6.146 20.44 51.51 1083 195.1 291.9 338.7 70 0.059756 108.0 273.4 479.6 607.1 5103 196.8 0.000 
10.0 0.049996 7.303 24.15 60.72 127.6 229.7 343.6 398.7 10.0 0.049996 154.8 390.8 685.8 867.5 729.1 281.1 0.000 
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Fic. 1. The phase shift Ko as a function of the energy E 
of the incident protons for square potential wells of radius 
0.75 e?/mc? and depths 19.690, 19.890 and 20.090 Mev with- 


out interior Coulomb potential. 
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Fic. 2. The phase shift Ko as a function of the energy E 


for square wells of radius e?/mc? and depths 10.305, 10.50 
and 10.60 Mev without interior Coulomb potential. 
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Fic. 3. The phase shift Ko as a function of the energy E 
for square wells of radius 1.25 e?/mc* and depths 6.2452, 
6.3452 and 6.4452 Mev without interior Coulomb potential. 
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E (Mev) 


Fic. 4. The phase shift Ko as a function of the energy E 
for square wells of differing radii (r is in e?/mc?). The depth 
(in Mev) is adjusted for each radius so as to give coinci- 
dence with the experimental Ko at approximately 1.4 Mev. 
There is no interior Coulomb potential. The broken line F 
is for the Gauss error potential well Ae“, where A = 51.44 
mc and a=21.59 Mmc?/h?, cut off at r=3e?/mc? and with 
interior Coulomb potential. 


that the high energy region will not be more 
useful than the low energy region in distinguish- 
ing between the two types of potential wells by 
means of the s wave phase shift. 

Figure 5 shows how the ratio ®& of theoretical 
scattering to Mott scattering varies with the 
energy and the scattering angle. The curves for 
the square well of depth 10.5 Mev and radius 
e?/mc? lie very close to those for the Gauss 
error well. 

The curves in Figs. 6 and 7 show the theo- 
retical number of proton counts as a function 

TABLE VI. Values of the coefficient of —sin Ke cos Kz 


for P/Py. Interpolations linear in E may be correct only 
to about +0.1 percent. 








@=15° 20° 25° 30° 35° 40° 45° 
6.177 3.390 — 5.764 —22.78 —43.35 —53.57 


E 
(MEv) n 
3.0 0.091279 5.241 





4.0 0.079050 6.060 7.150 3.929 — 6.686 —26.44 —50.33 —62.21 
5.0 0.070704 6.780 8.005 4.402 — 7.495 —29.65 —56.46 —69.78 
7.0 0.059756 8.029 9.486 5.221 — 8.896 —35.20 —67.05 —82.88 
10.0 0.049996 9.602 11.35 6.251 —10.66 —42.17 —80.37 —99.34 








TABLE VII. Values of the coefficient of sin? Kz for P/Pwx. 
Interpolations linear in E are very good. ( For the calculation 
of d scattering there is a third term, which is +(40/n*) sin Ko 
sin Ke cos[K2—Ko+202—200]. See BCP Eg. (6.7).) 








@=15° 20° 2s° Ss 35° 40° 45° 


14.03 525.9 738.7 
19.06 704.4 988.7 
24.08 299.4 882.8. 1238.7 
34.12 421.0 1239.6 1739 
49.19 603.4 1775 2489 


E (Mev) n 


3.0 0.091279 22.72 
0.079050 30.23 
0.070704 37.74 
0.059756 52.77 
0.049996 75.31 





7.085 

9.290 
11.498 
15.91 
22.53 


27.37 
36.31 
45.26 
63.14 
89.97 


Sam 
oocoo 

















-E 
th 
ci- 
Vv. 


th 


-2 


Oe! 


— —— 





PHASE SHIFT CALCULATIONS 1197 


of the energy of the incident protons. The 
number given is the number expected from the 
scattering chamber used by Herb, Kerst, Parkin- 
son and Plain* per microcoulomb of incident 
proton current per mm of oil hydrogen pressure 
in the chamber. Fig. 6 is for the square well of 
depth 19.690 Mev and radius 0.75 e?/mc? without 
interior Coulomb potential; added as broken 
lines are the curves for the Gauss error potential 
well Ae~*", where A=51.44 mc? and a=21.59 
Mmc*/h*, with interior Coulomb potential. Fig. 7 
is for the square well of depth 6.3452 Mev and 
radius 1.25 e?/mc? without interior Coulomb 
potential; the dotted lines are the same Gauss 
error curves as in Fig. 6. For a similar graph for 
square wells of radius 0.5 and 1.0 e?/mc*? see BTE 
Fig. 14. In each case pure s state scattering is 
assumed. 








$s 6 7 & 10 


4 
. Ejney 
Fic. 5. The ratio ® to Mott scattering, as a function of 
the energy E and the scattering angle 9, of the scattering 
calculated for: A, Gauss error well Ae-*”, where A =51.44 
mc, and a=21.59 Mmc*/h*?; B, square well of radius 0.75 
e?/mc* and depth 19.690 Mev, and C, square well of radius 
1.25 e?/mc* and depth 6.3452 Mev. The curves for a square 
well of depth 10.50 Mev and radius e?/mc* would lie very 
close to A. 





5 6 7 8 9 10 


2 4 
E (rev) 
Fic. 6. The number of proton counts expected plotted as 
a function of the energy E of the incident protons and as a 
function of the scattering angle @. The unbroken curves 
are for the square well of radius 0.75 e*/mc* and depth 
19.690 Mev without interior Coulomb potential. The 
broken-line curves are for the Gauss error well Ae-@”, 
where A =51.44 mc?, and a=21.59 Mmc*/h? with interior 
Coulomb potential. For curves for the square well of 
radius e?/mc? and depth 10.50 Mev, see BTE Fig. 14. The 
number plotted is the number of counts expected from the 
scattering chamber of HKPP* per microcoulomb of incident 
proton current per mm of oil hydrogen pressure in the 
scattering chamber (see BTE). Pure s state scattering is 
assumed. 


Figures 6 and 7 show that the sensitivity of 
the number of counts to a change in the range of 
nuclear force is no greater at higher energies than 
at energies below 3 Mev. (It must be remembered 
that the energy of the point of intersection of 
two phase shift curves for different ranges of 
force is determined by the depth of the wells.) 
For the larger values of the scattering angle 0 
experiments at energies from 1 to 3 Mev may 
be best for determining the range of the well. 
(See BTE Fig. 14 and discussion on p. 1059 
regarding small angle scattering.) Thus it appears 
that as far as the s state scattering is concerned 


TABLE VIII. Values of K, for the constant potential +10.5 Mev of radius e*/mc*. 








ENERGY OF INCIDENT PROTONS 





(MEv) 3 4 5 6 7 & 9 10 
Ki (attractive potential) 0.40° 0.63° 0.90° 1.1(9)° 1.5(4)° 1.8(7)° 2.2(3)° 2.6(1)° 
—0.67° —0.85° — 1.0(4)° —1.2(3)° —1.4(2)° 


Ki: (repulsive potential) —0.23° —0.36° —0.51° 
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Fic. 7. The same type of curve as Fig. 6. The unbroken 
curves are for a square well of radius 1.25 e?/mc? and depth 
6.3452 Mev. The broken-line curves are the same as in 
Fig. 6. 


the chief value of high energy experiments will 
be -to corroborate the results of the low energy 
experiments and to make possible better determi- 
nations of the range and depth of the potential 
well by increasing the length of the Ko, E curve 
to be fitted by theory. The determination of Ko 
may be more complicated because of the presence 
of phase shifts K, and Kz. 

The experiments at energies above 3 Mev are, 
however, expected to be very valuable in deter- 
mining the p state interaction (see BTE pp. 
1059-61). In Table VIII are given values of the 
phase shift K, calculated for a constant potential 
of +10.5 Mev with a range of e?/mc? and without 
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interior Coulomb potential. The second and 
third rows of the table correspond, respectively, 
to attractive and repulsive p interactions. In 
both cases it is seen that the absolute values of 
K;, increase with energy as would be expected. 
For the Gauss error potential +Ae-*” with 
A=51.44 mc and a= 21.59 Mmc*/h? superposed 
on the Coulomb potential, one finds, by inte- 
grating the wave equation numerically and join- 
ing to the Coulomb function at r=3 e?/mc?, that 
for a repulsive p state interaction K,= —0.74° at 
5 Mev, while for an attractive interaction 
K,=+1.23° at the same energy. These values 
are larger than the corresponding values for the 
constant potential because the Gauss error 
potential extends to larger values of r. Beyond 
showing that the values of K, are probably 
large enough to appear in high energy scattering 
experiments, these values of K, have little direct 
significance, since it is improbable that the 
singlet s state and triplet p state interactions are 
the same. 

If the Coulomb potential is considered to be 
effective inside a square potential well, the depth 
of the well must be increased in compensation 
if one wishes to obtain the same phase shift Ko 
as was obtained without the interior Coulomb 
potential. The necessary increase in depth of the 
square well of depth 10.5 Mev and radius e?/mc* 
was found to be 0.849 Mev at the energy 7 Mev, 
and 0.860 Mev at the energy 10 Mev. The calcu- 
lations were made by the use of BTE Eq. (11.2) 
and checked by BTE Eg. (11.5). It may be 
noted that, as a very rough empirical rule, 
6D =0.826+0.00333 (E—0.2), for this well, 
where 6D and E£ are in Mev. 

The authors are glad to express their thanks 
to Professor G. Breit for very helpful discussion 
and advice. 
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Measurements of the first and second Townsend coefficients for ionization by collision 
have been carried out for surfaces of Ni and Al in an atmosphere of hydrogen. The Al 
surface was evaporated in an atmosphere of argon and is estimated to be a few hundred 
molecular layers thick. The values found for the a- and y-coefficients are compared with the 
results previously reported for a Pt surface in hydrogen. The values of a/p, in pure hydrogen, 
plotted as a function of X/p show a decreasing characteristic between values of X /p from 900 
to 1400. The curves of values of y plotted as a function of X/p show a characteristic photo- 
emission peak at an X/p of 131, and a gradual rise at the higher values of X /p. Sparking poten- 
tial curves are calculated from the y-values and are compared with the experimentally deter- 


mined curves. 


INTRODUCTION 


EASUREMENTS of the Townsend ioniza- 

tion coefficients for Pt and NaH surfaces 
in an atmosphere of hydrogen have been pre- 
viously reported. These measurements have 
been extended to surfaces of Ni and Al. An 
improved ionization chamber has been con- 
structed with which it was possible to extend 
the measurements of a/p to an X/p of 1700 
volts/em per mm Hg. 

A diagram of the ionization chamber is shown 
in Fig. 1. The chamber is of glass 18 cm in 
diameter and 45 cm in length. Fitted inside the 
glass wall of the chamber is a grounded wire 
gauze which is not shown in the diagram. The 
electrodes are of Ni. The rod supporting the 
upper electrode is threaded and this electrode 
can thus be raised or lowered by turning the 
armature with an exterior electromagnet. The 
electrodes are 6 cm in diameter and the distance 
between them may be varied from zero to 6 cm. 
A quartz window is attached to a side arm by a 
graded seal. Light from a quartz mercury lamp 
passed through a quartz lens and the window and 
then fell upon a spot near the center of the 
cathode. This radiation gave the photoelectric 
current i at the surface of the cathode. From 2 
to 10 cm? of the cathode surface was illuminated, 
the larger area being used in the case of the Ni 
cathode. The data taken with the larger area of 
illumination showed some irregularities when the 


1 Donald H. Hale, Phys. Rev. 55, 815 (1939). 


plate separation was greater than 2 cm. This was 
taken to indicate that in these cases there was 
some field distortion at the greater values of the 
plate separation. 

The Al surface was evaporated directly upon 
the Ni surface of the cathode. This evaporation 
was carried out in an atmosphere of argon at a 
pressure of 0.8 mm of mercury. The side arm, 
shown sealed off in Fig. 1, was extended to 32 cm 
in length and carried a support for a filament and 
a soft iron armature. By means of an external 
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Fic. 1. The ioni- 
zation chamber. @ is a, 
asoft ironarmature, 


c the cathode and gq 
a quartz window. 
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Fic. 2. The values of a/p plotted as a function of X/p. 


electromagnet the filament could be moved to a 
position about 1.5 cm above the cathode and 
after the evaporation was complete the filament 
and its supports retracted into the side arm. The 
filament consisted of three 20-mil tungsten wires 
loosely twisted together. To this filament were 
fastened several small pieces of aluminum plate 
0.4 mm thick and about 0.25 cm? in area. The 
metal was analytical grade and contained 0.32 
percent of iron and 0.45 percent of SiOQe as 
impurities. The filament was first heated while 
retracted into the side arm until the sheets 
melted. This treatment removed the surface con- 
tamination which was deposited on the walls of 
the side arm. 

Tank hydrogen was used. The gas was purified 
by passing it through a purifying train in the 
conventional manner and then through two 
liquid-air traps and a fine leak before it entered 
the chamber. The photoelectric current was kept 
at a minimum at all times in order that the 
measurements might not be falsified by space 
charge effects. 

The ionization chamber was built into an oven 
and was baked out at a temperature of 425°C for 
16 hours each time air was admitted. A discharge 
in hydrogen was run for one hour after each 
baking and before any data were taken. The 
reason for this treatment of the cathodes was 
that in the presence of hydrogen neither the clean 
Ni surface nor the freshly evaporated Al surface 
was stable photoelectrically but the photoemis- 
sion changed continually with time and in these 
measurements it was essential to have a stabilized 
surface. Thus these surfaces were carefully 
formed probably as hydrogen coated surfaces of 
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some stability. In each case this treatment in- 
creased the photoelectric sensitivity of the 
surface. In the case of Al the treatment also 
apparently shifted the photoelectric threshold 
towards the long wave-length side of the spec- 
trum as it was found that this surface was sen- 
sitive to visible light after the discharge had been 
run. It was therefore necessary in all measure- 
ments on the Al surface to have the ionization 
chamber covered by a light-tight box. The Al 
surface was not baked after it was deposited 
upon the surface of the Ni cathode. 

In this type of measurement it is necessary 
that a fairly steady source of ultraviolet radiation 
be used as a source of the photoelectric current Zo. 
It was found that a lamp designed by Mr. E. H. 
Guyon of this laboratory was quite satisfactory 
when used with a stabilizer of the iron wire- 
hydrogen type.” 

All data given in this paper are reduced to a 
temperature of 22°C. 


EXPERIMENTAL RESULTS 


The Al and Ni surfaces gave the characteristic 
curves, as shown by Bowls,* when (1/ p) In (7/70) 
was plotted as a function of the plate separation 
x. The exact amount of Al evaporated in forming 
the Al cathode could not be determined but it 
was estimated that the surface was a few 
hundred molecular layers thick. The values of 
a/p plotted as a function of X/p are shown in 
Fig. 2. Here a is the first Townsend ionization 
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Fic. 3. The values of y plotted as a function of X/p for the 
Ni cathode. 


2 Francis A. Jenkins and Harvey E. White, Fundamentals 
of Physical Optics (McGraw-Hill Book Company, 1937), 


p. 248. 
3 W. E. Bowls, Phys. Rev. 53, 293 (1938). 
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coefficient and represents the number of electrons 
produced by each electron in 1 cm of path in the 
direction of the electric field. X is the field 
strength in volts/cm and p is the pressure in mm 
Hg. The Al cathode gave values of a/p which lay 
along this curve within experimental error. The 
values of a/p found for the Pt surface! also agree 
with those shown in Fig. 2. In the case of the Pt 
surface the measurements were reliable only to 
an X/p of about 900. With the new ionization 
chamber reliable measurements of a/p could be 
extended to an X/p of 1700. This curve shows a 
falling characteristic between values of X/p of 
900 to 1400. Above an X/p of 1400 the curve 
starts a second sharp rise. 

The values of the second coefficient y were 
computed from the equation of Townsend‘ which 
was also derived by Thomson,° 


{=ig—_—_—_—_. (1) 


This equation assumes a secondary emission 
process at the cathode. Loeb® has shown that 
under the conditions obtaining in the ionization 
chamber in this work it is likely that Eq. (1) 
properly describes the mechanism. This may be 
either a secondary emission of electrons at the 
cathode due to positive ion bombardment or a 
photoelectric mechanism. The values of the 
second coefficient plotted as a function of X/p 
are shown in Figs. 3 and 4. These curves show 
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Fic. 4. The values of y plotted as a function of X/p for the 
Al cathode. 


4J.S. Townsend, Electricity in Gases (Oxford University 


Press, 1914), p. 314. 

5J. J. Thomson and G. P. Thomson, Conduction of 
Electricity through Gases (Cambridge University Press, 
third edition, 1933), Vol. II, p. 518. 

6 L. B. Loeb, Rev. Mod. Phys. 8, 267 (1936). 
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Fic. 5. The sparking potential curve for the Ni cathode. 
The crosses are experimentally determined points. The 
circles are values of the sparking potential calculated from 
the y values determined in this work. 


the same general characteristics. It is quite prob- 
able that the values of y for the Ni surface are 
somewhat less accurate than those for Al. In the 
case of Ni the larger area of the cathode was 
illuminated and there were indications that at 
plate separations above 2 cm there was a small 
field distortion. Such a distortion would tend to 
lower the measured values of y but these values 
are of the right order of magnitude. This error 
would not alter the form of the curve. There is a 
pronounced peak at the lower values of X/p and 
at the higher values of X/p both curves show a 
general rise. The maxima, at the lower values of 
X/p, occur, within experimental error, at an X/p 
of 131. The maximum in the similar peak as 
shown by the Pt surface in hydrogen was at an 
X/p of 120 which is nearly the same. The maxi- 
mum value of y, at an X/p of 131, for the Al 
surface is greater than for the Ni surface. How- 
ever, this may be due to the field distortion which 
is suspected in the case of the Ni surface. The 
ratio of the value of y for Al to that for Ni, at an 
X/p of 131, is 1.25. 

In the range of values of X/p from 300 to 
1400 the value of y for the Al surface is smaller 
and the rise with increasing values of X /p is more 
gradual than is the case for the Ni cathode. 
However, above an X/p of 1400 the Al surface 
gives y values which increase more rapidly than 
is the case for the Ni cathode. This again may be 
due to field distortion as such distortion in 
falsifying a would have more effect on the value 
of y at the lower values of X/p. 
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The sparking potentials for both surfaces in 
hydrogen were experimentally determined. These 
data are shown in Fig. 5 for the Ni cathode. The 
crosses are the experimentally determined spark- 
ing potentials. The dotted line at the lower values 
of pé is taken from data published by Quinn.’ 
The circles in this figure are calculated break 
down potentials which are found by use of the 
values of y shown in Fig. 4. The agreement is 
good except at the minimum values of the spark- 
ing potential. The data for the Al surface gave 
about the same sort of agreement. 


DISCUSSION 


The a/p curve shows the same values for the 
Pt, Ni and Al cathodes as it should since this 
coefficient is a function of the gas between the 
electrodes. Above an X/p of 900 the values of 
a/p show a decline and then rise again starting 
at an X/p of 1400. The decline is analogous to 
the results of Huxford*® for argon-filled tubes. 
This behavior may logically be ascribed to a 
decrease in the cross section for ionization by 
electron impact at the higher energies as shown 
by the Compton-van Voorhis curves. The effect 
observed at the higher values of X/p is not, in 
any case, to be confused with the decline ob- 
served in hydrogen contaminated with a volatile 
sodium compound! which is caused by a second- 
ary ionization of the volatile compound by 
hydrogen photons. The subsequent rise, as 
shown by the curve of Fig. 2, above an X/p of 
1400 is unexplained. If as conjectured the decline 
is a reflection of the decline of probability of 
ionization with increasing electron energies, one 
must assume that the subsequent rise in a/p is 
caused by the failure of the electrons to gain 
their terminal energy in the gap at the pressures 
used and thus the increase in a/p is due to the 
fact that below the maximum terminal energy 
their chance to ionize is greater. The pressures 
of 0.1 mm used in determining the values of a/p 
in the range of values of X/p from 900 to 1400 
make this explanation plausible. However, at the 
rise in the curve the pressures were reduced by 
less than 10 percent below those near the 
minimum. Hence it is possible that the entire 
decline may be caused by failure of the electrons 


7 Robert B. Quinn, Phys. Rev. 55, 482 (1939). 
8 W.S. Huxford, Phys. Rev. 55, 754 (1939). 





H. HALE 


to achieve their terminal energy and that the 
hump corresponding to the Compton-van Voorhis 
curve may occur at much higher values of X/p. 

It is seen that, except for minor differences due 
in part to experimental errors, the values of y 
for the Ni and Al surfaces are nearly the same. 
There is for all three surfaces a definite photo- 
electric peak. These peaks occur at the same 
value of X/p for the Ni and Al surfaces and at 
nearly this value of X/p for Pt. However, it is 
to be remembered that a discharge was run in 
hydrogen, in the case of the Ni and Al surfaces 
before any data were taken. It must also be 
remembered that Pt is notoriously sensitive to 
absorption of hydrogen as shown by Langmuir 
in 1918.° The values found for y for Pt surfaces 
in nitrogen*® and for NaH surfaces in hydrogen! 
are not alike nor do they resemble the values 
found for the Pt, Ni and Al surfaces in hydrogen. 
This indicates the probability that for secondary 
emission by positive hydrogen ion bombardment 
surfaces of Pt, Ni and Al in hydrogen are essen- 
tially the same. 

In a measure the similarity of these curves for 
Pt, Ni and Al cathodes in hydrogen for y as a 
function of X/p is disappointing as it was hoped 
to find that different surfaces had markedly 
different properties. Since, however, Ni and Pt 
are noted catalysts for hydrogenation reactions, 
this result might have been expected. The be- 
havior of Al is interesting in that it indicates its 
affinity for hydrogen. Further less reactive metals 
should be deposited by the Beeck"” process and 
the y values determined for such cathodes. It is 
unfortunate that hydrogen, one of the few gases 
free from metastables, should be so chemically 
active with metals. This fact must lead one to 
seek another nonreactive vehicular gas free from 
metastables in which to test the intrinsic proper- 
ties of surfaces for secondary electron libration 
in the Townsend gap. 

The calculated values of spark breakdown at 
low values of pé agree satisfactorily with the 
experimentally determined values. 

It is a pleasure to acknowledge that Professor 
L. B. Loeb has contributed much to this paper 
by many helpful discussions of the interpretation 
of the experimental results. 

®1I. Langmuir, J. Am. Chem. Soc. 40, 1361 (1918). 


10Otto Beeck, A. Wheeler and A. E. Smith, Phys. Rev. 
55, 601(A) (1939). 
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The scattering of electrons in lead and in aluminum has been studied, mainly under the con- 
ditions used previously for the study of energy loss. The primary object of the investigation 
was to determine the extent to which the existing energy loss data (up to 13.5 Mev) are affected 
by multiple scattering. The conclusion reached is that the measurements so far made in lead 
with electrons below 9 Mev are so much affected by multiple scattering within the absorbing 
material that they are of very little use as a direct check on the theory for energy loss. In the 
measurements from 9 to 13.5 Mev the scattering is found to be small enough so that a com- 
parison of the energy losses with the theoretical values is possible. The observed losses are 
roughly 40 percent in excess of the theoretical. Some measurements on the multiple scattering 
of 0.9-Mev electrons in thin sheets of aluminum were made, and it was found that the most 
probable angle of scattering and the average angle of scattering were in good agreement with 


theoretical predictions. 





INTRODUCTION 


CONSIDERABLE amount of data is now 

available on the loss of energy suffered by 
electrons passing through solid materials such as 
lead, carbon and aluminum. During the past two 
years the work carried on in this laboratory has 
covered the energy range from 0.5 Mev to 17 
Mev.'-* Results from the University of Cali- 
fornia’ and the California Institute of Tech- 
nology® are in essential agreement with those 
obtained here. It may be said with some cer- 
tainty that, for the experimental method used 
at present, the results on energy loss are well 
established. The main uncertainties lie in the 
interpretation of the results, and one difficulty 
which overshadows all others is the estimation of 
the actual distance traveled by an electron inside 
the absorber. 

On the assumption that the measured thick- 
ness of a given absorber corresponds to the 
actual path traversed by an electron, the experi- 
mental results have consistently shown losses 
greater than those predicted by theory. The 
discrepancy is not the same over the whole range 
but varies from a factor of more than two at 


J. J. Turin and H. R. Crane, Phys. Rev. 52, 63 (1937). 
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J. J. Turin and H. R. Crane, Phys. Rev. 52, 610 (1937). 
~ J. Ruhlig and H. R. Crane, Phys. Rev. 53, 618 (1938). 


. R. Curtis, Phys. Rev. 53, 986 (1938). 

5L. J. Laslett and D. G. Hurst, Phys. Rev. 52, 1035 
(1937). 
(1938) A. Fowler and J. Oppenheimer, Phys. Rev. 54, 320 


the low energies to about one and a half at the 
high end of the range. The question is: How 
much of this discrepancy is due to the fact that 
the path taken by an electron does not corre- 
spond to the actual thickness of the absorber? 
The number we wish to find is called the path 
length-thickness ratio, or the ratio of the actual 
path length to the thickness of the absorber. 
This ratio is greater than unity for all absorbers, 
because the electron suffers an extremely large 
number of deflections, taking a zigzag path. In 
cases where there is much scattering the direction 
of emergence of a particular electron from the 
material tells little or nothing about its path 
length, since its final direction is a result of very 
many deflections in random directions. It there- 
fore seems to be of no advantage to restrict the 
data to those electrons which emerge in a direc- 
tion normal to the surface, or to try to make 
individual allowance for the effective thickness 
traversed by each electron, depending upon its 
angle of emergence. The best procedure seems to 
be to include all tracks in the data and to try to 
make a correction which applies to the average 
path length. 

The complexity of the theoretical treatment of 
the path length-thickness ratio may be indicated 
by listing some of the important facts to be 
taken into account. 

(1) The probability of a deflection varies with 
the size of the deflection. 

(2) The probability of a deflection of a given 
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Fic. 1. Graphical device which is used in obtaining the 
true (three-dimensional) angular distribution of scattered 
electrons, from the cloud-chamber data. 


size is a function of the energy of the electron, 
while the energy of the electron changes at an 
appreciable rate along its path due to two 
processes : ionization and radiation. 

(3) The electrons finally must be separated 
into two classes: those which lose their entire 
energy in the absorber, and those which succeed 
in passing through the absorber. It is the average 
path length for the latter group alone that must 
be found, if the result is to be compared with the 
experimental data that are available. 

(4) The situation has been made even more 
difficult to handle because the validity of the 
single scattering formulae has been placed in 
some doubt recently,’~" both in regard to large 
angles of scattering and in regard to the small 
angles which are important in multiple scattering. 
This is of direct concern because the theoretical 
estimates of path length are derived necessarily 
from the laws of single scattering. 

The calculations can be simplified if the scat- 
tering is so small that at all times the deviation 
from the initial direction is not great (sin @=6@). 
The treatment simplifies also if the scattering is 
so great that the mathematics of diffusion become 
applicable. It is desired, in the measurement of 
energy loss, to approach as nearly as possible the 
condition in which the scattering is small. Un- 
fortunately the choice of an absorber which is 
thin enough to insure the electron path is nearly 
straight, and yet thick enough to cause a loss of 
energy large enough to be accurately measurable 
seems to be difficult in the range of energy with 
which we have been concerned, at least when a 


7 A. Barber and F. C. Champion, Proc. Roy. Soc. A168, 
159 (1938). 

SE. J. Williams, Proc. Roy. Soc. A169, 531 (1939). 

*N.L. Oleson, K. T. Chao, J. Halpern and H. R. Crane, 
Phys. Rev. 56, 482 (1939); Phys. Rev. 56, 1171 (1939). 
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material of high atomic number is used. However, 
with the aid of the experimental results which 
we shall present, we have been able to arrive at 
some general statements which simplify the 
problem somewhat, at least in regard to the 
interpretation of the energy loss measurements 
already existing. 


PROCEDURE 


Electrons were allowed to enter the cloud 
chamber through a thin window in the side wall 
and to fall normally upon a piece of absorbing 
material placed across the center of the chamber, 
as described in previous papers.'~*. The energy 
and the angle of deflection in the horizontal pro- 
jection were measured for each track passing 
through the absorber. Although stereoscopic 
pairs of photographs were taken in all cases, it 
was found more satisfactory to measure only the 
component of deflection in the horizontal plane 
and to convert the distribution thus obtained 
into a three-dimensional one by a graphical 
method. This conversion can be made uniquely, 
because of the cylindrical symmetry about the 
initial direction of motion of the electrons. The 
reason for adopting this procedure rather than 
the direct three-dimensional measurement from 
the stereoscopic pictures is the following: 
Because of the limited depth of the chamber 
many tracks deflected through large angles away 
from the plane of the chamber are invisible or at 
least unmeasurable. Therefore an important 
correction would have to be made in this case, 
and the amount of the correction would be 
difficult to determine. By restricting the data to 
those tracks which lie within definite, small 
angular limits in the vertical direction (the plane 
of the chamber is horizontal), the correction to 
be applied becomes more definite and calculable. 
The data pertaining to the geometry are as 
follows. 


TABLE I. Area ratios for angles larger than those given in 











Fig. 1. 
ANGLE PERCENT ANGLE PERCENT ANGLE PERCENT 
26-28 42 36-38 28 46-48 22 
28-30 37 38-40 26 48-50 21 
30-32 35 40-42 24 50-52 21 
32-34 32 42-44 23 52-54 20 
34-36 30 44-46 23 54-56 20 
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(1) The depth of the visible region of the 
chamber is well defined by the light beam and is 
taken to be 3 cm. 

(2) The incident electrons are passed through 
a slit system outside the chamber, so that their 
paths lie within 5 degrees of the horizontal plane 
and so that they strike the center of the absorber. 

(3) The scattering is cylindrically symmetric 
about the direction of the incident electron. 

(4) Only those tracks were included in the 
data which, after passing through the absorber, 
did not pass out of the light beam. This means 
that only those tracks were included whose 
vertical component of scattering was smaller 
than 10 degrees. All angles in the horizontal plane 
(up to 90 degrees) were included. 

For the actual operation of making the con- 
version to the three-dimensional case the 
diagram in Fig. 1 is used. This shows the fraction 
of each scattering cone which falls within the 
horizontal and vertical angles defined above, and 
together with the supplementary data in Table I 
gives the numerical values necessary for con- 
version of the data to the three-dimensional case. 
The actual measurements yield the number of 
tracks in a particular vertical strip. This number 
then has to be distributed among those parts of 
the circular zones which lie within the vertical 
strip. It is necessary to start with the largest 
angle and work toward the center. An example of 
the procedure is as follows: Suppose we measure 
the number of tracks which fall in the rectangular 
column between 26 and 28 degrees. Since the 
angle is large, we can assume that the density of 
population in this rectangle is uniform, and 
therefore we can determine the number of tracks 
to assign to the shaded area N. Let us call this 


TABLE II. Ratio of the complete ring about the origin to that 
of the curved figure. 











ANGLE FACTOR ANGLE FACTOR ANGLE FACTOR 
10-12 $.2 30-32 18.5 50-52 28.0 
12-14 7.2 32-34 19.6 52-54 28.7 
14-16 9.3 34-36 20.6 54-56 29.4 
16-18 10.5 36-38 21.6 58-60 30.8 
18-20 11.7 38-40 22.5 62-64 32.1 
20-22 12.9 40-42 23.5 66-68 33.5 
22-24 14.1 42-44 24.5 70-72 34.4 
24-26 15.2 44-46 25.4 74-76 35.1 
26-28 16.3 46-48 26.3 78-80 35.6 
28-30 17.4 48-50 27.2 82-84 35.8 
86-88 36.0 

















number N. The density of population in the 
unshaded part of the 24-26-degree rectangle is 
the same as that of the area N, because of the 
cylindrical symmetry about the origin. The 
number of tracks in the unshaded part of the 
24-26-degree rectangle is therefore 0.47N, be- 
cause this is the ratio of the areas. The number 
of tracks which belong to the shaded area M is 
now determined by counting the number in the 
24-26-degree rectangle and subtracting 0.47N. 


TABLE III. Energy distribution of electrons striking 











aluminum, 
THICKNESS OF 

ALUMINUM NUMBER ENERGY MEV 
31 0.6 to 0.8 

0.0025 cm 87 0.8 to 1.0 
32 1.0 to 1.2 
45 0.6 to 0.8 

0.01 cm 78 0.8 to 1.0 
35 1.0 to 1.2 
43 0.6 to 0.8 

0.025 cm 99 0.8 to 1.0 
53 1.0 to 1.2 








The number in the shaded area L is equal to the 
number in the 22—24-degree rectangle, minus 
0.52M. By continuing this process to zero 
degrees, the density of population in all parts of 
the diagram is found. A continuation of these 
area-ratios for larger angles than those which 
appear on the diagram is given in Table I. 

To get the data finally in terms of the number 
of tracks in a given interval of solid angle or 
interval of angle, we make use of Table II. This 
table gives the ratio of the complete ring around 
the origin to that of the curved figure (for 
example the shaded one whose base lies between 
10 and 12 degrees in the diagram). 


EXPERIMENTAL RESULTS 


The experimental results are represented by 
the curves in Figs. 2 and 3. In these figures the 
points represent number of tracks scattered into 
an angular interval between ¢ and ¢+ Ad, in the 
plane of the chamber. In each case two curves 
are given. The full line curve is drawn smoothly 
through the experimental points, the dashed 
curve is the angular distribution which includes 
the entire cone of scattering, obtained by the 
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Fic. 2. Angular distribution of electrons scattered by aluminum. 


use of Fig. 1. In Fig. 2 an additional curve is 
given to represent a theoretical distribution of 
the form: y= Axe~8** with the parameters chosen 
arbitrarily. The reason for including the latter 
curve will be indicated in the discussion which 
follows later in the paper. The electrons which 
were incident upon the aluminum scatterers had 
an average energy of 0.9 Mev, and the actual 
distribution in energy is given in Table III. The 
electrons incident upon the 0.0038-cm and 
0.0066-cm lead scatterers had an average energy 
of 0.9 Mev and a distribution as given in 


Table IV. 


THEORETICAL CONSIDERATIONS 


There can be little doubt that in the experi- 
mental data we have presented, we are dealing 
mainly with the effects of multiple scattering, 
even for the fastest electrons and the thinnest 
absorbers. In such cases, where a large number of 
deflections occurs inside the absorber, the dis- 
tribution in direction of the emergent electrons 
must be found by compounding statistically 
effects of the large number of deflections. The 
result can be represented approximately by the 
simple error curve of the form 


N(0)d0= A 6e-2"d0. 


However, there is always a small number of 
large deflections due to single encounters, whose 
effect is to distort the distribution from the form 
mentioned above, especially at large angles. One 
must expect a noticeable deviation from the 
error curve at angles which are several times the 
most probable angle; this is the region of the 
so-called plural scattering. At extreme angles the 


distribution is expected to correspond closely to 
that predicted by the formula for single scatter- 
ing. The angle beyond which the single scattering 
formula applies can be estimated by the use of 
Wentzel’s criterion." The 4w used by Wentzel 
is approximately equal to the most probable 
angle, and it is customary to use 3 or 4 times 
this as the limiting angle. Referring to our 
curves, we may say that the scattering is essen- 
tially single at angles equal to about 4 times the 
angle at which the peak occurs in the dashed 
curve. The experimental distributions will be 
expected to conform to an error curve in the 
region of the maximum, and to have a plural or 
single scattering ‘‘tail’’ at the large angles, which 
will lie above the tail of the error curve. 

It is essential at this point to tell what happens 
as the most probable angle approaches 25 or 30 
degrees. In cases of thick absorbers, where the 
scattering is very great, the electrons will lose 
completely their original directions of motion, 
and take up random directions inside the ab- 
sorber. Thus the case will resemble diffusion. 
Obviously any further increase in the thickness 
of the absorber will not alter the distribution in 
direction of the emerging electrons, but will only 


TaBLE IV. Energy distribution of electrons striking lead. 
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11G, Wentzel, Ann. d. Physik 69, 335 (1922). 
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Fic. 3. Angular distribution of electrons scattered by lead. 


change the number emerging. We have found 
experimentally that in those cases in which dif- 
fusion seems to exist, the most probable angle 
of emergence of the electrons is 25 to 30 degrees 
(measured from the normal to the surface). This 
means that if we start with a very thin absorber 
and gradually increase its thickness, we will at 
first get an error curve for the distribution, whose 
peak will move toward the right with increasing 
thickness. As the peak approaches 25-30 degrees 
the form of the distribution will change over into 
one which is characteristic of diffusion. No 
matter how much more we increase the thickness 
after this point has been reached, the curve will 





not change in any way except in total intensity. 
The arguments of the preceding paragraphs, 
which are based upon the error curve, will 
therefore apply only if the most probable angle 
is considerably less than 25 degrees. 

Bethe, Rose and Smith" have treated the 
problem of multiple scattering and have given a 
theoretical expression for the “transport mean 
free path”’: 

1 2xrNZ*etW? 2ap 


= — log —— 


d (W2—mic!)? Zt 


12H. A. Bethe, M. E. Rose and L. P. Smith, Am. Phil. 
Soc. Proc. 78, 573 (1938). 
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TABLE V. Comparison of our observed values of O9max with those calculated by Bethe, Rose and Smith. 
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AVERAGE 


THICKNESS, ENERGY, Onax 
MATERIAL cM MEV CALC. 
Al 0.0025 0.9 3) 
Al 0.01 0.9 16 
Al 0.025 0.9 25 
Pb 0.0038 0.9 
Pb 0.0066 0.9 
Pb 0.05 3.0 
Pb 0.05 5.0 
Pb 0.05 7.0 


where W is the energy of the electron including 
the rest mass, a is the Bohr radius and p is the 
momentum of the electron. The quantity X is 
immediately connected with the most probable 
angle of scattering : 


Omax = (2¢, /r)). 


The above formula can be applied to our data 
only when the most probable angle is consider- 
ably less than 25 degrees. Under this condition 
Bethe, Rose and Smith give as the mean path 
length for the electron in the absorber, 


t=to(1 + 4G nex), 


where ¢ is the actual path length taken by the 
electron, ¢) is the measured thickness of the 
absorber and @max is the most probable angle of 
scattering. 

Williams" has recently treated both single and 
multiple scattering in detail, and has developed 
formulae for the angular distribution and also 
for the average angle of scattering. He has com- 
pared our results for the average angle of scat- 
tering in 0.0025 cm aluminum at 0.9 Mev with 
that calculated by means of his formula and finds 
fairly good agreement. Our experimental value 
for the average angle in the plane projection is 
5.5 degrees and his calculated value, after apply- 
ing the appropriate correction for our geometrical 
conditions, is 5.9 degrees. This indicates that 
theory cannot be seriously wrong for the scatter- 
ing in aluminum at this energy. 


DISCUSSION OF RESULTS 


Scattering 
A comparison between our results and the 
values of @nax calculated by the formulae of 


SE. J. Williams, Proc. Roy. Soc. A169, 531 (1939). 





PATH LENGTH ENERGY 
9unax THICKNESS LOSS, 
OBS. RATIO, CALC. EXP./THEOR. 

6 1.01 —_ 
11 1.04 — 
19 1.10 1.25 
22 3. 
25 2.7 

| 
29 Diffusion 18 
21 1.8 
21 1.8 


Bethe, Rose and Smith is shown in Table V. 
With the help of this the following can be said: 

(1) The angular distributions of the electrons 
are those which would be expected to result 
mainly from multiple scattering. The positions 
of the maxima and the shapes of the curves are 
in fair agreement with the predictions of the 
multiple scattering formulae, in those cases in 
which the theory is applicable, namely @mnax<25°. 

(2) The experiments on 0.0025 cm and 0.01 
cm aluminum lie within the region of applica- 
bility of the Bethe, Rose and Smith formula, and 
we see that here the experimental and calculated 
values of @nax are in fair agreement. The energy 
loss was too small to be measured in these two 
absorbers. 

(3) The energy loss discrepancy in the 0.025 
cm aluminum is somewhat larger than is ac- 
counted for by the increase in path length which 
is calculated, but this probably indicates that at 
this value of @max (19 degrees) we are already 
outside the region in which the formula gives the 
correct path length. 

(4) For all the lead absorbers listed in Table V 
the scattering is so great that it approaches the 
case of diffusion. It is important to keep in mind 
the fact that as the amount of scattering in- 
creases Omax approaches the stationary value of 
about 25 degrees, while the path length-thickness 
ratio continues to rise. Therefore in all cases in 
which @max is near 25 degrees, the formula gives 
only a lower limit to the path length. The value 
25 degrees for the limit which the angle ap- 
proaches was obtained by referring to the expert- 
mental data, and may therefore be to some extent 
characteristic of the particular set-up used for 
these measurements. 
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(5) Table VI gives a resumé of results from 
other papers. Some similar data exist for posi- 
trons,*: +: © but these have not been included in 
the table. In two of the measurements on lead 
(0.015 cm, 9 and 13.5 Mev) the results are 
reasonably free from the effects of diffusion so 
that the formula for the path length correction 
may be used. Although this statement is based 
upon the calculated values of @nax this seems 
safe, since the experimental values have been 
lower than the calculated, in cases in which both 
are available. 

(6) The fact that the observed values of @nax 
are smaller than the theoretical is worth noting, 
in view of the fact that others'*~"’ have recently 
reported a deficiency in scattering, both single 
and multiple. It is true that in the case of 0.0025 
cm Al the average angle was found to be in 
accord with Williams’ formula. This, however, is 
only a single case, and there is the possibility 
that it is in error. 

(7) Single scattering could not be checked 
against theory for any of the cases because the 
thickness of the foils and the energy of the elec- 
trons were such that only those electrons scat- 
tered at very large angles would satisfy the 
requirements. In the case of the thinnest ab- 
sorber (0.0025 cm Al, 0.9 Mev) Wentzel’s 
criterion requires that only those electrons 
scattered through ~20 degrees or more be con- 
sidered. Actually no electrons at all were found 
beyond this angle, for the 0.0025 cm Al. 


Energy loss 


In the light of the foregoing data, we may now 
make some statements about the interpretation 
of the measurements on energy loss. 

(1) When we use an absorber of high atomic 
number which is thick enough to cause an easily 
measurable loss of energy, the effect of multiple 
scattering is usually very large. It appears in 
Tables V and VI that the only values for energy 
loss in lead which can be accepted (on the basis 
of scattering) are those for 0.015 cm thickness and 
9 and 13.5 Mev. These are measurements by 


14 A. Barber and F. C. Champion, Proc. Roy. Soc. A168, 


159 (1938). 

1% F.C. Champion and A. Barber, Phys. Rev. 55, 111 
(1939). 

16 W, A. Fowler, Phys. Rev. 54, 773 (1938). 

17 N. L, Oleson, K. T. Chao, J. Halpern and H. R. Crane, 
Phys. Rev. 56, 482 (1939); Phys. Rev. 56, 1171 (1939). 





Fowler and Oppenheimer. The value for 0.038 
cm thickness and 12.7 Mev, by Ruhlig and 
Crane is on the borderline (calculated @msax=18 
degrees). In all the other measurements in lead 
the path in the absorber is so crooked that its 
average length cannot be estimated with reason- 
able accuracy. If we are willing to assume that 
positrons behave exactly as electrons, both in 
energy loss and in scattering, we may find some 
acceptable measurements on positrons in the 
papers already mentioned.*: * These values do 
not change our conclusions, but only support 
what has been said in regard to electrons. 

(2) The number of electrons which strike the 
lead absorber and do not emerge at all is found 
to be large in all those cases in which @max is 
large. This is understandable if we believe that 
in these cases the electrons penetrate the absorber 
by a process very much like diffusion. Those 
which fail to emerge are simply those which have 
traveled a distance equal to several times the 
thickness of the absorber and have thus ‘‘died”’ 
inside the material. The fraction of the electrons 
which fails to emerge from absorbers of various 
kinds has been measured and reported." 

(3) The fact that in some cases in which the 
scattering would be expected to be extremely 
great (for example 0.35 to 0.65 Mev in 0.0066 cm 
lead)* the specific energy loss becomes almost as 
low as the theoretical value, is understandable. 
The thickness of the absorber is not much less 


TABLE VI. Summary of results from other papers on energy 
loss of high speed electrons. 











PATH 
LENGTH 
Tuick- ENERGY 
TuickK- AVERAGE NESS LOSS 
NESS, ENERGY MAX RATIO, EXP./ REFER- 
MATERIAL CM MEV CALC. CALC. THOR. ENCE 
Pb 0.038 9.5 22 — 2 3 
Pb 0.038 12.7 18 1.05 1.3 3 
Pb 0.014 2.5 >25 -— 2.7 5 
Pb 0.015 9.0 15 1,04 1.7 6 
Pb 0.015 13.5 11 1.03 1.5 6 
Pb 0.05 10.0 24 — 1.5 2 
c 0.5 3.0 >25 -- ~1 1 
Cc 0.5 5.0 20 1.07 ~l1 1 
i 0.058 2.5 11 1.03 0.5 5 








18H. R. Crane and J. Halpern, Phys. Rev. 55, 838 (1939). 
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than the maximum range of the electrons in the 
material. Consequently only those which happen 
to have a fairly straight path will emerge and will 
be counted. Those whose paths are more crooked 
will not emerge at all and will not enter into the 
average. Thus there is a selection in favor of the 
electrons which are least scattered. This probably 
accounts for the fact that the energy loss for 
carbon, 0.5 cm, 3 Mev (see Table VI) is about 
equal to the theoretical, in spite of the fact that 
diffusion clearly exists. In this case 65 percent 
of the electrons failed to emerge from the 
absorber. 

After rejecting such a large part of the existing 
work on energy loss, it seems as though there 
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remains little upon which to base an opinion as 
to whether or not theory is in accord with experi- 
ment. If, in the three most acceptable measure- 
ments on lead, we apply the indicated path 
length corrections we find that, as an average, 
the experimental values are at least 1.4 times the 
theoretical. Aside from the question of whether 
or not one is inclined to place any faith even in 
this result, we believe that we have been able to 
expose some of the pitfalls in energy loss measure- 
ment, so that future experiments can be made 
upon a somewhat firmer basis. 

The authors are grateful for the financial 
support of this work which was made available 
from the Horace H. Rackham Fund. 
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The Spin of Carbon Thirteen 


C. H. Townes anp W. R. SMYTHE 
California Institute of Technology, Pasadena, California 


(Received October 16, 1939) 


Carbon thirteen was concentrated in a Hertz diffusion system to an abundance of 50 percent 

and enough 35 percent heavy carbon was obtained for a measurement of the nuclear spin of C®. 

The spin was determined from the relative intensities of the A-type doublets of two lines of the 

O—O Swan band of C“—C"*, These doublets are so close that a Lummer-Gehrke plate crossed 

with a 21-ft. grating was needed to obtain sufficient resolving power. An intensity analysis of 

the combined interference patterns of the two doublets proves that the C" nucleus obeys the 

Fermi-Dirac statistics and strongly indicates a spin of 3. This is in disagreement with a spin of 

} predicted for this nucleus from the Hartree nuclear model but is the value predicted by the 
alpha-particle model. 

HE determination of the nuclear spin of the CONCENTRATION OF CARBON THIRTEEN 

ground state of C™ has become of con- 

siderable interest lately because it is one of the 

few magnitudes concerning which the predictions 

from the alpha-particle nuclear model and from 

the Hartree model are in complete disagreement.! 

The former predicts a spin of $ and the latter a 

spin of 3. The most straightforward method of 

determining the spin is by measuring the relative 

intensities of alternate lines in the CC" molecu- 

lar spectrum. In ordinary carbon, the C"C" 

bands are only one ten-thousandth as intense as 

the C"C"” bands; so that for a measurement of 

the C™ spin, it is necessary to concentrate the 


C* isotope considerably. 
1R. G. Sachs, Phys. Rev. 55, 825 (1939), 


The enrichment of C™ was carried out in a 
34-member Hertz diffusion apparatus described 
by Wooldridge and Smythe. To obtain high 
concentration the diffusion was in two stages. 
For the first diffusion, pure methane was pumped 
slowly and continuously through the light end 
at 1.4 cm pressure to maintain the normal 
concentration of C™ at this end. Six liters of 
methane at 8 mm pressure and about 7 percent 
C® could be collected at the heavy end after a 
five-day run. Due to the large quantity of gas 
passed through the light end, considerable 
amounts of heavy impurities also collected in 


2—D E Wooldridge and W R Smythe, Phys Rev 50, 
233 (1936). 
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the heavy end. The CO and COs: were removed 
continuously by hot I,O; and a liquid-air trap. 
Nitrogen was removed from the heavy methane 
by a simple fractionation. Four five-day runs 
were made to collect enough methane enriched 
in C® to refill the diffusion system for a final 
enrichment. The second diffusion yielded a small 
quantity of methane containing 50 percent C™ 
and about one liter of methane at 1 cm pressure 
of 35-40 percent C". 


SPECTROSCOPIC PROCEDURE 


The nuclear spin and statistics may be deter- 
mined from the band spectrum of a homonuclear 
molecule by measuring the ratio of intensities of 
alternate lines, the statistics determining which 
lines are stronger, and the relative intensities 
being given by s/(s+1), where s is the spin. 
Thus for s=0, as in C®, alternate lines are 
entirely missing, while for spins of 3, 1, and 
the intensity ratios are 3, 3, and 3, respectively. 
Unfortunately, the only bands available in the 
case of carbon are 'II—'II and *II—*II, so that 
the lines whose relative intensities must be meas- 
ured are very close A-type doublets. The Des- 
landres-d’Azambuja (‘II—'II) bands are usually 
excited by the continuous breakdown of hydro- 
carbon gases at fairly high pressures in an 
electrical discharge, and since no method of 
exciting them which used a reasonably small 
quantity of carbon could be found, the Swan 
bands (*II—*II) were used. Each member of a 
single branch in a single band of this system 
consists of six lines. Near the head there are 
three well separated lines, each of which is a 
narrow doublet, the relative intensity of the two 
doublet components being the desired datum. 
Far from the head the triplet structure has 
drawn together into a single line and the doublet 
structure has spread out so that the lines appear 
double. In the intermediate region a combina- 
tion situation holds. The only place where the 
doublet structure is large enough to resolve and 
where the desired lines are not interfered with 
by some overlapping band of C¥®C®, CPC® or 
CXC seemed to be in the (0O—0) band (A5165) 
just before it reaches the (1—1) head. Here the 
triplets of the R branches of C8C*, CP®C™, and 
CC® lie side by side with little overlapping 
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Fic. 1. Spectrogram of carbon bands. 


except that one member, Ri, of each molecule is 
interfered with by the P branch of the next 
lighter molecule (see Fig. 1, a and d). Because 
the carbon molecule approximates Hund’s case 
“a” coupling for low rotational quantum numbers, 
the triplet members *IIy—*IIo(R3) show con- 
siderably wider separation of the A-type doublets 
than do the *II,—*II,(Re) and *Ile—*Ie(R;) the 
separation of the *ITp—*IT lines in question being 
about 0.14 cm. The resolving power required to 
resolve these doublets, assuming they are very 
sharp and have equal intensities is about 140,000. 
Actually the temperature broadening and _ in- 
equality of intensity of the two components 
demanded much higher resolving power. 

As no grating of sufficient resolving power was 
available, a 21-ft. grating was crossed with a 
Lummer-Gehrke plate having a theoretical 
solving power of 300,000. Through the kindness 
of Dr. Badger of the Chemistry Department an 
Kagle-mounted concave grating very fast in the 
second order of the green was used. The some- 
what complicated optical system required by 
the Lummer plate and for correcting the astigma- 
tism of the grating was accomplished by suitable 
lenses supplied by Dr. Bowen. 

Such an optical system requires a powerful 
light source. A discharge tube with a water- 
cooled quartz capillary 1 mm in diameter and 
one cm long was constructed. The heavy methane 
at about 3 mm pressure was admitted to the 
tube and a discharge run until most of the 
carbon was deposited in the tube, principally on 
the electrodes and in the capillary. The tube was 
then evacuated and 2 cm of pure argon admitted. 
When a discharge was run through this tube 
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from a 20,000-volt, 3-kw, transformer with bad 
voltage regulation, the Swan bands appeared 
brilliant for hours. To secure this long life from 
the small quantities of carbon used, it seemed 
necessary to exclude mercury vapor and to have 
the aluminum electrodes polished and with 
rounded edges. | 

Several plates of varying exposures were taken 
with ordinary carbon and with 35 percent C®, 
and step weakener photographs developed with 
each plate to make possible accurate intensity 
measurements. In Fig. 1 6 and c were taken on 
Eastman I-J plates from ordinary carbon and 
35 percent C® with exposures 7 min.—60 min., 
respectively. The *IIp—*IIp doublets of C?—C*, 
which have equal intensity, are not quite re- 
solved by the Lummer plate, and the correspond- 
ing interference fringes of C“°C™ may be seen 
to be somewhat asymmetrical, indicating com- 
ponents of unequal intensity. 


ANALYsIS OF SPECTRA 


There are three lines of measurable doubling, 
R;3(17), R3(18) and R;(19) which, for all the 
molecules, fall clear of other lines. From the 
wave-lengths of the C', lines given by Shea,’ 
the staggering of alternate lines and hence the 
doublet separation, can be computed roughly. 
This calculation indicates about 30 percent 
greater doublet separation for R3(17) and for 
R;(19) tha: for R3(18). This was confirmed 
qualitatively by a visual examination of the 
fringes and so only the first two were analyzed. 

A simple method of determining the statistics 
obeyed by the C™ nucleus is to compare the 
position of the center of gravity of an Rs; line 
relative to the corresponding Rs; line in the 
CP, C8C® and C, bands. Visual measurements 
of the intensity maxima on a comparator showed 
at once that the strong component in the C™, 
line was the one which is missing in the C, 
line, proving that the C™® nucleus obeys the 
Fermi-Dirac statistics. The position of this 
maximum was visually estimated to be very 
close to the position of a single component 
which indicated a very weak second component 
and a probable spin of 3. Subsequent micro- 
photometer measurements, however, make this 


3J. D. Shea, Phys. Rev. 30, 825 (1927). 
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conclusion very doubtful and indicate that it 
may have been an optical illusion caused by 
the tendency of the eye to judge density by 
contrast and thus set the cross hair too near the 
steep edge of the fringe. 

If we are given accurate intensity curves on 
the same horizontal scale of the R3(17) or R3(19) 
lines of C™,, C?C™ and C™, the most precise 
analytical procedure is as follows. Two equal C', 
profiles are drawn with various separations and 
summed. By adjusting the vertical scales of the 
resultant curves only one can be made to co- 
incide with the C"C*® profile which gives, quite 
accurately, the doublet separation. With this 
separation we now draw two CC", curves whose 
intensities are in the ratio 4, 3, 2, etc., and add 
them. The vertical scale of the resultant curves 
are adjusted to fit as closely as possible the C', 
profile. Only one can be made to fit and this 
gives the ratio of intensities and hence the spin. 

To obtain these intensity curves two plates 
were available, one being about three times as 
heavily exposed as the other. On both plates 
there was a uniform change in intensity across 
the plate, the intensity ratio of the outside 
orders being about 2 to 1. On the lighter plate 
eight or nine orders of each line of C". and 
CC® were of suitable density for microphotom- 
etering. On the heavier plate eight orders each 
of the C, lines were suitable but only from 3 to 6 
each of the C', and C"C® lines, the others being 
too dense. Microphotometer curves were run on 
all orders of the lines involved on both plates, 
the vertical height averaging about 2.5 cm and 
the width from 2.5 to 5 cm. From the appropriate 
step weakener picture appearing on each plate 
the line profile was converted, point by point, 
into an intensity profile. The resultant curves 
present a sawtooth appearance due to the grain 
of the plate, the average grain size being about 
1/15 the width of the line. To determine how 
completely this grain would average out a micro- 
photometer curve on the same scale was run 
along one step of the step weakener picture and 
several such traces were superimposed at random 
and averaged point by point. The average of 
eight such curves deviated from a straight line 
by a maximum of 1 mm in a 3-cm length and 
crossed a ruled line 20 times in this distance, 
indicating that the average of eight orders 
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Fic. 2. Microphotometer curves of carbon lines. 


should eliminate the grain structure almost 
completely. This was borne out by the fact that 
all average curves appeared smooth without 
artificial smoothing except such as inadvertently 
occurred by being unable to follow microscopic 
irregularities through the tracing paper. 

The unsmoothed means of 7 orders of R;(18) 
C®, and 8 orders of R3(17) CC", from the light 
plate are shown superimposed in Fig. 2a. It is 
evident that these two lines have the same shape 
and we will therefore be justified in using the 
mean of the same lines on the heavy plate to 
find the shape of the line to be used in analyzing 
the C; lines on that plate. From two curves of 
this shape with a separation of 1.30 cm we obtain 
the solid curve, which is shown in Fig. 2), super- 
imposed on the observed dotted curve which is 
the mean of six orders of R3(19) C"C® from the 
light plate. A 1.25-cm separation fits R;(17) 
CC® better on the sides but not as well on 
the top. 

The same shape for a single line cannot be 
assumed for the heavily exposed plate since 





changes of temperature of the Lummer plate 
might have produced a shift in position during 
exposure. A mean of 3 orders of R3(17) and four 
orders of R;(18) was therefore used for the 
standard profile on this plate. This turned out to 
coincide exactly with that of the light plate. 
The doublet separation was checked on the 3 
orders of R3(17) C"’C™ and 4 orders of R;(19) 
CC® available on this plate and gave 1.20 and 
1.25, respectively, the irregularities being notice- 
ably greater near the tops of the curves as might 
be expected from the smaller number averaged. 
The shapes of the peaks to be expected on this 
plate for R(19) C™, with spins of 3, 3 and 5/2 
were then constructed by adding peaks of the 
standard line shape spaced 1.30 cm apart with 
intensity ratios 1 to 3, 3 to 5 and 5 to 7, respec- 
tively. The resultant curves with the vertical 
scale suitably adjusted are shown superimposed 
on the observed dotted mean of 8 orders of 
R;(19) C™, in Fig. 2c. Corresponding curves for 
R;(17), using a 1.25-cm spacing, are shown in 
Fig. 2d. All curves fit almost exactly for two- 
thirds of the distance up the steep side. On the 
side of the weak component the $ curve fits well 
on the average but the bump which should be 
present is absent. No known smoothing process, 
such as the finite microphotometer slit width, 
could eliminate this bump completely. Acci- 
dental deviations of this order have been ob- 
served as shown in Fig. 2), but would be un- 
likely to occur in both lines. On the other hand, 
the observed curve is much too wide to agree 
with the calculated } curve. No accidental 
width deviations of this size have been observed 
even between single unaveraged curves. The 
standard shape used for the single line was well 
determined so the use of a wider standard cannot 
be justified. Use of a wider separation improves 
the } fit at the base but makes it worse in the 
middle. The spin of 1 is excluded by the 
statistics. 
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Wave Reflections from Diffuse Boundaries 


C. D. THoMAs AND R. C. COLWELL 
West Virginia University, Morgantown, West Virginia 
(Received September 18, 1939) 


The intensity of the waves reflected from various types of diffuse boundaries in which the 
dielectric constant changes only slightly is considered. The constant inside the reflection layer 
is assumed to vary as some power of the distance. It is further assumed that above and below 
this transition layer the media are uniform. The calculations show that the rate of decrease of 
intensity of the reflected wave with increasing layer thickness is more rapid for large values of 
uw the index of refraction, than for values close to unity. Fresnel’s equations are therefore 
applicable to layers of considerable thickness if u is very close to unity. 


T is customary to consider reflections from a 

boundary to be of considerable magnitude 
only if the boundary is very sharp compared to 
the wave-length of the incident wave. This is 
the case of importance in optics where the 
reflected light is usually near the order of 
magnitude of the incident ray. In the reflection 
of radio waves, however, the sensitivity of the 
receiver permits the detection of comparatively 
feeble waves, and it has been shown that de- 
tectable reflections may be produced by non-ionic 
layers if the boundaries are considered sharp." 
It does not seem probable, on the other hand, 
that the layers could have extremely sharp 
boundaries, relative even to radio wave-lengths. 
It is therefore desirable to consider whether or 
not one is justified in using Fresnel’s reflection 
laws here, as has sometimes been done. The 
important difference between these radio re- 
flections and those of light is that the radio 
reflections are from regions where the relative 
refractive index is very close to unity so that 
the coefficient must. be very small even though 
the boundary is truly sharp. 

For the purpose of testing whether or not 
Fresnel’s laws are a good approximation for 
thicker transition layers as the index approaches 
unity, a one-dimensional problem will be con- 
sidered, i.e., plane waves progressing in the 
direction of variation of the dielectric constant. 
The region near the origin will be taken as 
uniform up to some point x, above this the 
dielectric constant will be assumed to vary as 


1C. R. Englund, A. B. Crawford and W. W. Mumford, 
Bell Sys. Tech. J. 14, 369-387 (1935). 


some monotonic function of x up to x2, above 
this it is to be again uniform. 

The wave equation, after removing the time 
dependent part, becomes 


PX /dx?= —veX/c, (1) 


where X is the space dependent part of the 
electric vector, c is the velocity of light, and « is 
the dielectric constant which is itself a function 
of x. The general solution of this equation is 


X = Af(x)+Bg(x). (2) 


The boundary conditions which must be satis- 
fied by the solution are that x and its first 
derivative for the lower region must equal those 
of the transition layer at its lower boundary, 
and a similar set relating the upper region to 
the upper boundary of the layer. Thus there are 
four conditions for the six constants A, B (Ao, Bo 
in the lower space; A, B; in the layer; Ae, Be in 
the upper space). However, it is only the ratio 
Ao/By that is required if f(x) and g(x) are made 
the reflected and incident waves, respectively. 
Moreover, in the upper space one may assume 
that the backward component of the wave is 
zero (A2=0). Hence, there are left only four 
constants (Ao/Bo, Ai, Bi, Be), and there are 
four boundary conditions for determining these. 
The solutions of (1) in the three regions are now: 


X o=Ay’ exp (tke;'x)+Bo’ exp (—tke:'x), 
X,=Aif(x)+Big(x), (3) 
X2= By’ exp (—tkee!x), 


where k= v/c=27 frequency/c. For convenience, 
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let 


A’ exp (ike:'x)+Bo’ exp (—tke,'x) 
=A, exp [tke}(x—x) ] (4) 
+B exp [ —ike,}(x—x) ]. 


The boundary conditions then become 


Aot+Bo=Aif(xi1)+Big(x:), 5) 


ike} (Ao— By) = Af’ (x1) + Big’ (x). 


And similarly 
Bo=Aif(x2)+Big(xe), 
—ike' Be = A if’ (x2) + Big’ (xe). 


The primes on f and g of (5) and (6) indicate 
differentiation with respect to x before evaluation 
at x; and Xe. 

Upon solving Eqs. (5) and (6) for the re- 
flection coefficient when the thickness is ¢ 


(6) 


A Lf’ (xr) + ther f(x) ) JLe’ (x2) +tkes'g(x2) ]— Cf’ (x2) +ikeef(x2) Jeg’ (x1) +ike: *g(x1) J 


————. (7) 





~ Bo —Lf'(x1) —iker¥f(x1) Ie’ (x2) + ikeste( x2) +L’ (xa) + ikeatf(xs) Ie’ (xi) —ikes4g(x1)] 


This is applicable to any sort of transition layer 
of any thickness and for any value of the relative 
index w= (€2/e;)?. 
More specifically suppose that Eq. (1) inside 
the layer is 
d*X ; 
-+ F(x)X =0, (8) 
dx? 
where 
x” yp x” 


C F(x) }}=-e1! =-¢,!—. (9) 


¢ zor ¢ 





This is not a very serious restriction on F(x) 
when it is recalled that the layer is not very 
thick and that m may have any value whatsoever 
either positive or negative. Since x lies between 
x; and xe, it is evident that if (€2/e;)! is greater 
than unity, 2 must be positive, while if (€2/¢,:)! is 
less than unity, 2 must be negative. In general 
from Eq. (9) 

a= * (€2/e1)*= (%2/21)" _ (10) 


; 1k 
in| exp| 
n+1 





R.= 


(atrs—an) | 


or using 





Xe—X,=t, x,=t/(u"—1), (11) 


where ¢ is the thickness of the transition layer. 

In the cases under consideration yu is greater 
than one by a very small amount so that in 
Eq. (11), x: and xe are large except for layers 
which are extremely thin. Under these circum- 
stances, it is allowable to put 


f=e4, goe iD+5 (12) 
where 
ke,! yntl n 
A=——_——,, = —-— log », (13) 
n+1 (x,)” 2 


as may be verified by substitution into Eq. (8) 
with the aid of (9). 

If these values of f and g are put into Eq. (7), 
the coefficient R, takes the form 


1 x1 —1k 
x | (a's etn) || 
Me Xe n+1 


—, (14) 


——— CXp 








1 
4ke,'x, exp 
n 


If use is now made of the value of x; and xe 
from Eq. (11) along with the restriction that yp is 
always near unity, 


| ' utn 
n?(u'/™"—1)? sin? | ke, ) 
n+1 


R?= . (15) 
4ke,t? 











k 
sastie -_ ale) 


The value of R? reaches 2 zero at approximately 


n+l 
kt=—_ ——. (16) 
et utn 


It may easily be verified from this that the 
value of the thickness ¢ for the first zero is greater 
for a value of w nearer unity, since the factor 
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Fic. 1. The variation of the intensity coefficient with 
layer thickness for two sets of values of €; and es. The index 
of refraction u =(¢€2/e,)4 is nearer unity than wy’. 


(n+1)/(u+m) in Eq. (16) increases for smaller 
values of yw. Also from the value of dR?/dt it 
may be seen that R? decreases less rapidly for 
values of wu nearer unity. 

The result on an exaggerated scale is illus- 
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trated in Fig. 1. Thus, although the reflection is 
much greater for truly sharp boundaries if y 
differs considerably from one, the rate of decrease 
is also much greater. Hence, the Fresnel coeffi- 
cient (t=0) becomes a better approximation for 
a given thickness if the value of yu is closer to 
one. In order to obtain information as to the 
exact way in which R, decreases with ¢, it is 
necessary to make some assumptions as to the 
value of ¢« in the lower and upper regions and 
regarding the manner in which it changes in the 
transition layer (i.e., the nature of F(x) of 
Eq. (8). These calculations will be published 
later. 

The formulas developed here are applicable to 
the reflection of radio waves from different air 
masses in the troposphere. It is known from 
experiment that such reflections take place.’ 


2 R. C. Colwelland A. W. Friend, Nature 137, 782 (1936). 
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On the Interpretation of Atomic Distribution Curves for Liquids 


C. A. Coutson, University College, Dundee, Scotland 
AND 


G. S. RusHBROOKE, H. H. Wills Physical Laboratory, Bristol, England 
(Received July 24, 1939) 


It is shown that if we approximate to a liquid by an Einstein model, in which each atom has a 
restricted region of motion, wherein it moves independently of its neighbors, and is surrounded 
by coordination shells of other atoms; and if we denote the density distribution of the atoms by 
p(r), where r is the radial distance from any giver. atom; then the contribution to p(r) made by 
any coordination shell, e.g., the zth is a function p;(r) for which rp;(r) is symmetrical about its 
corresponding maximum value. The complete distribution curve, rp(r) against r, is the sum of 
peaks of equal width and similar shape. A semi-empirical application of this theory to liquid 
sodium (a reapplication of C. N. Wall’s theory) gives a latent heat of melting in fair agreement 
with experiment. The model suggests a change of structure on melting, since agreement with the 
experimental distribution curve is impossible if the number of atoms in the first coordination 
shell is that of solid sodium. It is also shown that if the parameters in a partition function 
developed in this way are chosen to give agreement with any one physical property, then it is 
incorrect to add to the partition function terms representing ‘‘communal entropy.” 


VOLUME 


§1. INTRODUCTION 


EVERAL experimental papers have appeared 
recently publishing atomic distribution 
curves for liquids.' These curves, found from 
1(a) F. H. Trimble and N. S. Gingrich, Phys. Rev. 53, 


278 (1938) ; (b) C. D. Thomas and N. S. Gingrich, J. Chem. 
Phys. 6, 411 (1938); (c) C. D. Thomas and N. S. Gingrich, 


Fourier analyses of the intensity distributions in 
X-ray scattering photographs, show p(r), or 
4rr*p(r), as a function of r. 42rr’p(r)dr is propor- 
tional to the probability that two unspecified 
atoms of the liquid are distant r to r+dr apart. 


J. Chem. Phys. 6, 659 (1938); (d) J. Morgan and B. E. 
Warren, J. Chem. Phys. 6, 666 (1938). 
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- Fic. 1. (Left) Typical distribution curves for a liquid. (a) No structural units 
; present. (b) With structural units. 

Fic. 2. (Right) (a) Cybotactic groups in a liquid. (b) Structural units in liquid 
yellow phosphorus. From the experimental value of the density (found from the 
parabolas in Figs. 4 and 5 of Thomas and Gingrich’s paper), i.e., po! = 29.1 A’ at 

48°C, and =31.4 A’ at 226°C, we can calculate an average distance apart for the 
centroids of two neighboring 'P, units. Since the data of Thomas and Gingrich 
suggest that each unit has a coordination of about 10, we have taken the mean of 
the results, which differ by only about 2 percent, found on the assumptions of 
body-centered and face-centered cubic packing. Thus the length OO’ is found to be 
5.42A at 48°C, and 5.55A at 226°C. It is then clear that the first, second and third 
peaks in the figures of Thomas and Gingrich refer to the mean values of PQ, OR 
and PR. On account of valence saturation the atoms Q and R try to get as far 
apart as possible; this is why the distances AA’, which are the shortest possible 
geometrically, do not occur in appreciable numbers. 


In each case curves appropriate to various tem-_ relation of this structure to the latent heat of 
peratures are given. melting (L=T7AS) which occurs at the phase 

It is important to consider with what other change; incidentally this entails enquiry into 
experimentally determinable properties of the the circumstances in which it is relevant to 
liquid we can correlate these curves; we may use_ introduce the phrase ‘“‘communal entropy.”’ We 
either the shape of a single curve, the effect of a__ shall be concerned mainly with the case of liquid 
change of temperature on the curve, or relations sodium, for which C. N. Wall? has already given 
between distribution curves for a family of sub- a theoretical discussion; indeed one object of 
stances such as the alkali metals. And the proper-_ this note is to reconsider Wall’s application of 
ties we deal with may be physical ones, e.g., his theory, for it seems to us that this can be 





density, viscosity, electrical resistance; or ther- 
modynamic ones, e.g., entropy, vapor pressure, 
specific heat. 

In this note we shall discuss the bearing of 
atomic distribution curves on the geometry of 
the atomic structure of a liquid, particularly in 
the neighborhood of the melting point, and the 





done in a more satisfactory way. 


§2. THE MEANING OF PEAKS IN A 
DISTRIBUTION CURVE 


A typical distribution curve for a liquid, 
4rr*p(r) plotted against r, has peaks as shown in 


?C. N. Wall, Phys, Rev. 54, 1062 (1938). 
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Fig. 1(a). The question arises whether these 
peaks have any meaning as far as the geometry 
of the structure of the liquid is concerned. 

At temperatures just above the melting point 
the peaks are sufficiently sharp and charac- 
teristically shaped for it to-be generally agreed 
that liquids have some sort of semi-crystalline 
structure.’ In particular two types of atomic 
arrangement have been proposed; one_ being 
characterized by molecular homogeneity, i.e., 
all atoms are thought to be equivalent in the 
structure of the liquid, and the other by the 
presence of cybotactic groups,’ i.e., it is supposed 
that there are groups of atoms, probably of 
varying sizes and shapes, which have a more 
definite crystalline structure within them than 
the liquid has as a whole. 

It is not, however, possible to decide, either 
from a single distribution curve, or from several 
curves for different temperatures, whether a 
liquid is molecularly homogeneous or consists 
mainly of small cybotactic groups, each con- 
taining less than, say, 2000 atoms: Fig. 2(a). 
For in the latter case the x-ray diffraction lines 
will be very much blurred by edge effects at the 
boundaries of the groups, irregularly spaced 
atoms between the groups and the unequalness 
of the spacings within the groups, due to their 
small linear dimensions. This blurring will make 
the distribution curve indistinguishable from 
that to be expected with a molecularly homo- 
geneous liquid. 

On the other hand, we could see from a dis- 
tribution curve whether there were only large 
cybotactic structures, each containing more than 
about 2000 atoms;> for then the peaks would be 
quite sharp and separate from each other, 
approaching the lines we should get if we 
analyzed the intensity distribution in a powder 
photograph. We can also tell whether the atoms 
of the liquid form smaller groups, which we may 
call structural units, as, for example, with liquid 
yellow phosphorus in which the atoms are 


3(a) P. Debye and H. Menke, Ergebnisse d. Tech. 
Réntgenkunde II, 1 (1931); (b) N. F. Mott and R. W. 
Gurney, Reports on Progress in Physics (Physical Society 
of London) 5, 46 (1938). 

‘J. D. Bernal, Trans. Faraday Soc. 33, 27 (1937). 

5 It is interesting in this connection that Riley, Chemistry 
and Industry 58, 391 (1939) has been able to determine by 
use of x-rays the mean size of small graphite crystallites, 
which have about this limiting size. 
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grouped, both in the liquid and in the vapor, as 
P, molecules. In this case the fact that the first 
peak is detached, p(r) vanishing, or sensibly 
vanishing, between the first and second peaks, 
shows that there are structural units com- 
paratively far apart from each other; Figs. 1(b) 
and 2(b). We shall indeed have valence satura- 
tion in each structural unit and therefore repul- 
sion of any atom not in the same unit; see note 
under Fig. 2(b). The breadth of the first peak 
will be due to fluctuations in the atomic distances 
within the P, groups, as a result of vibrations, 
rotations and occasional collisions. If only a 
small fraction of the atoms are grouped in 
structural units, then the detached peak due to 
them will be masked and not separate from the 
rest of the atomic distribution curve. 

On the assumption of molecular homogeneity 
the function 47r*p(r), which gives the distribu- 
tion of interatomic distances in the liquid, will 
give the distribution of the distances of the 
atoms from any one particular atom. If, for 
clarity, we suppose that there are no small 
structural units, then according to Bernal* we 
can define round any atom first, second... 
coordination shells as in a crystalline solid except 
that now the atoms in these are distributed, 
somewhat randomly, about mean distances from 
the central atom. For the shells to be properly 
defined a suitable probability function for the 
distribution of the radial distances about the 
mean radii of the coordination shells ought to be 
prescribed. Bernal, following earlier work of 
Prins,* assumes a Gaussian distribution, so 
that in a shell of mean radius Ri, pi(r)a 
x<exp [—A(r—R;)? ]. On the other hand, Wall 
finds that, to his approximation,? rp;(r) and not 
pi(r) or even r’p;(r) is, for each shell, symmetrical 
about the mean radius. There seems to us to be 
insufficient experimental data for deciding be- 
tween these three possibilities ; we shall, however, 
show theoretically in the following paragraph 
that for any Einstein model of a liquid, rp;(r) for 
a single shell, is symmetrical about the mean 
radius of that shell. Wall’s theory, which we are 
going to reapply to liquid sodium, is based on a 
particularly simple Einstein model. 

An Einstein model is characterized by the 
assumption that each of the atoms of a substance 

6 J. A. Prins, Physica 3, 147 (1936). 
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Fic. 3. Mean positions O and O’ of the atoms P and Q. 


moves independently of its neighbors, vibrating 
to and fro about a fixed center in some kind of 
average field which is determined by the rest of 
the substance. (See papers by Lennard-Jones and 
Devonshire.’?) Thus, in Fig. 3, O and O’ are the 
mean positions of two atoms P and Q. Ri, i.e., 
the length OO’, is then the mean radius of the 
shell in which Q is found relative to P. We denote 
the average potential field in which each atom 
is assumed to move independently of its neigh- 
bors, by ¢(s), and suppose that it is spherically 
symmetrical. Then the probability that P lies ina 
small volume element d Vj, distant s =r, from O, is 
proportional to exp [—¢$(r1)/RkT ]dVi; similarly 
the probability that Q lies in dV2, distant s=re 
from O’, is proportional to exp [ — ¢(re)/RT ]d V2. 
And in fact, the joint probability that P lies in 
dV, and Q in dV 2 is 


1 
—e—lo(ri)+o(r2)1/kTd Vid Vo, (1) 


v* 
where 


y= fever MdVi= fee kTd Vo. 


The probability that the length PQ lies between 
r and r+dr is the integral of (1) taken over all 
those coordinates of P and Q for which r& | PQ) 
€r+dr: thus, if p;’(r) denotes the density con- 
tribution per atom in the ith coordination shell, 
it follows that 


1 
dertpi(r)dr=— f e-ternssenn ‘Td V dV», 


ud 


v* 
subject to 
r&|PQ| <r+dr. 
Now if there are N; atoms in the ith shell, 


7 J. E. Lennard-Jones and A. F. Devonshire, Proc. Roy. 
Soc. A163, 53 (1937); A165, 1 (1938). 
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pi(r) = Nip,’ (r): further, if we denote the length 
O’P by x, and angles about the lines R; and x by 
x1 and x2, then we may write 


dV,=xr,dxdridx;/Ri, dV2=rredrdredx2/x. 


Introducing these results into (1), and taking the 
integration over dV, to precede that over dV, 
we find that 


rN; x Rit+ri r+z 
ro(r)=——f f f rye o'r) kT 
9 
Rw? =|Ri-—r,| 


ry=0"% x ra=|r—z| 


“roe? ("2/kT. drodxdr). 


Now owing to the shape of $(s)—see Fig. 1 of 
Lennard-Jones and Devonshire’7—we can assume 
that there is a length o such that ¢(s) is effectively 
infinite for s2o. The integration over re gives a 
factor 


¥(r+x)—y(|r—x]), 
where ¥(t) = /o's exp [—$(s)/RkT lds. Now r+x 


is necessarily greater than o, and clearly y(t) 
=y(c) if t20, so this factor is equal to 


¥(0)—¥(|r—x]). 


For the integration over 7; also, only the range 
from 7;=0 to r;=¢ is significant, and so we find 


rN; a 
rp;i(r) =— | reo (rier 
0 


iv" 
Ritri 


x| arto) - v( r—x|)defars (2) 


Ri-ri 


This is a function of |y—R;| and vanishes if 
|r—R;| 2 2c. It therefore represents a contribu- 
tion to the total function rp(r), in the form of a 
symmetrical peak with a spread at the base given 
by 4c, i.e., independent of Rj, and therefore the 
same for each coordination shell. Wall’s Eq. (10) 
is a particular example of our equation for the 
case in which ¢(s) is constant for s<o. 
According then to Wall’s theory, or to any 
theory based on an Einstein model for a liquid, 
4mrp(r) is resoluble into peaks of equal width 
(4c) and similar shape. We give below, Fig. 4(a), 
the result of an empirical attempt which we have 
made to separate into symmetrical peaks the 
curve for 4xrp(r) which pertains to the experi- 
mental data of Trimble and Gingrich" for liquid 
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Fic. 4. Distribution curves for liquid sodium at 100°C. 
(a) Empirical separation into symmetrical peaks. (b) 
Comparison of predicted distribution curves with the 
experimental curve. 


sodium at 100°C, and we see that the peaks are 
approximately of equal width and similar in 
shape. This is as satisfactory as we can reason- 
ably expect, bearing in mind that the peaks are 
successively more sensitive to the precise shape 
of the first peak, i.e., to the position of its 
maximum and its form for small values of 7, for 
which we cannot be certain of the experimental 
values. 

The first two or three peaks do not overlap 
much and therefore may be supposed to cor- 
respond to coordination shells which have some 
physical meaning. But subsequently there is very 
considerable overlapping, and so there is little 
geometrical basis for assigning atoms, at these 
greater distances from the central atom, to 
definite coordination shells about that atom. 
This is even more the case at higher tempera- 
tures, and eventually the peaks can have little 
significance as far as the geometrical structure of 
the liquid is concerned; for whether this is 
interpreted in terms of coordination shells or 
small cybotactic groups, which will subdivide as 
the temperature is raised, it loses its significance 
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as the structure of the liquid approaches that of 
a very imperfect gas. 

In liquid sodium there will be no structural 
units other than atoms (or ions) and possibly a 
small percentage of diatomic molecules if the 
temperature is great enough. At 100°C, ie., 
2.4° above the melting point, we can ignore the 
possibility of the latter. We shall make the 
assumption, which in this case seems to us most - 
reasonable, of molecular homogeneity, and 
interpret the peaks in the 4mrp(r) curve as cor- 
responding to coordination shells. Unfortunately 
the small, 2.5 percent, volume change of sodium 
on melting does not rule out the possibility of 
cybotactic groups in the way that the com- 
paratively large changes of —9 percent and 11 
percent for water and argon, respectively, seem 
to do in those cases: for large volume changes 
must imply a radical change of structure, but 
small volume changes leave the question open. 


§3. REAPPLICATION OF WALL’s THEORY TO 
Lieuip Sopium AT 100°C 


Wall’s theory is based on the assumption that 
for each atom (or ion) of the liquid there is a 
sharply defined, spherical volume of radius ¢, in 
which its center can move without increasing the 
potential energy of the liquid as a whole. This 
implies among other things, that the atoms (or 
ions) behave as hard rigid spheres. On this 
assumption Wall finds a theoretical form for the 
shape of an isolated peak, viz: 











3N; 5(R;-7r)? 
y=4arp,(r) = | i— 
5cR; 4c? 
|Ri-r| |R:—r|? 
fry oy 
20 400° 


This depends on three unknown parameters oc, 
N; and R;. o has been defined already ; N; is the 
number of atoms in the coordination shell 
corresponding to the peak for which R; is the 
mean radius. To get the complete distribution 
curve, 4zrp(r) against 7, it is necessary to add 
together the effects of successive peaks; o will 
be the same for each of them but their N’s and 
R’s are further unknown parameters. Wall deter- 
mines No, Ro and o, where No and Ro appertain 
to the first shell, by making the left-hand side 
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of the first peak, r< Ro, coincide as closely as 
possible with the experimental curve. He then 
uses some semi-empirical calculations of Bernal 
to find, from No and Ro, the N’s and R’s which 
correspond to more distant shells. Knowing 
these, the complete distribution curve can be 
constructed, and Wall gets a curve which agrees 
fairly well, but by no means exactly, with the 
experimental one. Unfortunately, in drawing a 
curve through the experimental points a small, 
but quite definite, peak (strictly speaking the 
second), was smoothed out; this makes the 
agreement Wall has obtained (reference 2, Fig. 2) 
seem rather better than it really is. 

Now the values of ¢ and Np found in this way 
are very sensitive to the precise shape of the first 
peak in the experimental curve. Moreover, we 
must not expect Wall’s theoretical peaks to fit 
the experimental curve exactly, since that would 
imply that the atoms (or ions) behaved as rigid 
spheres. Finally it seemed to us better to try to 
test Bernal’s formulae for N; and R; in terms of 
No and Ro, rather than to assume them in the 
course of this theory. So to find ¢ and No, Ni:-- 
we proceed rather differently. 

An essential feature of Wall’s theory, and 
indeed of any theory based on an Einstein model, 
is that, for any shell, N; is proportional to R; Yi, 
where Y; is the value of y in Eq. (3) which cor- 
responds to r= R;. We find the correct values of 
R; and Yj, for each peak, from our empirical 
resolution of the experimental distribution curve 
into separate peaks. We then determine No and 
N, from the two conditions 


No/Ni=Ro Yo/Ri Y;, (4) 
1+No+N,/2=4/3: x: Ri®po, (5) 


where po is the density of the liquid in atoms per 
unit volume. We use the experimental value for 
sodium at 100°C,® which is pop=0.0247A~*. The 
first formula (4) is essential to Wall’s theory and 
the second (5) is the result of assuming that if we 
take a small imaginary sphere of radius R, then 
it will enclose the same number of atoms 
wherever it be situated in the liquid. This is 
approximately true for a crystalline solid body, 
and we may suppose that it will be even more 
so for a liquid. 


SE, Rinck, Ann. de Chemie 18, 395 (1932). 





We find 
Ro=3.78A, R,=4.96A, 
Yo= 2.45, Y,=1.21, 


and with these values, (4) and (5) then give 
No=8.78, N,=5.69. 


The empirical resolution of the distribution 
curve into peaks, as shown in Fig. 4(a), suggests 
that the third peak is specified by R2=6.18A, 
Y:=1.74. But when we come to use Wall's 
formula to construct the curve theoretically, we 
find that we get a much better fit if we move the 
third peak a little to the left, and have R.=6.06A, 
Y2=1.74. This shift of Re lies within the uncer- 
tainties unavoidable in any graphical separation 
of the peaks. The formula N2/Ni=R2Y2/Ri Vi 
then gives N2=10. The results, are summarized 
in Table I. We give Wall’s values for comparison. 
In Fig. 4(b) we show the distribution curve ob- 
tained on using our values for the parameters in 
Wall’s theoretical formula. o is given by 
3No/5RoYo and is found to be 0.57, whereas 
Wall’s value was 0.60. Wall’s curve is shown by 
the broken line in Fig. 4(b). 

We shall not give the details of a similar 
analysis of Trimble and Gingrich’s experimental 
curve for liquid sodium at 400°C, because we 
do not suppose that, at that temperature, a 
resolution into peaks will have much meaning. 
The calculations give, at 400°C: 


No=8.50, N,=6.81 ° 
Ro= 3.80, R,=5.14. 


TABLE I, Comparison of constants. 








No Ni Ne Ro Ri R: R:i/RoR:/Re 


o 
Values; this note 8.78 5.69 10 3.78 4.96 6.06 1.31 1.60 0.57 
Wall's values 9.28 6 13 3.79 5.1 6.2 1.34 1.63 0.60 











Finally the entropy change on melting (AS) 
and the associated latent heat of fusion (L) are 
estimated as in Wall’s paper; the results are 
shown in Table II. 


§4. DIsCUSSION OF THE RESULTS 


We do not wish to attach any importance to 
the precise numerical values which we have 
found and given above; indeed it would be un- 
reasonable to expect better agreement than that 
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TABLE II. Melting of sodium. 

















[AS] L (KILOJOULES/MOL.) 
Present paper 0.747 R 2.31 
Wall’s paper 0.906 R 2.8 
Experimental 0.861 R 2.65 





shown in Table II from so simple a theory, 
but we believe that certain conclusions can, 
legitimately, be drawn from them. 

In the first place, at 100°C, i.e., only 2.4° 
above the melting point, No differs quite appreci- 
ably from its value, which is 8, for sodium in the 
solid state. The first two shells together contain 
nearly the same number of atoms (14.4) as in 
the solid state (14). We have tried to obtain a 
fit with the experimental curve on the assump- 
tion that No=8, but could not obtain even 
reasonably good agreement. We conclude there- 
fore that there is quite a definite change of 
structure on melting. Moreover No (400°C) is 
less than Np (100°C) while N, (400°C) is greater 
than N, (100°C), so that it appears that at 
higher temperatures the number of atoms in the 
second shell increases at the expense of those in 
the first. The liquid, then, does not tend to 
become close-packed at high temperatures. 

Secondly, on comparing the predictions of 
Bernal’s paper, i.e., that for No=8.78 we shall 
have N,=8, Ne=12, Ri/Ro=1.33 and R2/Ro 
= 1.63, with the results of Table I, we find that 
Bernal’s values for Ri/Ro and Re/Ro are in 
good agreement with those we have found 
empirically; but his values for Ni, Ne --+ are 
not those which we have found. The only agree- 
ment here is that if No, when in the neighborhood 
of 8, decreases, then JN, increases. Bernal himself 
states that his estimates of successive N’s are 
less reliable than those of successive values 
of R/Ro. 

Next we observe that using the same N’s and 
R’s as in Table I, we could have got a better, 
and almost perfect, fit with the experimental 
curve if- Wall’s theoretical shape for a single 
peak had been rather fatter at the base. This 
would have been so if the atoms (or ions) had 
not been treated as hard spheres, with sharply 
defined free volumes, but their repulsive inter- 
actions had been taken into account, as in 
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Lennard-Jones’ picture of a liquid.’ This is an 
essential defect of Wall’s theory. 

Finally, since our value for the latent heat of 
fusion is too low, we are led to enquire whether 
this also is due to the essential limitations of a 
model in which the atoms are treated as hard 
spheres, or whether it is due to our having 
neglected a communal-entropy term (see, in 
particular Lennard-Jones and Devonshire, 1938). 
We are driven to the former conclusion and end 
with a discussion of this point. 


§5. THE MEANING OF COMMUNAL ENTROPY 


The introduction of the notion of a free- 


volume for each atom of a liquid, in which, 
effectively, it can move without influencing the 
other atoms, at once elicits the idea of communal 
entropy.® 

Free volume may be defined in two ways. 
We shall say that it is defined thermodynamically 
if it is defined purely formally from the partition 
function. The partition function F(T, V, N) is 
supposed known (though in practice only ap- 
proximations to it are known) as a function of 
T, V and N, where T is the temperature, V the 
volume and JN the total number of atoms of the 
liquid; the thermodynamic free-volume 2n, is 
then defined by the equation: 


F(T, V, N) =(2armkT/h?)**"*-on, (6) 


where m is the atomic mass, and the other 
symbols have their usual meanings. 

On the other hand, we shall speak of the 
geometrical free-volume Ygeom When we approxi- 
mate to the liquid by an Einstein model. decom is 
then v of Eq. (1) and replaces 1%, in Eq. (6). 
It has a well-defined meaning geometrically in 
the simple approximation in which it is supposed 
that each atom may be treated as a rigid sphere 
imprisoned in a cell whose size is fixed by the 
density of the liquid; this is the approximation 
on which Wall’s theory is based. 

If we restrict ourselves to Einstein models for 
liquids, v4, is an ideal quantity to which we try 
to approximate by assuming suitable forms for 
Vgeom as a function of JT, V and N. vm is an 


® (a) O. K. Rice, J. Chem. Phys. 6, 476 (1938); (b) R. W. 
Gurney and N. F. Mott, J. Chem. Phys. 6, 222 (1938) ; (c) 
L. Tonks, Phys. Rev. 50, 955 (1936). 
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average determined by all the positions of all 
the atoms of the liquid, though only those that 
are at all likely will count significantly. vgeom, in 
which each atom is confined to the neighborhood 
of a certain site, may then be too small, since 
by confining the atoms to particular regions of 
space, large fluctuations in the interatomic 
distances are ruled out; but it need not be too 
small because exact agreement with %, may be 
obtained, at any one temperature, by suitable 
choice of the potential field ¢. Of course it may 
not be possible to get a good approximation 
to v, at all temperatures by choice of a potential 
field ¢, which must depend explicitly only on 
the density and be otherwise independent of the 
temperature. For this reason an approximation to 
the partition function is sometimes made (e.g. 
Lennard-Jones and Devonshire, 1938) in which a 
reasonable form for @ is chosen, and then the 
resulting ‘‘Einstein” partition function is multi- 
plied by another factor, frequently e’, to take 
account of large fluctuations in the atomic 
distances. It is this extra factor which is called 
the communal entropy term and e% is chosen 
because that is the difference which would be 
made to the partition function of a perfect gas 
of N molecules if its volume were to be sub- 
divided into N equal cells and each molecule be 
confined to one cell. But of course the factor 


ought strictly to be temperature dependent and ° 


its nature must depend on the initial choice of ¢. 
Now at any temperature, the measurable 





properties of a liquid depend on all the possible 
positions of all the atoms of the liquid, and to 
predict them theoretically we must use a good 
approximation to vx. If we include parameters 
in the potential field ¢ of Ugeom—for example in 
Wall’s case the radius of his spherical free volume 
—and so choose these that some physical 
properties of the liquid have their correct experi- 
mental values, then we are not at liberty still 
to include a communal entropy term in the 
partition function. If the resulting formulae do 
not give the correct experimental values of other 
measurable properties, this is because we do not 
have sufficient parameters in @¢, i.e., to an 
essential limitation of the model used. 

This is the case in the present instance. We 
have chosen ¢ in Wall’s theoretical formula for an 
atomic distribution curve, Eq. (3), so as to get 
as good a fit as possible with a distribution 
curve found experimentally, and the fact that 
Wall’s theory does not then give a good value 
for the entropy of melting is due, we suggest, 
not to omission of a communal entropy term— 
inclusion of which would mean that there were 
positional distributions of the atoms which 
affected the entropy but not the atomic dis- 
tribution curve—but to the essential limitations 
of a model which treats the atoms (or ions) of 
the liquid as hard spheres. 

The writers would like to acknowledge helpful 
discussion of the ideas in this paper with Mr. J. 
Corner. 
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In evaluating the anisotropy constants of ferromagnetic 
materials a common method is to measure the torque 
exerted on a specimen of circular section by a uniform 
magnetic field, to determine the maximum torque by 
extrapolation to infinite field strength, and to multiply 
this by an appropriate number. The law of approach of 
the torque, 7x, to its limiting value, T.., is now found to be 


Tu =T..(1—Ho/Ha) 


for a number of disks and ellipsoids of single crystals of 


iron-silicon alloys containing up to 7.4 percent silicon, and 
for polycrystalline disks of iron and of silicon steel. Here 
H, is the strength of the applied field uncorrected for 
demagnetizing field and Ho is a constant which depends 
mostly upon the dimension-ratio of the disk. This law is 
discussed in relation to a similar law proposed by 
Schlechtweg, and to the familiar law of approach of mag- 
netization to saturation. It is concluded that for disks the 
approach to saturation of the torque peaks is connected 
with the lack of magnetic saturation at their edges. 





T is a well-established fact that ferromagnetic 
single crystals generally show considerable 

anisotropy in their magnetic properties. The 
simple theory of ferromagnetic anisotropy and 
the several methods of measuring it have recently 
been reviewed by Bozorth.' The simplest of 
these methods involves the use of a magnetic 
torquemeter to measure the torque exerted by a 
strong field on a suitably held single crystal 
disk, the torque arising from a lack of parallelism 
between the direction of the magnetization in- 
tensity in the disk and the direction of the 
applied field. This deviation between J and H 
occurs because of the anisotropy of the material. 
If the orientation of the disk with respect to the 
crystal axes is known, the anisotropy constants 
can easily be calculated from the torque data.?~* 
It is usually preferable to have some simple 
crystallographic plane lying in the plane of the 
disk because the resulting torque curve must 
then be symmetrical, unless the crystal is im- 
perfect; in this case, the degree of symmetry 
throws some light on the trustworthiness of the 
calculated anisotropy constants. 

In order to calculate the anisotropy constants 
from the torque data, it is necessary to use an 
applied field strong enough to make the torque 
independent of it. It has generally been thought 
that a field of several thousand oersteds is 
sufficient for a disk with not too large a ratio 


1R. M. Bozorth, J. App. Phys. 8, 575 (1937). 

?R. M. Bozorth, Phys. Rev. 50, 1076 (1936). 

3H. Schlechtweg, Ann. d. Physik [5] 27, 573 (1936). 
‘L. P. Tarasov and F. Bitter, Phys. Rev. 52, 353 (1937). 


of thickness to diameter and some evidence to 
this effect was found? when a moderately accu- 
rate torquemeter was used. On the other hand, 
some later work by Schlechtweg? indicated that 
the torque did not reach a saturation value in 
fields of the same order of magnitude. Unfortu- 
nately, his highest field strength was only 2400 
oersteds, and under the conditions of the experi- 
ment this was too low to prove conclusively that 
the torque would not have reached a limiting 
value at a somewhat higher field. By using a 
torquemeter of very high sensitivity, the author 
has been able to show in the present paper that 
a saturation value of the torque cannot be 
obtained experimentally, but can be calculated 
by extrapolating to 1/H=0. 


APPARATUS 


The torquemeter used in this work was of the 
familiar torsion type, but built on a larger scale 
than usual. The electromagnet was about the 
same size as the one illustrated elsewhere,‘ 
mounted so as to rotate around a vertical axis. 
With the flat faces of the pole pieces 1¢ inches 
apart, a field of 4000 oersteds could be obtained, 
but measurements ordinarily could not be made 
above 3500 oersteds because of overheating. The 
field was calibrated in terms of the magnetizing 
current in the usual manner with the aid of a 
search coil and a ballistic galvanometer. The 
fields could be measured with an accuracy of one- 
tenth of a percent relative to each other, but the 
uncertainty involved in reading the ammeter 
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during the actual use of the torquemeter brought 
the possible error in H to about one-half percent. 
The field was uniform to at least this extent over 
the area of the single crystal disks used in this 
investigation, which were mostly half an inch in 
diameter; however, the field probably varied 
somewhat more over the area of the one-inch 
polycrystalline disks which were also used. The 
uncertainty in the field was in all cases too small 
to affect the results appreciably. The highest fields 
used were too low to require any correction for 
the image effect described by Weiss and Forrer.® 

The specimen was held by the pressure of a 
thumbscrew in an aluminum holder which was 
suspended between the pole pieces. Phosphor 
bronze wires 0.032 inch in diameter were used for 
the suspension. The upper wire, which was 
fastened securely to the frame of the torque- 
meter, was 10 inches long; the lower wire, held 
at the other end in a wire chuck attached to an 
accurately machined gear, was 24 inches long. 
This gear was meshed with a worm gear, a 
complete revolution of which corresponded to a 
12-degree twist of the lower end of the suspension. 
On the same shaft with the worm gear was a 
large dial graduated into 100 divisions. The 
mechanical torque required to balance the mag- 
netic torque exerted by the field on the specimen 
could be reproducibly measured to a tenth of a 
division, and since the torque peaks of most of 
the specimens examined amounted to one or 
two complete revolutions of the dial, an accuracy 
of one part in a thousand or better could be 
reached. A galvanometer mirror was mounted on 
the upper suspension wire just above the alumi- 
num disk-holder, and the image of a cross-hair 
thrown on a scale clamped to the torquemeter 
frame showed accurately when the specimen was 
in its zero position, i.e., the position in which the 
magnetic and mechanical torques were properly 
balanced. To convert the arbitrary dial readings 
to absolute units it was found best to calibrate 
the suspension by means of a specimen whose 
torque peaks had been measured on a balance- 
type torquemeter, which gave the results in 
absolute units directly.* It was thus found that 
— and R. Forrer, Ann. de physique [10] 12, 297 
; ‘This was done through the kindness of Dr. A. R. 


Kaufmann at the Massachusetts Institute of Technology, 
who used the torquemeter described in reference 4. 





one dial division corresponded to 65.4 dyne-cm. 
In making the measurements, it was necessary to 
introduce a slight correction for the torque 
exerted by the aluminum disk-holder which 
probably arose from ferromagnetic impurities. 
Because the precision was so high it was also 
necessary to correct for any temperature change 
during a set of readings. This correction was 
found experimentally to be one part in 700 for 
1°C, a value very close to the known tempera- 
ture variation of the anisotropy constant of iron 
at 25°C. The temperature was read on a ther- 
mometer suspended from the frame of the torque- 
meter in such a way that the bulb was only a 
millimeter or two above the specimen and 
slightly to one side of it. 


MEASUREMENTS ON SINGLE CRYSTAL DISKS 


The first set of measurements was made on.a 
group of fourteen well-annealed single crystal 
disks of iron containing from 2.4 to 7.4 percent 
of silicon by weight,’ as found from subsequent 
chemical analyses on the disks. Eleven of these 
had the (100) plane in the plane of the disk, 
the others had (110). Back-reflection Laue pat- 
terns of the disks proved that the deviations from 
these ideal orientations were a few degrees at 
most, not enough to affect the validity of the 
results to be described. After the measurements 
on these disks were complete, some of them were 


TORQUE IN SCALE OrvISiONS 





OM, 


Fic. 1. Variation with 1/H, of the magnetic torque 
peaks of a (100) disk. The successive peaks are numbered 
in order. H, is the applied field. 


7 These disks, to be described in greater detail in a subse- 
quent paper on the anisotropy of tron-silicon alloys, were 
prepared several years ago by Dr. K. J. Sixtus, at that 
time of the General Electric Research Laboratory. 








1226 L. P. 


polished and etched deeply to reduce the thick- 
ness and were remeasured. There were thus 
available data on disks ranging from 0.37 to 
0.75 inch in diameter, and from 0.008 to 0.035 
inch in thickness. The ratio ¢/D of thickness to 
diameter varied from 0.015 to 0.090. 

For the present the discussion will be restricted 
to the (100) disks, whose four torque peaks are 
of equal magnitude. In Fig. 1 the torque peaks 
found for one of these disks are plotted indi- 
vidually as a function of 1//7,, where H, is the 
applied field. This was varied between 2000 and 
3500 oersteds. It is clear that at least in this 
range the magnitude of the torque peak changes 
linearly with 1/77, to a high degree of accuracy. 
The results can be expressed by the equation 


Tu=T,(1—Ho ‘i ), (1) 


where 7y and 7, are the magnitudes of the 
torque peaks for an applied field H, and for an 
infinite one, respectively, and Ho is a constant. 
This linear relationship was found to hold only 
above 2000 oersteds or thereabouts, the curve 
breaking suddenly downward for lower fields. 
Above 2000 oersteds the law held all the way to 
1/77,=0, as was demonstrated by measurements 
on an oblate ellipsoid of revolution to be de- 
scribed later. For disks or ellipsoids cut parallel 
to the (100) plane, the anisotropy constant is 
given by K,=27,.. 

At this point, it is desirable to mention the 
effect of rotational hysteresis upon the results. 
It is known that the two torque curves obtained 
by rotating the field in opposite directions are 
separated vertically by an amount proportional 
to the rotational hysteresis. The data plotted in 
Fig. 1 were obtained by always rotating the field 
in the same direction through all four peaks. 
Hence, two of the peaks were too small by an 
amount equal to the torque arising from rota- 
tional hysteresis, while the other two were large 
by the same amount. Since rotational hysteresis 
tends toward zero as magnetic saturation is 
approached, the error becomes less in higher 
fields with the result that the slope of two of 
the lines, such as are drawn in Fig. 1, may 
differ slightly from the slope of the other two 
lines. In some cases, as in Fig. 1, the error arising 
from rotational hysteresis was so small that all 
four lines were parallel within experimental 
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Fic. 2. Ho as function of t/D, the ratio of the thick- 
ness of a disk to its diameter. The straight line on the right 
was calculated by least squares using only the circles. The 
three straight lines at the left, drawn in directly, connect 
the three sets of data for cold-rolled and strain-relieved 
materials. 


error ; in other cases, especially for polycrystalline 
disks, the error was large enough to give a 
perceptible difference in slopes. Were it not for 
imperfections in the specimen, such as a slight 
deviation between the (100) plane and the plane 
of the disk, the four lines of Fig. 1 would be 
coincident. Even where rotational hysteresis was 
comparatively large, it was far too small to 
account for the observed differences in the sizes 
of the supposedly equal torque peaks. To elimi- 
nate the effect of rotational hysteresis, the results 
were always plotted to give a single line showing 
the variation of the sum of the four torque 
peaks with 1/H,, for in this case the errors 
automatically canceled. 

It was noticed that Ho increased with the 
ratio of the thickness of a disk to its diameter. 
Several ways of plotting results were tried, but 
the only one that gave a reasonably linear 
variation was the one in which Hp was plotted 
against log (t/D), where ¢t and D are the thickness 
and diameter. The results for all the (100) disks 
are indicated by circles in Fig. 2. A dashed 
circle shows the result calculated from similar 
data obtained by Schlechtweg’ on a (100) disk. 
This point is in reasonably good agreement with 
the others. 

It is obvious that there is considerable scatter 
in Ho, when only the points for the (100) disks 
are considered. An examination of the original 
data showed that this scatter could not be 
accounted for on the basis of variations in the 
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silicon content of the disks which, as has already 
been mentioned, contained from 2.4 to 7.4 
percent, nor because of slight deviations from 
the ideal orientation. The only explanation of 
the scatter is that it arises from certain peculiari- 
ties of the specimens whose nature is unknown 
at present, for the amount of scatter is larger 
than can be ascribed to errors in the measure- 
ments. The straight line drawn through these 
points was calculated by least squares, and its 
equation is 


HT) = 637 +246 logio (t/D), (2) 
or in a slightly different form, 
FH )= 246 logio (385¢/D). (2a) 


It is obvious that this equation is valid only 
between certain limits of ¢/D; moreover, it 
applies only to the torque peaks of (100) single 
crystal disks. The line corresponding to this 
equation in Fig. 2 is drawn somewhat beyond 
the experimental points to indicate the probable 
range of validity. 

The discrepancy between the experimental 
values of H» and those corresponding to the line 
was generally less than 5 percent. If it is assumed 
that H» can be calculated for a (100) single- 
crystal disk from its dimensions with an error 
not exceeding 5 percent, then the torque peaks 
can be extrapolated to 1/H,=0 according to 
Eq. (1) with an error of less than 0.5 percent, 
using the reasonable values of H,=3000 and 
H,=300. Thus this equation should be useful 
when experimental circumstances make impos- 
sible the accurate determination of Ho. 

We now turn to the (110) disks, which have 
two major and two minor peaks instead of four 
of the same size. Similar measurements on the 
(110) disks showed that H for the minor peaks 
is considerably larger than for the major ones. 
The two sets of points are plotted in Fig. 1 as 
crosses. The scatter is considerable because few 
disks were available and also because they were 
smaller than the (100) disks so that the variation 
of the peaks with 1/H, was harder to determine 
accurately. There can be no doubt, however, 
that Hp» for the (100) disks is smaller than for the 
minor peaks of the (110) disks and larger than 
for their major peaks. At the same time, it 
appears that HH, for both sets of (110) peaks does 





not increase as rapidly with ¢/D as does Hy for 
the (100) peaks. More data on (110) disks are 
needed before it will be possible to calculate Hy 
from the dimensions with any reasonable ac- 
curacy. 


MEASUREMENTS ON POLYCRYSTALLINE DISKS 


Some 3 percent silicon steel was cold rolled 
from 100 mils down to various thicknesses 
ranging from 20 to 3 mils. Disks one inch in 
diameter were punched from this material, which 
is described in detail elsewhere,* and were 
measured both in the cold-rolled condition and 
after being annealed just below the recrystalliza- 
tion temperature to relieve the strains as much 
as possible without changing the texture. The 
predominant texture of severely cold-rolled 
silicon steel is one having the (100) planes in the 
rolling plane or deviating by various amounts 
from this ideal orientation. The torque curve of 
this material is similar to that of a (100) single 
crystal disk in that the four torque peaks are of 
approximately the same size; furthermore, the 
torque curve remains practically the same after 
the disks have been annealed for strain relief. 
Some work was also done on three disks of fairly 
pure iron, all of which were cold rolled to a 95 
percent reduction in thickness. 

The three sets of results are shown, dis- 
tinctively marked in Fig. 2 on the left. Although 
the points are very scattered, it is possible to 
draw three straight lines, one through each set. 
The same disks were used before and after the 
strain relief so that the decrease in Ho following 
the anneal is known to be a real effect. The values 
of Hy for the three iron disks are even smaller, 
but not enough specimens were available to make 
sure that this was not an accidental variation. 
Although the lines drawn through the three sets 
of points have roughly the same slope as the 
line for the (100) single-crystal disks, the lines 
themselves are all somewhat higher. One reason 
for this is the difference between the strain 
conditions in the two cases. It may also be that 
instead of drawing a straight line for the single- 
crystal disks and another one for the strain- 
relieved disks, there should be drawn a slightly 
curved line passing through both sets of points. 


8’ L. P. Tarasov, Trans. A. I. M. E, 135, 353 (1939). See 
Table I. 
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Fic. 3. Sum of four torque peaks of the (100) ellipsoid 
as function of 1/H,. This ellipsoid was filed from the (100) 
disk referred to in Fig. 1. The straight line is drawn to have 
a width of 0.1 division, the limit of accuracy in the meas- 
urements, 


The considerable scatter that occurs in Hp» for 
polycrystalline material shows that it is best not 
to depend too much upon the values shown in 
Fig. 2. 


MEASUREMENTS ON SINGLE CRYSTAL ELLIPSOIDS 


After the measurements were finished on the 
(100) single-crystal disk for which the results are 
shown in Fig. 1, this disk was filed into a highly 
accurate oblate ellipsoid of revolution.* The 
torque peaks of this ellipsoid, which contained 
2.4 percent silicon, were measured with the 
highest possible precision after it had been etched 
and annealed to remove the distorted surface. 

The results, shown in Fig. 3, indicate that an 
ellipsoid follows the same type of law as do the 
disks. However, Ho is only 21.1 as compared with 
328 for the same single crystal when it was still 
a disk. The accuracy with which the 1/H, law 
is followed can be seen from the thickness of the 
line which corresponds to 0.1 division, the limit 
of accuracy of a single reading. The linear 
relationship extends to considerably lower fields 
than for disks, entirely on account of the relative 

*The method is essentially the one described by L. W. 
McKeehan, Rev. Sci. Inst. 5, 265 (1934). This and another 


ellipsoid were prepared at M. I. T. through the kind cooper- 
ation of Dr. A. R. Kaufmann. 
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ease of “saturating” an ellipsoid and a disk. 
The deviations from linearity are so negligible 
that the relatively short extrapolation to1/H,=0 
needs no, further justification. 

The work on the ellipsoid made it possible to 
check the validity of the linear extrapolation to 
1/H,=0 in the case of the (100) disks. The 
effective anisotropy constants were calculated 
(K,=2Ty) for various fields for both the disk 
and the ellipsoid into which it was made, and 
the values were plotted as in Fig. 4. The 
anisotropy constant extrapolated to 1/H,=0 
was the same within experimental error in the 
two cases. This establishes the validity of Eq. (1) 
for the (100) single crystal disks, and it is 
reasonable to assume that the same is true of 
polycrystalline disks. 

A (110) disk containing 5.9 percent silicon was 
similarly filed into an ellipsoid and it was found 
that Hp for the major peaks was very close to 20. 
The variation of the minor peaks with 1/H,, 
which was hard to determine with any certainty, 
gave Hy equal to about 15. The validity of the 
linear extrapolation to 1/H,=0 for (110) disks 
was also checked and Eq. (1) was again found to 
be satisfactory in the extrapolated range. 


TORQUE BEHAVIOR IN LOWER FIELDS 


It has already been mentioned that the linear 
relation given by Eq. (1) was followed only for 
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Fic. 4. Apparent anisotropy constants calculated for 
various fields for a (100) single crystal disk and for the 
ellipsoid into which it was filed. The true anisotropy con- 
stant is the value of K, for 1/H,=0. 
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fields above a certain minimum value, which in 
the case of disks generally was in the neighbor- 
hood of 2000 oersteds. Disks with a high ratio 
of thickness to diameter required a stronger 
field for the 1/H, law to hold than did disks 
with a low ratio. This behavior is clearly con- 
nected with the relative ease of approaching 
magnetic saturation in the disks in the two 
cases. 

For the (100) ellipsoid, the minimum satis- 
factory field was considerably lower than for a 
disk of the same dimension ratio, as would be 
expected on the basis of ease of saturation. 
The deviations from the 1//7, law were negligible 
down to 1100 oersteds, but just below this there 
was an abrupt drop in the magnitude of the 
torque peaks. Had the ellipsoid been thinner 
with respect to its diameter, the break would 
not have occurred until even lower fields were 
reached. 

The major peaks of the (110) ellipsoid be- 
haved exactly like the peaks of the (100) ellipsoid, 
but it was a different story for the minor peaks. 
The variation of a major and a minor peak of 
the (110) ellipsoid is shown in Fig. 5, where it 
can be seen that at 1000 oersteds or so the minor 
peak is 25 percent larger than the value it 
attains at much higher fields. Below 800 oersteds, 
the major and minor peaks are of practically 
the same magnitude. A satisfactory explanation 
of this phenomenon has not yet been given. 


DISCUSSION OF THE 1/H, LAW 


The only theoretical work on the approach 
of torque measurements to their saturation value 
is that of Schlechtweg,'® who predicted that the 
torque peaks of a homogeneously magnetized 
disk should approach saturation according to a 
1/H, law. The restriction to homogeneously 
magnetized disks is nothing more than a re- 
striction to the use of ellipsoids instead of disks. 
Calculation of Ho for the torque peaks of (100) 
and (110) ellipsoids according to Schlechtweg’s 
equations gave H)=0 in both cases, although for 
other orientations it was not zero. This indicated 
that the experimentally observed values of Ho 
for the two ellipsoids had to be explained by a 


10 Reference 3; also H. Mussmann and H. Schlechtweg, 
Ann. d. Physik [5] 32, 290 (1938). 





different theory than the one proposed by 
Schlechtweg. 

It is interesting to compare the approach to 
saturation of the torque peaks with the similar 
approach to saturation found for the intensity 
of magnetization in the case of iron and iron- 
silicon alloys.®: “ One highly important difference 
between the two cases is that the applied field 
is used in connection with the torque measure- 
ments while the effective field is used for the mag- 































































































ae enenreenseresiny meee eee -—— 
a — a 
120 MAJOR PEAK 
| 
= a a 
| | 
| 
100 fl T = es 
| 
| 
o 
FS | 
5 80 ; -- = 
2 | —-— EXTRAPOLATED VALUE 
~ OF TORQUE PEAK 
a 
3 | ACCORDING TO 7b LAW 
| Q 
Ses _ 
w 
> } 
3 AN | 
« | 
°o T ee — 
. ry 
een ee ee — “> “> 
| Ll R 
40 , | a | MINOR PEAK 
| 
| | | 
| | 
| a ee nn 
= = i ad 
] 
2000 3000 





HAIN OERSTEOS 


Fic. 5. Major and minor torque peaks of (110) ellipsoid 
plotted against magnetic field. 


netic intensity. The coefficient @ appearing in 
place of Hy in the similar equation of approach 
to saturation of magnetization has been measured 
by Fallot," who found that a varied between 5 
and 24 for iron-silicon alloys. The values ob- 
tained depended entirely upon some unknown 
property of the specimen measured and not upon 
the silicon content. Hence, all that can be said 
about Hy and a is that they are of the same order 
of magnitude. 

Before trying to reconcile this apparent simi- 
larity of the coefficients in spite of the difference 


11 A. Kussmann and E. Schoen, Physik. Zeits. 38, 777 


(1937). 
2 M. Fallot, Ann. de physique [11] 6, 305 (1936). 
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between the meanings of //7, an attempt was made 
to replot the torque data for one of the (100) 
disks on the basis of the effective field rather 
than the applied one. In order to get a linear 
approach to saturation using the effective field 77, 
it was necessary to use a law.of the form 


Ty=T,(1—Hy'/H"), (3) 


where 7 was very nearly equal to §. This equation 
is useless because neither m nor 7, can be ob- 
tained with any certainty from the data. Thus 
the only simple and easily obtainable equation 
that fits the torque data is the one involving 
the applied field. An examination of Fallot’s data 
shows that it is equally hopeless to replot them 
on the basis of the applied field rather than the 
effective one and still retain a linear law; such 
a change, moreover, could not be _ justified 
theoretically. 

The best explanation of the apparent similarity 
in the two laws of approach seems to be the 
following: The 1/H law for the approach to 
magnetic saturation holds only for polycrystal- 
line material, where the direction of magnetiza- 
tion is averaged over all orientations, because it 
is known that the magnetization curves of a 
single crystal along the directions of easy mag- 
netization do not conform to this law; instead, 
saturation is reached at very low fields for the 
[100] direction in iron, and in somewhat higher 
fields for the [110] direction. With the torque 
peaks, on the other hand, the direction of mag- 
netization is along the same crystallographic 
direction within a few degrees, so that the 
magnetization along this direction may well 
approach saturation according to some other 
law than 1/H. Since the magnitude of the 
torque peak is, among other things, proportional 
to the intensity of magnetization, it follows that 
the approach to saturation in torque will also 
follow this law. It seems, then, that the difference 
in crystallographic conditions between the two 
cases is accidentally balanced by the difference 
in the nature of the magnetic fields, i.e., applied 
or effective, leaving the laws formally the same 
and with coefficients of the same order of 
magnitude. 
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As has been mentioned, Schlechtweg’s torque 
equations in 1/H, predict a zero value for H for 
the torque peaks of the (100) and (110) ellipsoids, 
and nonzero values for other orientations. For 
the two special orientations, /7) was experi- 
mentally found to be in the neighborhood of 20, 
which appears to be an empirical constant of 
the material. It is probable that the torque 
peaks of an ellipsoid of some other orientation 
than (100) or (110) would have values of Ho 
somewhat higher than this, the excess being 
given by Schlechtweg’s theory. 

Thus far, the discussion has been concerned 
only with ellipsoids. For disks, it is clear that 
the very much larger values of 7» are associated 
with the difficulty of saturating their edges, i.e., 
the material outside the inscribed ellipsoid. This 
is in line with the results shown in Fig. 2. 
Probably a very small part of Ho for a disk-is 
to be attributed to the approach to magnetic 
saturation of the specimen as a whole, as in 
the case of the ellipsoids, but the important part 
is due to edge effects. 


CONCLUSIONS 


There are really three distinct 1/H, laws for 
torque peaks. The first is Schlechtweg’s theo- 
retical one, which predicts an increase with 
field in the peaks of single crystal ellipsoids 
unless they are of the (100) or (110) orientation. 
Next, there is the experimentally determined 
1/H, law for (100) and (110) ellipsoids, which is 
attributed to the slow increase in magnetic 
saturation noticeable even in very high fields. 
And finally, there is the 1/H, law observed for 
disks, for which Hy is quite high because it is 
connected with the approach to saturation of the 
edges. 

The usefulness of the 1/H, law for disks lies 
in the fact that it makes possible the calculation 
of a unique value of the anisotropy constant, the 
one corresponding to 1/H,=0, which does not 
depend upon the field used nor upon the dimen- 


‘sions of the specimen. The validity of the 


extrapolation to 1/H,=0 has been satisfactorily 
checked by the measurements on two ellipsoids. 
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The ferromagnetic anisotropy constant K, has been 
accurately measured by the magnetic torque method for 
a group of iron alloys containing up to 13.7 atomic percent 
silicon. The points fall on two intersecting straight lines 
given by the equations: 

K,X10-5=5.29—0.279A A <9.86 
and Ki, X10-°=4.43-—0.1915A A>9.86 


K;, is expressed in ergs/cm' and A is the atomic percentage 
of silicon. The value for iron is considerably higher than 
the one due to Akulov which has often been quoted in the 
literature. It is shown that the latter is unreliable because 
it is an average value based on the highly discordant 
results of several investigators. Torque measurements 
were made on the identical ellipsoid of iron for which 


Piety measured the magnetization curves along the simple 
crystallographic directions and a large discrepancy was 
found between the values of K, calculated by these two 
methods. This leads to the belief that the present-day 
theory of ferromagnetic anisotropy is in need of some 
revision, possibly in connection with the higher order 
terms in the magnetic energy equation. When K, is plotted 
against the atomic percentage of silicon, the change in 
slope occurs at very nearly the same concentration as the 
similar change in slope found by Jette and Greiner for the 
lattice parameter do. An attempt is made to explain the 
occurrence of these breaks in terms of a change from a 
hypothetical superlattice with a silicon concentration of 
3/32 to a partially formed Fe;Si superlattice. 





INTRODUCTION 


KNOWLEDGE of the ferromagnetic ani- 

sotropy of iron-silicon alloys is useful from 
both the practical and the theoretical points of 
view. The study of the preferred orientations of 
the grains of cold-rolled strip has its greatest 
importance in connection with transformer steels, 
which contain up to 6 percent by weight of 
silicon. Any thorough investigation of this 
texture by the magnetic torque method requires 
some knowledge of the anisotropy constant, to 
which the measured torque is directly propor- 
tional.! 

The variation of the anisotropy constant with 
silicon content is also of interest in connection 
with the known anomalous behavior of the alloys 
in the a@ solid solution range.? Several investi- 
gators have measured the anisotropy of iron, the 
most recent work being done on material of high 
purity, and others have measured the anisotropy 
of individual iron-silicon alloys. The results on 
iron have been in poor agreement with each 
other. As for the iron-silicon alloys, all that 
could be said was that the anisotropy constant 
decreased rapidly with the addition of silicon. 


1L. P. Tarasov, J. App. Phys. 9, 192 (1938); and Trans. 
A. I. M. E. 135, 353 (1939). 

2E. R. Jette and E. S. Greiner, Trans. A. I. M. E. 105, 
rg and M. Fallot, Ann. de physique [11] 6, 305 


Four methods of determining the anisotropy 
are available, three of which are described in a 
single article by Williams.* These are magnetic 
torque studies, the measurement of the areas 
between magnetization curves along various 
simple crystallographic directions, and the com- 
parison of the shapes of the magnetization curves 
with the theoretically predicted ones. The fourth 
method involves the direct measurement of the 
normal component J, of the magnetization by 
means of a suitable search coil, as was done by 
Honda and Kaya.‘ The derivation of the anisot- 
ropy constant from such measurements is 
described by Bitter.® 

For simplicity, it has been tacitly assumed in 
the foregoing that there is only a single anisot- 
ropy constant. Actually, the experimental data 
may require more than one such constant, as has 
been shown by Bozorth.* For cubic lattices, the 
first two anisotropy constants, K, and Ko, are 
defined by the equation 


E=K ot Ki(ara?+a°as+a3*a;") 
+Koea;*ar’a3*, (1) 


where E is the energy of magnetization in a given 





3H. J. Williams, Phys. Rev. 52, 747 (1937). 

*K. Honda and S. Kaya, Sci. Rep. T6hoku Imp. Univ. 
15, 755 (1926). 

5F. Bitter, Introduction to Ferromagnetism (McGraw- 
Hill, 1937) p. 201. 

®*R. M. Bozorth, Phys. Rev. 50, 1076 (1936). 
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direction, and the a’s are the cosines of the angles 
between this direction and the cubic axes; along 
these axes, Ky is the energy of magnetization. In 
the case of iron and iron-silicon alloys, it is 
fortunate that at room temperature the effect 
of Ke upon the results is quite small compared 
with that of K;. Were this not so, it would be 
very difficult to use torque data in studying the 
texture of polycrystalline material. When torque 
measurements are used to determine the form 
of the anisotropy, it is found from considerations 
of symmetry that for a disk having (100) lying 
in its plane only K, appears in the equation of 
the torque curve; but if (110) lies in the plane 
of the disk, K2 also has to be taken into account. 
Since the measurements reported below were 
mostly on (100) disks, information was obtained 
in those cases only about K,, which will generally 
be referred to as the anisotropy constant without 
further qualification. 

The rest of this paper can be conveniently 
divided into three parts. The first of these 
describes the experimental work on iron and 
iron-silicon alloys done by the writer. The next 
part is a critical discussion of the discrepancies 
between these results and the ones: that have 
been previously reported, especially for unalloyed 
iron. In the last part, the anisotropy data ob- 
tained for the iron-silicon alloys are compared 
with other types of data in an attempt to find 
some reasonable interpretation of certain discon- 
tinuities that have been observed in the a@ solid 
solution range. 


Part I 


Specimens 


The specimens were single crystal disks of iron 
containing various percentages of silicon and cut 
so as to have either a (100) or a (110) plane 
parallel to the plane of the disk. The single 
crystals were prepared’ by melting the required 
proportions of electrolytic iron and the 90 percent 
grade of ferrosilicon in small beryllia crucibles in 
an atmosphere of hydrogen. The ingots were 
made to solidify as large crystals by lowering the 
crucibles at a suitable rate through a vertical 
tube furnace. Slabs of the desired orientation, as 


7The work up to and including the grinding of the 
properly oriented disks was done by Dr. K. J. Sixtus. 
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determined by x-rays, were cut from the large 
crystals and from these slabs disks were ground 
about } inch in diameter and from 0.020 to 0.035 
inches in thickness. The disks were polished and 
slightly etched to make certain of their mono- 
crystalline nature. 

After a set of torque measurements had been 
made on the disks in this condition, they were 
etched more deeply to remove any surface 
material that was deformed during the grinding 
and polishing operations; following this, they 
were annealed for 2 hours at 1100°C in pure dry 
hydrogen and were then furnace cooled. A surface 
etch showed that the disks were still single 
crystals and did not contain any small included 
grains. The anisotropy data on the iron-silicon 
alloys were obtained, with one exception, from 
the torque measurements made on these annealed 
disks. The exception was a disk that was filed 
into an ellipsoid of revolution’ and then annealed 
like the others. It may be mentioned here that 
the effect of annealing was usually, but not 
always, to increase K, by 1 or 2 percent. 

When all the torque measurements had been 
completed, the disks were cut into several pieces 
which were analyzed individually for silicon with 
a probable error of 0.02 weight percent. This was 
necessary because in some instances the silicon 
had not diffused uniformly throughout an ingot, 
so that considerable variations in the silicon 
content existed within a single disk. Only those 
disks were used in the final calculation of the 
variation of K, with silicon content for, which 
the analyses of the separate portions ‘agreed 
within 0.1 percent of silicon. These individual 
analyses were suitably averaged to give the 
effective composition of the specimens. 

The carbon content of the original ingots 
varied from 0.02 to 0.04 percent, a range that 
can be accepted as substantially correct for the 
disks themselves since the only two disks that 
were anlyzed for carbon gave results in good 
agreement with those for the corresponding 
ingots. The 2-hour anneal at 1100°C probably 


8 This was done at the Massachusetts Institute of Tech- 
nology through the kindness of Dr. A. R. Kaufmann; the 
same holds for the measurement of the torque peaks of this 
ellipsoid for the purpose of calibrating the writer’s torque- 
meter suspension in absolute units. 

®Unless otherwise stated, weight percent is meant 
throughout the paper. 
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Fic. 1. Anisotropy constant K, plotted as a function of 
atomic percentage of silicon. 


did not decrease the carbon that was present 
because the hydrogen was dry whereas moist 
hydrogen is required for the reduction of the 
last traces of carbon. In addition to carbon, 
other impurities must have been present since 
the 90 percent ferrosilicon from which the alloys 
were made is known to contain 1 to 3 percent of 
total impurities other than iron, mostly alu- 
minum together with some calcium and a little 
titanium.!® Unfortunately, these impurities were 
not analyzed for and all that can be said is that 
the disks contained these additional impurities 
to the extent of 1 to 3 percent of the silicon 
content. 

Besides these iron-silicon disks, two ellipsoids 
of high purity iron were measured on the mag- 
netic torquemeter. They were the identical 
ellipsoids used by Piety" in 1936 in studying the 
anisotropy constants of iron at several tem- 
peratures, and were measured in the condition 
in which they were left at the conclusion of his 
work. 


1 According to Table 4, p. 19, Alloys of Iron and Silicon 
by E. S. Greiner, J. S. Marsh and Bradley Stoughton 
(McGraw-Hill, 1932). 

1 R. G. Piety, Phys. Rev. 50, 1173 (1936). One of the 
ellipsoids measured by the writer was the one identified as 
P10, having (110) in the equatorial plane; the other one 
was the (100) ellipsoid, which Dr. Piety did not test suff- 
ciently to enable him to calculate K,. The second (110) 
ellipsoid used by Dr. Piety, P12, had been destroyed prior 
to the present investigation. 





Magnetic torque measurements 


The magnetic torque curves from which the 
anisotropy constants were calculated were ob- 
tained on the torsion-type torquemeter described 
elsewhere by the writer." The measurements 
consisted of the four peaks of the torque curve, 
each taken at six different values of the applied 
magnetic field H,, which were suitably spaced 
between 2000 and 3500 oersteds. The peaks were 
plotted against 1/H, and extrapolated to 
1/H,=0, for it has been established” that the 
correct values of the anisotropy constants are 
obtained in this manner, regardless of whether 
the specimen investigated is an ellipsoid or a 
disk. The calibration of the torquemeter sus- 
pension in absolute units and the necessary 
corrections to the measured values of the peaks 
are described in this same reference; and the 
curves there shown for the variation of torque 
with 1/H, are representative of the ones from 
which the anisotropy data were derived. When 
the iron ellipsoids were measured, it was neces- 
sary to replace the usual 32-mil phosphor bronze 
suspension wire with one 20 mils in diameter, in 
order to increase the sensitivity sufficiently to 
compensate for the extremely small size of the 
specimens. 

The torque peaks were accurately reproducible 
and the values obtained by a linear extrapolation 
to 1/H,=0 are unlikely to be in error by more 
than 0.5 percent. When K;, is plotted as a func- 
tion of the silicon content as in Fig. 1, the possible 
errors in the two coordinates are of the same size. 


Calculation of anisotropy constants 


It should be mentioned here that the speci- 
mens were not all perfectly oriented, i.e., the 
(100) or (110) plane did not lie exactly in the 
plane of the disk. Since the orientation was 
known to within a degree from x-ray pictures 
obtained by the back-reflection Laue method, it 
was possible in each case to calculate the cor- 
responding change in the magnitude of the 
torque peaks by using the equations for the 
coefficients of the torque equation as given by 
Tarasov and Bitter." 

#21. P. Tarasov, Phys. Rev., this issue, preceding paper. 
(‘Dependence of Ferromagnetic Anisotropy on the Field 


Strength.”’) 
183. P. Tarasov and F. Bitter, Phys. Rev. 52, 353 (1937). 
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The torque equation for disks which are nearly 
of the (100) type reads in the notation of this 
reference, 


T=(K,/2)(A; sin 20+Az sin 46), (2) 


the term in Ke being negligibly small, as are the 
coefficients B, and Bz of the cosine terms. A, is 
zero for the perfect (100) orientation and quite 
small for orientations close to it. In the latter 
case its only effect is to increase slightly two of 
the four torque peaks and decrease the other 
two, so that the average of the four peaks 
remains exactly equal to }K,A>2. The error intro- 
duced by the uncertainty in Ag is negligible in 
comparison with other errors. 

The situation is considerably more complicated 
for disks of the (110) type because K, affects the 
torque curve most for just this orientation, and 
lengthy computations are necessary if the disks 
deviate from this ideal orientation. It is experi- 
mentally impracticable to extrapolate the whole 
torque curve to 1//7,=0 because the nature of 
the extrapolation has been determined only for 
the peaks. For this reason, instead of finding the 
values of K, and Kz to give the best fit for the 
entire curve, they were chosen to agree with the 
extrapolated values of the peaks. The results in 
the two cases would be the same only if the 
anisotropy can be represented sufficiently well 
by two constants; otherwise, a slight difference 
in K,, and a considerable one in Ke, could arise. 
These constants were calculated from the extra- 
polated values of the torque peaks with the help 
of these equations: 


K,=1.155C,—0.693C, 
and 


Here C; is the sum of the magnitudes of the two 
major peaks and C2 of the two minor peaks. Eq. 
(3) is valid only for K, positive and Ke not too 
large; with certain changes of sign, it would 
hold for K, negative. 

The volumes of the disks, required for express- 
ing the anisotropy in the customary ergs/cm’*, 
could not be determined with sufficient accuracy 
by direct measurement; instead, they were ob- 
tained from the weights and densities of the 
disks, the latter having been calculated from the 
lattice parameters of the iron-silicon alloys. The 
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Fic. 2. Same points as in Fig. 1, replotted against the 
weight percentage. Results of previous investigations are 
included. 


reason for not using measured densities was the 
possibility of fairly large errors arising from 
voids and imperfections likely to occur in cast 
material that had never been subjected to any 
mechanical treatment. The lattice parameters as 
given by Jette and Greimer? fall quite well on 
two intersecting straight lines when plotted 
against the atomic percentages A of silicon. The 
equations of the two lines as calculated by least 
squares are: 


ayo= 2.8610 —0.000723A, A<9.55, 
and 


ay = 2.8744—0.002130A, A>9.55, (4) 


where do is the lattice parameter in angstrom 
units. When 28.06 and 55.84 are taken as the 
gram-atomic weights of Si and Fe, the densities 
are also found to fall on two intersecting straight 
lines, which are given by 

p=7.864—0.0332A, A<9.55, 
and 

p=7.767—0.0232A, A>9.55. (5) 


The volumes of the disks were derived according 
to Eq. (5) from the densities corresponding to the 
proper compositions. 
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Variation of K, with silicon content 


The anisotropy constants of a number of 
iron-silicon single crystal specimens and of two 
high purity iron ellipsoids were calculated in the 
foregoing manner. The values of K, at 25°C are 
plotted in Fig. 1 against the atomic percentage 
of silicon. Of the fourteen disks originally avail- 
able, five were rejected for too great a variation 
in silicon content, for too great a deviation from 
the ideal orientation, or for surface or edge 
defects. Including the two ellipsoids of iron, 
there are eleven experimental points, of which 
six fall on one straight line, three on an inter- 
secting one, while the other two points fall 
below the solid lines. The equations of the two 
lines, as calculated by least squares, are: 


K,X10°=5.29—0.279A, A<9.86, 
and 
K,X10°=4.43—0.1915A, A>9.86. (6) 


The atomic rather than the weight percentage is 
theoretically preferable for the usual physico- 
chemical reasons when dealing with solutions. A 
linear extrapolation, not to be taken too seri- 
ously, shows that K, becomes zero at 23.1 
atomic percent or 13.1 weight percent of silicon. 

For general use, it is more convenient to plot 
K, against weight percent, as is done in Fig. 2. 
It happens that the points still lie on two straight 
lines because of the near-linearity existing 
between the atomic and weight percentages in 
the two ranges of composition separated by the 
intersection of the straight lines. The equations 
of the lines in terms of the weight percentage W 
are: 

K,X10°=5.29-—0.532W, W<5.20, 


and 
K,X10%=4.21—0.324W, W>5.20. (7) 


Returning to Fig. 1, it is important to note 
that the break in slope there shown cannot be 
ascribed to the break which occurs in the density 
curve at about the same composition. If the 
original torque data are used to get the anisot- 
ropy constant in terms of ergs/g so that density 
values are not required, the results still plot as 
two intersecting straight lines and the break is 
of very nearly the same magnitude as that shown 
in Fig. 1. This, of course, arises from the very 


much slower decrease of density than of anisot- 
ropy with increasing silicon content. 

Those points in Fig. 2 which were measured 
in the present investigation are only for the 
etched and annealed specimens since the results 
obtained before such treatment were not as 
consistent. The two points (unmarked circles) 
that do not fall on the curve are shown because 
there was no experimentally justifiable reason for 
rejecting them. The higher of the two points for 
iron is thought to be more likely the correct one 
since it lies exactly on the extrapolated portion 
of the straight line passing through the other 
points. The definitely lower value for the (100) 
ellipsoid of iron must be attributed to some 
unknown factor and not to experimental error, 
and the same is probably true of the other point, 
for which the discrepancy is much less. 

The anisotropy of iron has been investigated 
by others much more thoroughly than that of 
the iron-silicon alloys, and it will be discussed in 
considerable detail in Part II of this paper. This 
is necessary because the values in the literature 
differ greatly among themselves and the generally 
quoted value of 4.3X10° ergs/cm*® for K, does 
not agree at all with 5.29X10° as found in the 
present study. 


Part II 


Previous work on iron 


The first work from which ferromagnetic ani- 
sotropy can be calculated was that by Beck," 
whose experiments were carried out on iron 
containing 1.62 percent silicon. Webster" made 
torque measurements on (100) single-crystal 
disks whose composition approached a little 
closer to that of pure iron. Sufficient amounts of 
C, S, Si, Mn and P were present to classify the 
specimens as low-carbon steel rather than iron. 
Disk A contained 0.220 percent P, an abnormally 
high amount; while disk B, which was annealed 
in hydrogen so that the carbon was reduced to a 
small but unknown percentage, had a typical 
low-carbon steel composition, containing about 
0.4 percent Mn. Shortly afterward, Honda and 


4K. Beck, Vierteljahrsschrift der naturforschenden 
Gesellschaft in Ziirich 63, 116 (1918). The data are dis- 
cussed in a more available publication by L. W. McKeehan, 
Trans. A. I. M. E, 111, 11 (1934). 

4% W.L. Webster, Proc. Roy. Soc. A107, 496 (1925). 
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Fic. 3. Webster’s original data for single crystals of iron 
plotted against 1/H, to show the extent to which his results 
can be reliably used in deriving the value of K. 


co-workers‘: '® prepared specimens of iron of 
somewhat greater purity, the analysis being 0.40 
percent Mn, 0.021 S, 0.021 P and a trace of Si. 
Carbon was practically eliminated by decar- 
burizing four days in moist hydrogen. The ellip- 
soids prepared by Piety" were of very much 
greater purity than any of the preceding, con- 
taining less than 0.1 percent of total impurities. 
The provisional limit of 0.03 percent Si given by 
Piety was used in plotting the points for iron in 
Figs. 1 and 2. 

Let us first consider the data obtained by 
Webster using the torque method. These data 
are listed by Mahajani,'’ who averaged the 
results to get K,=4.66X10° for disk B and 2.98 
10° for disk A. If the same points are plotted 
against 1/H, as in Fig. 3, it is evident that disk 


B obeys the 1/H, law™ to a high degree of pre-' 


cision except for one point which is clearly in 
error. On the other hand, disk A does not follow 
this law, possibly on account of its abnormal 
phosphorus content. The values of K, obtained 
by extrapolating to 1/H,=0 are 5.26X10* for 
B and 4.60 X 10° for A. The extrapolation for the 
latter is, of course, merely suggestive and it is 
possible that the broken line should rise even 
higher. The curves of Fig. 3 prove beyond ques- 
tion that Mahajani’s averaging process leads to 
low values of K,. The good agreement between 
the recalculated value of K, for disk B and the 
value found by the author for the (110) ellipsoid 
of iron is probably fortuitous, in view of the 
differences in composition. Still, this agreement 
is of interest in that it furnishes evidence for a 

16K. Honda, H. Masumoto and S. Kaya, Sci. Rep. 


Té6hoku Imp. Univ. 17, 111 (1928). 
17G.S. Mahajani, Phil. Trans. Roy. Soc. A228, 63 (1928). 
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considerably higher value of K, than has gener- 
ally been assumed to occur. 

A value calculated by Akulov,'® 4.28 10°, 
has been repeatedly cited in the literature. 
Taking the magnetization curve data published 
by Honda and Kaya,‘ Akulov obtained five 
separate values of Ki, varying from 4.84 to 3.34 
and whose average was 4.22. (For convenience, 
a factor of 10° has been omitted.) In a separate 
calculation, some data from the same source 
involving the maximum value of J,, the normal 
component of magnetization, gave values of 4.94 
and 4.48; and when these were combined with 
the previously mentioned calculations by Maha- 
jani'? (which have been shown to be much too 
low), the much quoted value of 4.28 was ob- 
tained. The large discrepancies between indi- 
vidual values, the experimental uncertainties 
inherent in much of the original data, and the 
error introduced by using Mahajani’s average 
values are the factors that make these values 
given by Akulov unacceptable. But the fact that 
some of the individual values are quite high 
compared to the averages again shows that the 
value of K, found in the present investigation 
should not be considered unusually high. 


K, as measured by two methods 


Piety" made some very careful measurements 
of the magnetization curves of two (110) ellip- 
soids and from these calculated K, to be 4.23 and 
4.38X10° ergs/cm* for specimens P10 and P12, 
respectively. (His original results have been 
corrected slightly to bring them to the standard 
temperature of 25°C used in this paper.) The 
previously described torque measurements on 
P10 showed that K, was 5.27 instead of 4.23, a 
difference far greater than could be attributed 
to any experimental errors. P12 was not available 
for a similar comparative study. 

A comparison of the two methods of evaluating 
K, was also made with the (100) iron ellipsoid, 
for which K, was 4.74 according to torque meas- 
urements and 4.66+0.30 when derived from the 
area between the magnetization curves.'® The 


18N.S. Akulov, Zeits. f. Physik 57, 249 (1929). 

19 These magnetization curves were measured by Dr. 
E. M. Grabbe at Yale University in the same manner as 
Dr. Piety had previously measured the (110) ellipsoids. 
The considerable uncertainty in K, was due to the in- 
stability of the magnetic field at high intensities in the 
apparatus used by Dr. Grabbe. 
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large probable error unfortunately makes it 
impossible to say whether or not the two values 
of K, are in good agreement. 

The only other work in which K, has been 
determined with satisfactory accuracy by more 
than one method on a single material is that 
reported by Williams,? who worked with iron 
containing 3.85 percent Si. He obtained 2.72 
and 2.80X10° for K, by two different methods 
of utilizing magnetization curve data. The torque 
data for a (110) disk gave 2.87, but this value is 
raised to 3.05 by a rough extrapolation to 
1/H,=0. Here again K, is appreciably higher 
when obtained from torque measurements in- 
stead of from magnetization curves. 

This suggests that the present-day theory of 
ferromagnetic anisotropy is in need of some 
revision because the magnitude of K, seems to 
depend upon the method by which it has been 
determined. One possibility is that the dis- 
crepancies would disappear if higher order energy 
terms were used in Eq. (1). 

In one respect, it seems quite reasonable that 
K;, should be at least slightly different in the two 
cases. When it is deduced from the area between 
the magnetization curves along the (100) and 
(110) directions, the two curves come together 
for an applied field of perhaps 2000 oersteds, 
taking the case of the iron ellipsoids used by 
Piety. The area from which K;, is determined is 
thus definitely fixed. The torque data taken on 
the same specimen at this field strength must be 
increased by about one percent to take care of 
the necessary extrapolation to 1/H,=0, so that 
a discrepancy of this amount can be expected 
between the values obtained by the two methods. 
This increase in the torque for an ellipsoid, which 
is probably connected with the observed slow 
increase of the saturation intensity of mag- 
netization,” cannot explain the much larger 
discrepancies that have been mentioned. 


K;, for iron-silicon alloys 


The values of K, for four iron-silicon alloys 
that have been recorded in the literature*: ™. ™, 2° 
are included in Fig. 2. Where necessary, the 
values have been extrapolated as well as possible 
to 1/H,=0. Although Beck" made _ torque 
measurements on disks, these data were not 


20H. Schlechtweg, Ann. d. Physik [5] 27, 573 (1936). 





available in suitable form for the calculation of 
K,, and the point shown is that calculated by 
McKeehan from Beck’s magnetization curves. 
The other three points were derived from mag- 
netic torque studies. These points fall from 5 to 
13 percent below the straight line representing 
the results of the present work. In view of the 
discrepancies between the present results and 
the previous ones, it appears that although the 
anisotropy constants depend mostly upon the 
composition, certain unknown factors have an 
appreciable influence. The very good linear 
relationship that was found to hold for the iron- 
silicon alloys should be ascribed to the identical 
history of the specimens involved. 

It may:be well to mention here Zaimovskii’s*! 
indirect determination of K, for a group of iron- 
silicon alloys, since this reference is the only one 
listed in Gmelins Handbuch.” His values are very 
much lower than those shown in Fig. 2, but they 
can be omitted from consideration because his 
indirect method was later shown* to be without 
foundation. 


Significance of K; 


A few scattered values of Ke were available 
from torque measurements on (110) specimens 
and some results have been reported in the 
literature. They are tabulated in Table I, and it 
is clear that no relationship can be set up between 
Kz and the composition. The lack of any regu- 
larity either is a further indication that the form 
of the correction to the fourth power term of the 
anisotropy is more complicated than has been 
assumed, or it shows that K, is very sensitive to 
the condition of the specimen. 


TABLE I. Values of K2 from torque measurements. 








PERCENT 





Si K:X107% COMMENTS 
None 1.4 P12, measured by Piety. 
None —1.7 P10, measured by Piety. 
None —30 P10, torque data. 
3.0 —0.13 Si content in doubt. 
3.85 1.0 As measured by Williams. 
4.6 —0.9 Not exactly a (110) disk. 
5.9 —0.3 Defective edge, Si only approximate. 








21 A. S. Zaimovskii, Katchestvennaia Stal 1, No. 6, p. 35 
(1935). 

22 Gmelins Handbuch (8th ed., 1936), Eisen, Part D 
p. 15. 

24 W. S. Messkin and B. E. Somin, Zeits. f. Physik 98, 
610 (1936). 








1238 L. P. 


Conclusions 


The results show that under carefully con- 
trolled experimental conditions it is possible to 
obtain values of K,; which vary in a regular 
manner with silicon content, but at the same 
time occasional discrepancies are to be expected, 
whose cause is unknown at present. Measure- 
ments of the anisotropy using two different 
methods on the same specimen indicate that the 
theory may require some revision, and until this 
is done, it may be necessary to distinguish 
between the anisotropy constants obtained from 
torque studies and those obtained from mag- 
netization curves. 


Part III 


Discontinuities in properties in the a-phase iron- 
silicon alloys 


Although the phase diagram of the iron-silicon 
system shows that the a@ solid solution phase 
extends at room temperature to about 26 atomic 
percent of silicon, a number of puzzling discon- 
tinuities in physical properties have been ob- 
served in this range which are difficult to inter- 
pret in terms of a single phase. 

The discontinuity in the slope of the lattice 
parameter curve at 9.55 atomic percent of silicon 
and the similar change in slope for K, at 9.86 
atomic percent have already been described. 
Fallot? made a comprehensive study of the 
variation with silicon content of certain mag- 
netic properties and correlated his discontinuities 
with various superlattices which he proposed. He 
found that the curve of magnetic moment against 
atomic percentage showed sudden changes in 
slope at fractional silicon concentrations of 75, 
+ and }; that the same was true of the Curie 
temperature; while the curve of oo/a290, ¢ being 
the saturation intensity corresponding to the 
absolute temperature indicated by the subscript, 
consisted of two distinctly separated parallel 
parts, with the break occurring at 12.5 atomic 
percent. However, the Curie temperature can 
just as well be replotted to show other discon- 
tinuities and hence cannot be used as the basis 
of any argument. Finally, from a practical point 
of view, it has long been known that silicon steels 
containing more than 5 weight percent Si are 
brittle and that discontinuities in other mechan- 





TARASOV 


ical properties are observed at lower percentages,”4 
but so many variables are always present that a 
really satisfactory explanation has never been 
developed. 

As to the relative trustworthiness of the 
various phenomena, the most dependable discon- 
tinuities are those shown by ay and by ao/e290, but 
to make matters difficult, each of these has only 
one discontinuity between 0 and 25 atomic 
percent and these occur at unmistakably different 
concentrations. The break observed in the slope 
of the K, curve is pronounced enough to cor- 
roborate the existence of some peculiarity in the 
iron-silicon system in the neighborhood of 9.5 
atomic percent. 

Since the behavior of do is accurately known 
to the limit of solid solubility at 26 atomic per- 
cent while K, has been studied only half as far, 
and because it is easier to picture changes in do 
than in K, when changes occur in atomic dis- 
tribution, the subsequent discussion will center 
on the lattice parameter. In view of the very 
nearly identical compositions at which the dis- 
continuities occur in do and Ky, it appears 
reasonable to assume that the same explanation 
holds for both. A more detailed interpretation of 
the observed behavior of K, must wait until 
more is known both about the atomic arrange- 
ments in this range and about their possible 
effect upon the anisotropy. 


Interpretations of the discontinuities in a, 


Phragmén* found that the silicon atoms in a 
thoroughly homogenized alloy containing 25 
atomic percent silicon were substituted in a 
regular manner for the iron atoms in the body- 
centered cubic lattice of the latter. The super- 
lattice cell is a cube, 2ao on the side, where do 
always refers to the size of the unit cell of the 
iron lattice. The silicon atoms are arranged 
tetrahedrally at the alternate corners of a small 
cube, as shown in Fig. 4, reproduced from Jette 
and Greiner’s paper.’ The superlattice has been 
detected by x-rays in alloys containing as little 
as 13 atomic percent silicon. Jette and Greiner 
interpret the data as showing that at 12.5 atomic 
percent every superlattice cell has two of the 


*4E. S. Greiner, J. S. Marsh and Bradley Stoughton, 
Alloys of Iron and Silicon (McGraw-Hill, 1932), p. 75. 
% G. Phragmén, J. Iron and Steel Inst. 114, 397 (1926). 
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Fic. 4. Fe;Si superlattice. (Jette and Greiner, reference 2). 


tetrahedron corners occupied by the available 
silicon atoms, and that as the concentration 
increases, the iron atoms at the other corners of 
the Fe;Si tetrahedron are progressively replaced 
by the silicon atoms. 

Below 12.5 atomic percent, they offer two 
related hypotheses in an attempt to account for 
the sudden change of slope in their lattice param- 
eter data. According to hypothesis A, up to 6 
atomic percent the alloy would be a solid solution 
of silicon in iron; from 9 to 25 atomic percent, 
it would be a solid solution of iron in FesSi; and 
in between, it would be a mixture of the two solid 
solutions. According to hypothesis B, there is 
only a single phase in the whole range, the alloys 
being disordered below 9 atomic percent and 
ordered above. The original article must be 
referred to for a complete discussion of the two 
hypotheses, neither of which could be proved or 
disproved. The important thing is that according 
to either hypothesis, the break in the vicinity of 
9 atomic percent is associated with the formation 
of a superlattice. The picture that is offered is 
that at a silicon concentration of 35, i.e., 9.375 
atomic percent, one-half of the eventual super- 
lattice cells have only one corner of the tetra- 
hedron occupied by a silicon atom, and the other 
half have two corners thus occupied. For lower 
concentrations, not enough silicon atoms are 
present to create this partial superlattice. This 
picture is not very convincing because it is 
assumed that two forms of an incomplete super- 
lattice can coexist; moreover, this picture does 
not even attempt to “explain’’ why the super- 





lattice should first appear at such an unorthodox 
concentration. 


An alternative interpretation 


Since a superlattice has to be admitted into 
any discussion of the a solid solution range, a 
somewhat different hypothesis was developed to 
serve not only as a picture, but also as an 
“explanation” of the breaks in a) and K;,. 
Although this hypothesis, like those already 
described, is of an admittedly ad hoc nature, it 
is offered because it seems plausible and can 
perhaps be subjected to experimental check. 

It is first necessary to describe the superlattice 
which is supposed to exist at a silicon concen- 
tration of ;3;. There is no x-ray evidence for such 
an atomic arrangement, but that can be over- 
looked for the present since the other hypotheses 
were built on an equally weak foundation. We 
first consider the usual body-centered arrange- 
ment of iron atoms and take as the unit a cube 
2a on a side. If the iron atom at the center of 
this unit is replaced by a silicon, and the same 
is done with respect to the iron atoms at the 
alternate corners of the unit, the silicon concen- 
tration will be 4; and the silicon atoms will be 
arranged in a large tetrahedron with a silicon 
atom at its center. The closest Si—Si distance is 
4/3a, and the next closest 2a». Eight such units 
are required for a complete unit cell of the super- 
lattice, which thus contains 116 Fe and 12 Si 
atoms and is a cube 4a» on a side. The projection 
of this superlattice cell on a cube plane is shown 
schematically in Fig. 5. 

The observed breaks in the slopes of ao and K, 
are then explained in the following manner: For 
lower concentrations than 35, the silicon atoms 
are assumed to be arranged at random, but as 
their concentration is increased the lattice posi- 
tions taken by the silicon atoms tend to approxi- 
mate more and more the 3 superlattice just 
described. At 3%;, the superlattice is completely 
formed, but at a higher concentration too many 
silicon atoms are present for this superlattice. 
Possibly the interatomic forces between the 
silicon atoms become great enough, on account 
of the excess of such atoms, to bring them closer 
together than 1/3a» allowed by the large tetra- 
hedron. The atoms then go into lattice positions 
consistent with the Fe;Si superlattice shown in 
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Fic. 5. Projection on cube plane of unit cell of super- 
lattice, proposed to account for sudden changes in slope of 
the curves of lattice parameter and of K, occurring near a 
3s atomic concentration of silicon in iron. The unit cell of 
this superlattice, whose edge is four times that of the iron 
lattice, contains 12 atoms of silicon and 116 of iron. z is 
the third coordinate, expressed in units of do, the lattice 
parameter of iron. 


Symbol Interpretation 
3 z=}, 3, } ---; always Fe. 
P 4 z=0, 1, 2 ---; always Fe. 
0 Si if z=0, 4, 8 ---; otherwise Fe. 
1 Si if z is odd; otherwise Fe. 
2 Si if z=2, 6, 10, ---; otherwise Fe. 


Fig. 4. In other words, the large tetrahedron of 
silicon atoms containing a silicon at its center 
collapses into a much smaller and incomplete 
tetrahedron, the shortest Si—Si distance chang- 
ing from /3do to +/2do. 

This picture, although artificial, is quite con- 
sistent with the observed variation in the lattice 
parameter. Up to a concentration of 3%, or pos- 
sibly slightly higher, a) decreases slowly with the 
addition of silicon since the silicon atoms, which 
are somewhat smaller than the iron ones, are far 
apart. Above this concentration, the shortest 
Si—Si distance suddenly drops and the small Si 
atoms become much more effective in reducing do. 
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The objection to this picture, that a super- 
lattice has not been detected by x-rays for alloys 
containing less than 13 atomic percent silicon, 
can now be briefly considered. Since at such low 
concentrations as 3°; the silicon atoms are spaced 
fairly far apart and since the unit cell of the 
superlattice is 4d9 on a side, it is easily conceiv- 
able that only short-distance order exists, such 
as cannot be detected by x-rays. A similar state- 
ment would hold for concentrations between 3; 
and 3, where only the beginnings of a Fe;Si 
superlattice are supposed to exist. Another pos- 
sible objection to the hypothetical 3; superlattice 
is that the arrangement of the silicon atoms is 
not in accord with the accepted ideas of the 
silicon bonding. In answer to this, it should be 
pointed out that silicon in a metal lattice may 
behave entirely differently than when in a dif- 
ferent type of compound. Incidentally, the same 
objection was unsuccessfully raised to the Fe;Si 
tetrahedron when it was proposed. 

As has already been pointed out, Fallot’s 
results are irreconcilable with the ones obtained 
for do and K;. Since nothing is known about any 
relationship between superlattices and the mag- 
netic properties reported by him, Fallot’s results 
cannot at present be considered as contradicting 
the explanation offered for the sudden changes 
in slope in the case of the lattice parameter and 
the anisotropy constant. 
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Fast Neutron Measurements with Recoil Counters 


For certain purposes it is convenient to have counters 
which respond to fast neutrons only. Proportional count- 
ers,' with gases other than BF;, may be used to count 
recoil nuclei. A fast neutron may collide with one of the 
nuclei of the gas in the counter, and will then produce a 
recoil nucleus which will in general ionize heavily along a 
comparatively short range. If the counter be operated at a 
voltage below its threshold then the size of the pulse ar- 
riving on the central wire will be proportional to the 
number of ions formed in the counter by the passage of the 
ionizing particle. By adjusting the amplifier to discriminate 
against pulses below a certain size, or by adjusting the 
voltage applied to the counter and hence the amount of 
additional ionization produced by each ion as it travels 
toward the central wire, it is possible to count only particles 
producing more than a certain number of ions. In this way 
a counter may discriminate against pulses produced by 
beta- or gamma-rays. Such a counter will have a back- 
ground due to the alpha-rays emitted by the natural 
contamination, which can easily be measured and allowed 
for. 

The counting rate, m counts per second, of a fast neutron 
recoil counter will be a measure of the current i(v) of fast 
neutrons per cm? per sec. of velocity v passing through the 
counter, as given by: 


n=VNMp{ ~ ,,i(o)er(v)de, (1) 


where N is the Loschmidt number, p is the gas pressure in 
the counter in atmospheres, M is the number of nuclei per 
molecule, V is the volume of the counter, ¢; is the cross 
section of the nuclei of the gas in the counter for the 
production of recoils of energy E>Emin, where Emin is the 
smallest amount of energy which, if liberated as ionization 
in the counter, will produce a count, and vmin is the ve- 
locity corresponding to Emin. This latter quantity can be 
determined experimentally for a given counter at a given 
voltage by observing the response to rays of different 
known ionizing powers. The above relation neglects wall 
effects and end corrections, which for certain cases can be 
shown to be unimportant. 

In a specific example, Fig. 1 shows the counting rate of 
such a counter, with a hydrogen-methane filling when 
operated in the presence of a Ra-Be neutron source. It will 
be seen that at low voltages (1450 volts) the most heavily 
ionizing of the recoils produces pulses just big enough to 
record. As the voltage is increased the number of additional 


ions produced in the counter is increased, and recoils of 
less and less energy are counted, until at 1770 volts the 
largest of the pulses produced by beta-rays is big enough 
to be counted. At still higher voltages (1850 volts), the 
counting rate rises rapidly as all the beta- and gamma-rays 
are counted. 

It will be noted that this counter has no plateau, but 
that the counting rate rises steadily with the voltage. 
This absence of plateau is due to the continuous distribu- 
tion of energies of the recoiling nuclei, and is to be con- 
trasted with the behavior of the BF; counters! in the pres- 
ence of slow neutrons. The latter exhibit a plateau since the 
pulses produced by the alpha-particles resulting from the 
boron disintegration are all greater than a given size. 

This type of counter has been used in studying the 
neutrons in the cosmic radiation. A preliminary balloon 
flight has been made with such a counter. The data were 
transmitted by short wave radio to the ground station. 
Interpretation of the results in terms of a neutron flux, 
by use of Eq. (1), is difficult, because of our lack of knowl- 
edge of the energy distribution. However, by making a 
simplifying assumption, a rough calculation can be made. 
Let us assume that protons are the only recoiling particles, 
and further, that ¢, is a constant and about 2 10~* cm*. 
(This corresponds? to 4-Mev neutrons; at lower energies ¢ 
is larger and at higher energies smaller.) The observed n 
was 6X 10-3 per sec. (at 1 meter of water), V was 60 cc, 
MNp was 8X10"*, and hence 4 was 6.9 neutrons per cm? 
per sec. This may be compared with a current of 1.5 re- 
ported by v. Halban, Kowarski and Magat* at a lower 
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Fic. 1. Counting rate of recoil counter in presence of Ra-Be neutron 
source, as a function of voltage. No beta- or gamma-rays are counted 
below about 1770 volts. 
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elevation (3 m water). It must be remembered that this 
estimate is very rough. Further experiments and calcula- 
tions are in progress. 

This work has been partly supported by the Carnegie 
Institution of Washington. The author is indebted to M. D. 
Whitaker and New York University for the loan of the 
neutron source, and to E. T. Clarke for assistance with the 
observations. 

S. A. Korrr 

Bartol Research Foundation of the Franklin Institute, 


Swarthmore, Pennsylvania, 
October 30, 1939. 


1S. A. Korff and W. E. Danforth, Phys. Rev. 55, 980 (1939). 

2C. Kittel and G. Breit, Phys. Rev. 56, 747 (1939). 

3H. v. Halban, L. Kowarski and M. Magat, Comptes rendus 208, 
572 (1939). 





The Absorption of Mesotrons in Air and in 
Condensed Materials* 


It has been pointed out by several authors that the ab- 
sorption of mesotrons in air is considerably larger than the 
absorption by equal masses of condensed materials. This 
fact has been interpreted as evidence for a spontaneous 
decay of the mesotron. A lifetime of about 2X10~® sec. 
is required in order to account for the difference. 

The great theoretical importance of this conclusion 
justifies a careful investigation of possible alternative ex- 
planations of the observed difference in abs~rption. I have 
therefore considered the following effect which seems to 
explain the observations, at least to some extent, without 
assuming a decay of the mesotron. 

The ionization loss of energy by a fast particle passing 
through matter is partly due, as is well known, to close im- 
pacts between the particle and the material electrons; a 
large fraction of the loss, however, is due to impacts at 
distances greater than the atomic radius. For a mesotron 
with energy of the order of some billions of ev the ionizing 
effects of the particle can reach to distances of over 10‘ times 
the interatomic distances. 

In a rarefied gas the action of the field of the passing 
particle on every molecule is independent of the perturba- 
tion due to the surrounding molecules. This is no longer 
true for a condensed material in which the electric field of 
the passing particle is largely affected by the reaction of 
the electric polarization of the substance. A calculation of 
this effect based on the classical theory of the ionization 
loss shows that it is by no means negligible. 

Simple formulas can be obtained if the dielectric proper- 
ties of the medium are schematized by assuming all the 
electrons (7 per unit volume) to be elastically bound with 
the same frequency v. The dielectric constant for low fre- 
quencies would then be e=1-+-ne?/rmy?. With these 
assumptions it can be proved that the energy loss per unit 
path for a particle moving with velocity v is less than the 
loss calculated with the ordinary theory by the following 
amount: 
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The result is not essentially dependent on the special 
assumption as to the dielectric properties of the substance. 

While these formulas give a relatively unimportant 
change in the stopping power of gases and solids for slow 
particles like protons or a-particles, the difference in be- 
havior becomes rather large when the velocity approaches 
that of light. Let us consider for example the energy loss of 
a mesotron of 310° ev in two different media: A con- 
densed material for which we take «=2; and air for which 
we take «= 1.00054. Neglecting the present effect one would 
expect the energy loss to be approximately 2.3 Mev-cm?/g 
for both media. The reduction of loss due to the inter- 
action is negligible in the case of air; it is instead about 
0.5 Mev-cm?/g for the condensed substance. This reduces 
the loss in the latter case to only 1.8 Mev-cm?/g. 

The effect of this difference on the absorption of cosmic 
rays can be estimated if we assume the number of meso- 
trons with energy >W to vary as W~-*, The ratio of the 
mesotron intensity observed under equal amounts of air 
and of condensed materials should then be (2.3/1.8)!-°= 1.6. 
According to Ehmert the experimental value of the ratio is 
about 2. 

It seems therefore that the theory accounts for the order 
of magnitude of the effect even without any contribution 
from the decay of the mesotron. 

Whether all the effects, and especially the somewhat 
greater differences of absorption reported as results of 
observations with relatively thin absorbers, can be inter- 
preted on the outlined basis is doubtful. Indeed the the- 
oretical result seems to be near one-half of the experimental 
difference. But in any case the interactions between atoms 
represent an important factor to be taken into account in 
the interpretation of experiments of this type. 

I hope to be able to give soon the details of the theory and 
of its applications in a more extensive publication. 

ENRICO FERMI 

Pupin Physics Laboratories, 

Columbia University, 
New York, New York, 
October 1, 1939. 


* Publication assisted by the Ernest Kempton Adams Fund for 
Physical Research of Columbia University. 





Helium of Mass Three as an Agent in Nuclear Reactions 


Recently Alvarez and Cornog! have demonstrated the 
existence of stable He*. It seems worth while to call at- 
tention to the important applications this isotope may 
have in extending the scope of nuclear physics to the study 
of nuclear types not previously obtainable. 

One observes first that the average binding energy per 
particle in He® is very low (2.72 mMU), so that as a 
general rule processes induced by it will be exothermic. 
In addition, with large cyclotrons it is convenient to ob- 
tain beams of He***+ ions having exceedingly high energies 
so that the penetration of Coulomb potential barriers up to 
large atomic numbers is possible. 

Recent work at Princeton? on the systematic production 
and study of light odd nuclei having one more proton than 
neutron was undertaken because of the usefulness of these 








cial 
ice. 
ant 
low 


hes 
s of 
on- 
ich 
uld 
er- 
ut 


ces 


nic 


ler 
on 


or 


we 








LETTERS TO 


TABLE I. Some (He', ) reactions. Energy release in 
thousandths of a mass unit. 














TARGET PRODUCT Q Ilrancer PRODUCT Q 
H? Li* —4.1 || Al p29 6.1 
Lié Bs* —2.4 || Sits S30 0.3 
Li? Be 9.7 || Size 3! 5.5 
B10 N! 0.8 || Ps» CRs 5.1 
Cr On —1.5 || S2 A 0.6 
N' Fé —1.9 | $8 AS 4.6 
ow Nes —3.3 || CI Ks? 5.4 
Ne? * Mg 0.7 |} A%* Cas 0.1 
Mg Si 1.2 | 








* These nuclei are probably unstable against disintegration into 
heavy particles. The masses have been estimated by H. A. Bethe, 
Phys. Rev. 55, 434 (1939). 





nuclei in the understanding of the beta-decay process,’ and 
for the accurate estimate of the Coulomb energy difference 
between isobars. The available reactions leading to mem- 
bers of this series (the second column of Table IV in 
reference 4) are severely limited by the proton energies 
available. The use of He*, however, would seem to offer the 
possibility of extending the work far beyond the present 
limit. 

Of probably even greater interest is the feasibility of 
making members of the series of even nuclei having two 
more protons than neutrons (listed in the first column of 
Table IV of reference 4). No nucleus of this type has as yet 
been reported, although the sequence of these nuclei, 
along with the two isobaric series having, respectively, 
equal numbers of protons and neutrons, and two more 
neutrons than protons, is of the greatest importance in the 
study of the finer details of the like and unlike particle 
interactions. 

Of the possible reactions (He*, m), (He*®, H') (He’, He‘), 
(He’, H?) and (He’, y), the first three types may with con- 
siderable probability be expected to occur. Of these the 
most useful for the production of nuclei at present not 
known is (He?, 2), It has been rather well established from 
the reactions (He‘, 7), (H?, 2), (H!, ), (y, 2) and (mn, 2) 
that when sufficient energy is available for the emission of 
a neutron from the compound nucleus, the probability of 
this process will be large. It is difficult therefore to think of 
a reason why reactions of the type (He*, m) should not be 
highly probable. In Table I are listed some of the more 
interesting of such reactions together with the antici- 
pated‘ release of energy. Many of these product nuclei 
cannot be formed by any other process. 

It may be pointed out, too, that the reaction 
H*(He’, H')Het is expected to yield protons of extraordi- 
nary energy since the “Q”’ for this reaction is 19.73 mMU. 

It is of course unfortunate that the light isotope of 
helium exists with only a small abundance, but pre- 
sumably the methods of isotope concentration now in use 
will help to increase the obtainable beam intensities. 


WALTER H. BarKAS 
Institute for Advanced Study, 
Princeton, New Jersey, 
December 1, 1939. 


ae) W. Alvarez and R. Cornog, Phys. Rev. 56, 379 (1939); 56, 613 
(1939). 

2? White, Delsasso, Fox and Creutz, Phys. Rev. 56, 512 (1939). 

3 E. P. Wigner, Phys. Rev. 56, 519 (1939). 

4W. H. Barkas, Phys. Rev. 55, 691 (1939). 
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The Use of Radioactive Elements as Tracers in Physiology 


In a Letter to the Editor of The Physical Review' Dr. A. 
Barnett has raised a question as to whether the radiation 
from a radioactive tracer may exert a disturbing influence 
upon the physiological process which is being studied. 
In particular he refers to the measurement of permeability 
of the red blood cell, by means of radiosodium. Fortunately 
there exist some data on the effects of radiation on perme- 
ability, and an estimate of the magnitude of the effect can 
be made. It hardly need be remarked that it is unnecessary 
to consider the local effect of a particular radioactive atom 
upon a cell which happens to be close to it, first because 
the range of the beta-ray (or gamma-ray) is much greater 
than the size of the cell; second because if the atom emits 
radiation before entering the cell it is thereafter “dead” 
as far as detection is concerned. The effect upon permeabil- 
ity will therefore be the result of the generalized dose of 
ionization which is delivered to the material throughout 
its volume, and this may be spoken of in terms of r units. 
A large number of experiments upon induced changes in 
permeability by exposure to gamma-rays and x-rays have 
been summarized in a review article by Heilbrunn and 
Mazia.? Although the results vary considerably, it is to be 
noticed that in general no clear-cut effect occurs for doses 
less than a few hundred r units. What concentration of 
radiosodium will give a dose of, say, 500 r in 10 hours? If 
it is assumed that the beta-rays alone are absorbed in the 
material, the answer is about 10* disintegrations per cc per 
second, or the equivalent of about 10°? milligram of 
radium per cc, in terms of gamma-ray strength. It is 
ordinarily not necessary to use tracers in such high concen- 
tration as this. However, in exceptional cases a correction 
may have to be applied for the effect of volume irradiation 
of the material. 

As to the ability of biological material to discriminate 
between the radioactive and nonradioactive isotopes of a 
given element, a remark may be made. In the case of 
hydrogen it is quite clear that the isotopes behave differ- 
ently, chemically, and this is best illustrated by the fact 
that they can be separated by chemical means in the 
laboratory. Obviously this is the extreme case, in that the 
isotopes differ in weight by a factor two, or even three, if 
radiohydrogen is included. For the heavier elements the 
isotopic difference in chemical behavior is agreed to be 
negligible, as far as any simple reaction is concerned, but 
biochemists might not allow us to generalize this to in- 
clude the complicated and delicately balanced systems of 
reactions which-eccur in living material. One way in which 
we may hope to find a rather conclusive answer to this 
question is as follows. So far as the author is aware, no 
appreciable difference has been noted in the behavior of 
biological material toward the various stable isotopes of 
any element, with the possible exception of hydrogen. 
If such a difference does exist, a variation in the isotopic 
ratio in an element should be found in going from one 
compound of biological origin to another. An accurate test 
of this point is possible by means of the mass spectrograph. 
It follows, of course, that the difference, biochemically, 
between radioactive and nonradioactive isotopes will 
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certainly not be greater than the difference between the 
stable isotopes of a given element. Finally, it may be 
pointed out that the use of radioactive isomers as tracers 
would be free even from the above question, since they 
differ neither in weight nor in charge. It is almost certain 
that the number of known isomers will grow, and they 
may find a unique application in certain tracer problems. 
H. R. CRANE 


University of Michigan, 
Ann Arbor, Michigan, 
November 27, 1939. 


1A. Barnett, Phys. Rev. 56, 963 (1939). 
*L. V. Heilbrunn and Daniel Mazia, Biological Effects of Radiation, 
edited by B. J. Duggar (New York, 1936). 





Vitrification of Water 


The main reason for the general failure of attempts to 
obtain water in the vitreous (amorphous) state by a rapid 
cooling from the liquid state, seems to be the high rate of 
formation of crystallization nuclei and the high velocity of 
growth of the ice crystals. With the idea that a sufficiently 
rapid cooling would finally prevent the formation of ice, 
we tried various procedures which consisted essentially in 
immersing suddenly in liquid air preparations in which 
droplets of water were either sprayed with an atomizer on 
thin glass or mica sheets or deposited on glass plates by the 
condensation of steam. An examination of these droplets 
between crossed Nicols revealed a crystalline structure. 
The cooling velocity had evidently not been high enough. 

We finally succeeded in obtaining amorphous solid water 
by flattening small quantities of liquid water between two 
metal plates, cooled in liquid air, one of which was thrown 
toward the other at a high speed. A cylindrical stream of 
water 1 mm in diameter was made to flow vertically from 
a pipette. A brass disk 12 mm in diameter and 2 mm thick, 
mounted on a rod perpendicular to the plane of the disk, 
and cooled in liquid air immediately before being used, was 
fastened about 2 mm behind the stream, the direction of 
the latter being parallel to the plane of the disk. Another 
brass disk of the same size, mounted on a rod 12 cm long, 
and cooled in liquid air, was thrown against the first disk 
by the propelling action of the spring of a toy pistol. 
The pistol had previously been fastened in a steady posi- 
tion and adjusted so that the disk that it was to throw was 
parallel to the other disk and about 13 mm away from it. 
The disk thrown by the pistol was, then, at a distance of 
1 cm from the stream. 

The water flattened between the two disks takes the 
shape of a film a few microns thick. Pieces of this film, 
which sometimes have an area of 10 mm?, were transported 
on a cooled glasc slide into a specially built desiccated 
chamber placed on the stage of a polarizing microscope. 
These pieces are dark between crossed Nicols and stay 
dark for about 30 seconds. Then, with the rise in tempera- 
ture, one can observe a gradual reestablishment of light 
and, finally, when the temperature approaches the melting 
point of ice, the preparation is full of crystals large enough 
to be well observed under the microscope. Water has 
evidently undergone vitrification and devitrification. 

This experiment, of course, does not allow one to con- 
clude anything as to the possible formation of crystalline 
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nuclei smaller than the wave-length of the light used to 


detect them. 
B. J. Luyet 
Department of Biology, 
Saint Louis University, 
Saint Louis, Missouri, 
November 16, 1939. 





Radioactive Isotopes in Biology 


In a recent note! it has been suggested that the general 
body of biological work with artificially radioactive indi- 
cators, which is currently appearing in journals, is subject 
to the criticism that apparently no one has made any 
measurements of the effects of these radioactive isotopes on 
biological materials. As regards the general systemic effects 
of radioactive isotopes, Hamilton and Alles? have data to 
show that at the concentrations with which they worked, 
there was no effect on electrocardiographic tracings, 
respiratory rate, pulse, and blood pressure. Greenberg and 
Glazko’ also found no difference between active and in- 
active salts when used on isolated hearts, in vitro. Work 
from this laboratory‘ has shown that in the case of the 
unicellular alga Nitella, high activities (1-50 millicuries/ 
liter) of radioactive Na will produce a decrease in the rate 
of uptake of sodium by the cell. This rate of penetration is 
inversely proportional to the logarithm of the millicurie 
dose, and is almost certainly due to the radiation from the 
artificial isotope 1:Na* and not due to some differential 
penetration of inactive Na rather than active Na. That 
this is so is shown by the fact that the above results could 
be quantitatively duplicated with x-rays as a radiant 
source. Radioactive 19K“ could also be shown to produce 
the effect as stated above. 

When doses of below 1 millicurie per liter of radioactive 
Na were used, there was no longer any change in the rate of 
penetration with dosage, and it was concluded that below 
this value there was no effect. In most of the published 
experiments so far done the dosage has been below this 
limit, although it is true that there are papers describing 
work in which undoubtedly the radiation was producing 
an effect. In the case of the experiments of Cohn and 
Cohn’ who used radioactive Na, the writer is acquainted 
with the dosage used and it is certainly below the activity 
to produce any so far known biological effect. In regard to 
the criticism in the previous note! it appears there are 
several points on membrane structure and permeability 
in which more recent work has shown him to be in error. 
Limited space does not permit of their discussion. In con- 
cluding it seems desirable to emphasize that the radiations 
emitted from the various radioactive isotopes important 
in biology vary so markedly in their character that a 
different limit of biological effect will have to be made for 
each one. 

Lorin J. MULLINS 

Department of Zoology, 

University of California, 


Berkeley, California, 
November 17, 1939. 


1A. Barnett, Phys. Rev. 56, 963 (1939). 
2 J. G. Hamilton and G. Alles, Am. J. Physiol. 125, 2 (1939). 
3A. J. Glazko and D. M. Greenberg, Am. J. Physiol. 125, 2 (1939). 
4L. J. Mullins, J. Cell. Comp. Physiol. 14, 3 (1939). 
cnees) E. Cohn, and E. T. Cohn, Proc. Soc. Exp. Biol. Med. 41, 445-49 
9). 
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On Explosion Showers 


Let us consider a high energy collision between two 
particles of initial impulses | p:|=|p2|>>po=hc/l (where 
pi, P2 are referred to a barycentric frame: f1+2=0, and / 
is a universal length ~10-* cm). 

Heisenberg,’ starting from a half-classical theory of 
explosion showers, showed that the number of particles 
created by collision must be sensibly proportional to the 
initial energy E=2|pi|c~nhc/l. 

Our purpose is to examine whether it is permissible to 
apply the second quantization method and to assume the 
existence of a suitable relativistic interaction-operator in 
order to obtain some additional information about this 
kind of shower. We arrive at the same description of 
showers obtained by Heisenberg, and this seems to support 
the consistency of our assumption. Moreover our method 
leads to the result that the cross section for explosion 
showers decreases for n>1 and tends to 0 when E>, 

In analogy to Heisenberg’s Lagrangian we shall consider 
here, as a plausible and typical example, an interaction- 
operator G() which is a function of 


&=PO. 


We suppose also that this operator corresponds to some 
kind of short range forces, e.g. we consider G(é) of the 
type exp [— |£| ] or (+1) (for Yukawa forces), or other 
types used for “cutting off” purposes. 

The general term of interaction, which in the formalism 
of second quantization, corresponds to an explosion shower, 
is: 

bki*bke* - « -bkatanraps f els*(x,1)- *+ pkn®G (E)powte(x,t)dx, 
where gk, represent plane waves with impulses &,p, 
and bk,*ap, are the usual operators representing the crea- 


tion and the annihilation of particles. gt.p, are normalized 
in a volume Q, and thus assume the form: 


Q-hu, exp (i/h)(px— Et), 
where u, are only spin dependent factors. We have: 


G(E)Yowp2=G {| — po *[ (pit po)? —c- (E+ E2)*)} 
Xupiup, exp (t/h)[(pitp2)x— (Ei: +E2)t]. 


Applying the usual perturbation method in order to ob- 
tain the Fourier coefficients of the function x which 
represents a state: 


x(NkiNke- . -No1 e+ of) 
=c(Nk- ° -Nor- - +f) exp (-( h)( Ney E+ - )t] 


(with obvious signification of the symbols) we obtain: 


|c(1ky+ + 1knOp1Opat) |? 








2 
(t=) 


2n 1 ~ o 
wT Joe ekn*Gypwodx 


1 (=)'(=) (*) 
AE h' v1 L Un 
is the number of quantum states referred to unit energy 


interval. In a barycentric frame we have p:1+:=0 
=k, +k2.+---+k,. An easy calculation shows that: |c|? 


where 
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results independent of 2 and proportional to: 


k,tk:?- --k,2G [ (V7 
CPo 


In order to study how this probability varies with m and 
with the spatial distribution of vectors k,, we make use of 
the conservation law: 


Epi t+ Epo = Eki + Ekr+ +++ +Ekn 


and find that this probability is the maximum in the case:* 
[Ril ~|R2|~-++~| Rn] ~do. 

Thus is proportional to the initial energy and the 
frequency of explosion showers results reduced by the 
factor G for incident energies greater than cho. 

GLEB WATAGHIN 

Department of Physics, 

Sao Paulo University, 


Brazil, 
November 16, 1939. 


1W. Heisenberg, Zeits. f. Physik 110, 251 (1938); 113, 61 (1939). 

? Following a remark of Mario Schénberg, the operators of the type 
exp (A/Ko?) applied to a function give in the corresponding Fourier- 
development a cut-off factor. 

ea; Wataghin, Nature 142, 393 (1938); Comptes rendus, August 
(19. ‘ 





Ferromagnetic Anisotropy of Low Nickel Alloys of Iron 


Although the ferromagnetic anisotropy constant K; has 
been determined accurately for iron-nickel alloys contain- 
ing 50 percent or more of nickel, nothing at all is known 
about the low nickel alloys. The marked discontinuity in 
magnetic properties of 30 percent nickel makes it impossible 
to interpolate between iron and the high nickel alloys. 
Until the time when suitable low nickel single crystals of 
iron are grown, the only way of obtaining any idea as to 
the behavior of the anisotropy constant K, in the low 
nickel range is to make use of cold-rolled material. This is 
possible because the texture of cold-rolled iron-nickel al- 
loys has been carefully investigated by McLachlan and 
Davey,' using standard x-ray diffraction technique. They 
found that for a given total percentage of cold reduction 
in thickness, the texture was independent of the composi- 
tion up to at least 20 percent nickel. 

The magnetic torque curve of cold-rolled iron is similar 
to that of an iron single crystal disk having (100) in its 
plane.2 The exact shape and the amplitude of the poly- 
crystalline torque curve depend upon the nature of the 
deviations from the ideal (100) orientation and upon the 
fraction of material having a random texture. It is con- 
venient to speak of a texture factor, which is defined as the 
ratio of the amplitude of the predominant 40 term in the 
torque curve of the cold-rolled material to the corresponding 
amplitude for a single crystal of the same composition 
and volume For all practical purposes, the texture factor 
can be taken as the ratio of the average of the torque peaks 
of the polycrystalline material to that of the single crystal. 
In view of the work of McLachlan and Davey, it can be 
assumed that if iron and iron-nickel alloys are cold-rolled 
in the same manner, their texture factor will be the same, 
at least as a first approximation. The average values of the 
torque peaks will then be directly proportional to the 
corresponding values of K,. 
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Fic. 1. Ferromagnetic anisotropy constant for various low nickel alloys 
of iron. 


Pressed powder ingots of iron containing various amounts 
of nickel were hot-rolled from a thickness of 1 in. to 0.60 in., 
and then cold-rolled 97 percent to 0.020 in. Several 1-in. 
disks were stamped from each strip for the magnetic 
torque measurements. The torque peaks were measured at 
several different fields to permit of extrapolating them to 
their limiting values. The results, calculated on the as- 
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sumption that for iron K,=5.29X10° ergs/cm* and that 
the texture factor is constant for all the compositions 
studied, are shown in Fig. 1. A texture factor of 0.335 was 
found for the iron specimens. 

The probable accuracy of the results depends entirely 
upon the constancy of the texture factor. It is reasonable 
to assume that if this factor varies by 20 percent or there- 
abouts, the texture must have changed enough to affect 
the x-ray diffraction patterns. Since no such effect was 
found by McLachlan and Davey, a conservative estimate 
of the probable error of K, for the higher nickel range is not 
more than 20 percent. The error should be proportionately 
less for the alloys poorer in nickel, if we assume that, as the 
nickel content increases, the texture factor departs in a 
uniform manner from its value for iron. According to the 
straight line drawn through the points of Fig. 1, the aniso- 
tropy constant of iron containing 12 atomic percent of 
nickel is half that of pure iron. 

L. P. TARASOV 

Research Laboratory, 

General Electric Company, 
Schenectady, New York, 
December 1, 1939. 


1D. McLachlan, Jr., and W. P. Davey, Trans. Am. Soc. for Metals 
25, 1084 (1937). 

2R. M. Bozorth, Phys. Rev. 50, 1076 (1936). 

3L. P. Tarasov, Phys. Rev. 56, 1231 (1939). 
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Dielectric constants (see Dielectrics and dielectric proper- 
ties) 
Dielectrics and dielectric properties 
Absorption of charge, J. W. Simmons—848(A) 
Of NHs3, Nz and CO, at high frequency, G. W. Fox, A. H. 
Ryan—1132 
Breakdown in ionic crystals, H. Fréhlich—349; R. J. 
Seeger, E. Teller—352; R. C. Buehl, A. v. Hippel—941 
Permanent electrets, behavior of, W. M. Good, J. D. 
Stranathan—810 
Of Rochelle salt with D,O, A. N. Holden, G. T. Kohman, 
W. P. Mason, S. O. Morgan—378(L) 
Of water vapor at 42 megacycles, A. C. Tregidga—856(A) 
Diffusion 
Coefficient of thermal diffusion, A. O. Nier—1009 
Self-diffusion of Cu, J. Steigman, W. Shockley, F. C. Nix 
—13 
In solid state, R. P. Johnson—814 
Discharge of electricity in gases 
Anode spots in glow discharge, S. M. Rubens, J. E. 
Henderson—854(A) 
Arcing phenomena in Hg switches, C. G. Suits—215(A) 
Cathode sputtering, C. Starr—216(A) 
Chamber geometry; h. f. discharge, S. Githens, Jr.— 
845(A) 
Effect of pressure on positive point-to-plane discharge, 
H. C. Pollock, F. S. Cooper—170 
Formation of ions in cyclotron, R. R. Wilson—459 


Humidity and sparkover voltage of sphere gap, A. B. 
Lewis—846(A) 

Initial recombination, N. E. Bradbury—849(A) 

Magnetic effects in plasma, L. Tonks—360 

Temperatures in electrodeless discharge, M. S. McCay, 
H. P. Knauss, B. Archdeacon—847 (A) 

Townsend coefficients for Hz, D. H. Hale—858(A), 1199 

Transitions between glow and arc, W. G. Cady—123(L) 


Discharge of electricity in high vacua (see Thermionic 


emission) 


Disintegration of nucleus (see also Radioactivity) 


Of Be by electrons, G. B. Collins, B. Waldman, E. Guth 
—876 

Of B by deuterons, E. Pollard, W. L. Davidson, Jr., 
H. Schultz—206(A) 

Cross section for reaction H?+H*—+H!+ Hs; gas target, 
J. A. Van Allen, A. Ellett, D. S. Bayley—383(L) 

Delayed gamma-rays from U, J. C. Mouzon, R. D. 
Park—238 

Of deuterium by protons and p-n reactions, W. H. Bar- 
kas, M. G. White—288(L) 

Disintegration induced by cosmic radiation, I. Zlotowski 
—484(L) 

Dispersion formula for nuclear reactions, A. J. F. Siegert 
—750 

Energies released in Li? (p,a)He* and Li*(d,a)Het, S. K. 
Allison, L. S. Skaggs, N. M. Smith, Jr.—288(L), 
860(A); N. M. Smith, Jr.—548 

Energy released in Be®(p,d) Be’, L. S. Skaggs—24 

Fission products of Th, A. Langsdorf, Jr.—205(L) 

Fission of protactinium, A. v. Grosse, E. T. Booth, J. R. 
Dunning—382(L) 

y-radiation from F+H, C. C. Lauritsen, W. A. Fowler, 
T. Lauritsen—858(A) 

y-radiation from Li, W. A. Fowler, C. C. Lauritsen— 
841(L) 

y-radiation from N+H?, L. A. Pardue, E. R. Gaerttner 
—855(A) 

Heat of fission of U, M. C. Henderson—703(L) 

Mechanism of nuclear fission, N. Bohr, J. A. Wheeler— 
426 

Neutron production and absorption in U, H. L. Ander- 
son, E. Fermi, L. Szilard—284 

Neutrons from F by deuterons, T. W. Bonner—207(A) 

Of N™ and N® by deuterons, M. G. Holloway, B. L. 
Moore—705(L) 

Pair emission from F, W. A. Fowler, C. C. Lauritsen— 
840(L) 

Of protactinium, N. Bohr, J. A. Wheeler—1065(L) 

Protons from V, Cu, Mn and Sc, W. L. Davidson, Jr.— 
1061(L); 1062(L) 

Radiation from deuteron-deuteron reaction, M. H. 
Kanner, W. T. Harris—839(L) 

Reactions involving He*, W. H. Barkas—1242(L) 

Recoil of nucleus in beta-decay, H. R. Crane, J. Halpern 
—232 

Resonance in a-particles from F!*+H!, W. B. McLean, 
R. A. Becker—859(A) 

Short range protons in D-D reaction, F. E. Myers, R. D. 
Huntoon, C. G. Shull, C. M. Crenshaw—1104 
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Symmetrical distribution in B™(n,a), W. E. Good, R. D. 
Hill—288(L) 

Of U by neutrons; products, P. Abelson—1 

U and Th fission by fast neutrons, R. Ladenburg, M. H. 
Kanner, H. Barschall, C. C. Van Voorhis—168 

Dissociation 

Of hydrocarbons by electron impact, J. 
Bleakney—208(A), 256 

Of H.O vapor by electron impact, A. Hustrulid, M. M. 
Mann, J. T. Tate—208(A) 


Delfosse, W. 


Elasticity 
Adiabatic elastic moduli of NH;Cl, A. W. Lawson, R. 
Scheib—211(A) 
Compressibility of Al+Cu by x-ray measurements, R. 
B. Jacobs—211(A) 
Constants of single crystals of CusAu, S. Siegel—211(A) 
Internal friction; intercrystalline thermal currents, R. H. 
Randall, F. C. Rose, C. Zener—343 
Internal friction of Zn crystals, T. A. Read—211(A) 
Photoelasticity in three dimensions, R. Weller—216(A) 
Young’s modulus; new method, A. King—216(A) 
Electrical circuits 
Quenched primary of induction coil, W. P. Boynton— 
853(A) 
Electrical conductivity and resistance 
Change in a magnetic field, L. Davis, Jr.—93 
Of Mn compounds; effect of temperature, C. F. Squire— 
960(L) 
Metallic bridges between separated contacts, G. L. 
Pearson—471 
Of thin Pb films, M. G. Foster, A. Goetz—849(A) 
Electromagnetic theory 
Diffraction and radiation, S. A. Schelkunoff—308; J. A. 
Stratton, L. J. Chu—99, 316 
Energy fluctuations in field, P. Morrison—937 
Reflections from diffuse boundaries, C. D. Thomas, R. C. 
Colwell—1214 
Space charge and field waves in electron beam, S. Ramo 
—276 
Wide-angle interference; nature of elementary light 
sources, P. Selényi—477 
Electrons 
Proper energy of electron G. M. Volkoff—858(A) 
Structure of electric particles and the number 137; 
Born’s reciprocity, A. Landé—486(L) 
Electrons, scattering of (see Scattering of electrons, neu- 
trons and ions) 
Electrons, secondary 
By high energy electrons, W. A. Fowler, C. W. Sheppard 
—859(A) 
From metals; positive ion bombardment, J. B. Fisk— 
846(A) 
From Ni and Co, D. E. Wooldridge—1062(L) 
Theory of, D. E. Wooldridge—562 
Electro-optical effects (see also Photoelectric effect and 
properties) 
New effect, H. Mueller, B. W. Sakmann—615(L) 
Energy states of atoms (see Atomic structure) 
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Energy states of nucleus (see Nuclear structure) 
Errata 
Interaction of configurations: sd-p*, R. F. Bacher 
385(L) 
M, and M, x-ray absorption edges of lead, J. W. Mc- 
Grath—1068(L) 
Paschen-Back effect. VI. The spectrum of neon, J. B. 
Green, J. A. Peoples, Jr.—385(L) 


Velocity of radio waves in air, G. H. Brown—385(L) 


Field currents 
Fields involved, J. E. Henderson, K. V. MacKenzie— 
857(A) 
Range and validity of equation, F. R. Abbott, J. E. 
Henderson—113, 853(A) 
Field theory 
Self-energy and electromagnetic field of electron, V. F. 
Weisskopf—72 
Films, properties of 
Conductivity of thin Pb films, M. G. Foster, A. Goetz— 
849(A) 
Structure of very thin films, L. H. Germer—58 
Fine structure (see Spectra, etc. and Hyperfine structure) 
Fluorescence 
Of OH molecules, E. R. Lyman—466 


Gravitation 
Continued gravitational contraction, J. R. Oppenheimer, 
H. Snyder—455 


Heat of combustion 
Mileage of gasoline of different grades, J. B. Derieux— 
845(A) 
High voltage tubes and machines (see Methods and 
instruments) 
Hydrodynamics 
Streamline flow, D. Kirkham—852(A) 
Hyperfine structure (see also Nuclear moments and spin) 
Of iridium lines, L. Sibaiya—768 
Metastable *P3,2 state of In™, D. R. Hamilton—30 


Instruments (see Methods and instruments) 
Ions, mobility 
Of atmospheric ions, E. A. Yunker—855(A) 
Hg ions in Hg vapor, K. H. Kingdon, E. J. Lawton 
215(A) 
Of K ions at high fields, A. V. Hershey—908 
Theory of; for high velocity ions, A. V. Hershey—916 
Isomers, nuclear (see Radioactivity) 
Isotopes 
Concentration of C by thermal diffusion, W. W. Wat- 
son—703(L); A. O. Nier—1009 
Concentration by centrifuge, J. W. Beams, C. Skarstrom 
—266 
He’ and H?, L. W. Alvarez, R. Cornog—379(L), 613(L) 
Mass of Li®, Li?, Be*, Be®, N. M. Smith, Jr.—548 





1260 


Isotopes (continued) 
Separation of Br isotopes by centrifugation, R. F. 
Humphreys—684 
Vapor pressures, heats of vaporization and melting points 
of N“ and N® ammonia, H. G. Thode—214(A) 


Kinetic theory of gases 

Convection and conduction, I. Brody, F. de Korésy— 
217(A) 

Fugacity of CO, W. E. Deming, L. S. Deming—108 

Transport phenomena in mixture of gases, E. J. Hellund, 
E. A. Uehling—818, 851(A) 

Turbulence in convection between concentric cylinders, 
L. Onsager, W. W. Watson—474 


Liquids (see also Hydrodynamics) 
Atomic distribution curves, C. A. Coulson, G. S. Rush- 
brooke—1216 
Contact electric phenomena in oils, S. Kyropoulos— 
853(A) 
Structure of liquid K, N. S. Gingrich, C. N. Wall—336 
Velocity of sound in liquid He, J. C. Findlay, A. Pitt, H. 
Grayson Smith, J. O. Wilhelm—122(L) 
Luminescence 
Of sulfide phosphors, R. P. Johnson—213(A) 
Of synthetic halite crystals, B. E. Cohn—856(A) 


Magnetic properties 
Anisotropy of Fe and Fe-Si alloys, L. P. Tarasov—1231 
Anisotropy of Ni, H. J. Williams, R. M. Bozorth— 
837(L) 
Anisotropy of Ni-Fe alloys, L. P. Tarasov—1245(L) 
Antiferromagnetism in manganous compounds, C. F. 
Squire—922 
Of electrolytically polished silicon-iron crystals, W. C. 
Elmore—210(A) 
Ferromagnetic anisotropy, L. P. Tarasov—1224 
Ferromagnetism in austenite, L. P. Tarasov, E. R. 
Parker—379(L) 
Hysteresis losses in steel, L. P. Tarasov—1170(L) 
Of solid solutions, F. Bitter, A. R. Kaufmann—1044 
Surface effects of impurities, F. W. Constant, J. M. 
Formwalt—846(A) 
Susceptibility of Cu-Ni and Zn-Ni alloys, M. A. Wheeler 
—1137 
Test for superstructure in permalloy, S. T. Pan—933 
Transition in paramagnetic salts, C. F. Squire—211(A) 
Magneto-optical effects 
Observations at low light intensities, W. N. Latimer, 
H. A. Young—963(L) 
Mass defects (see Isotopes) 
Mechanics, quantum-atomic structure and spectra 
Wave functions of relativistic hydrogenic atom, L. 
Davis, Jr.—851(A) 
Mechanics, quantum—general 
Dispersion formula for nuclear reactions, A. J. F. Siegert 
—750 
Double beta-disintegration, W. H. Furry—1184 
Fluctuations in electromagnetic field, P. Morrison— 
849(A), 937 
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Klein’s fifth dimension, W. Band—204(L) 
Proper energy of electron, G. M. Volkoff-—-858(A) 
Weinstein’s theorem, G. Horvay—214(A) 
Mechanics, statistical 
Fluctuations and entropy, E. C. Kemble—1013 
Quantum-mechanical basis, E. C. Kemble—1146 
Mesotrons 
Absorption in matter, E. Fermi—1242(L) 
8-ray selection rules; theory of mesotron, M. Schénberg 
—612(L) 
Disintegration of, B. Rossi, H. Van N. Hilberry, J. B. 
Hoag—837(L) 
Intensity of mesotrons; function of altitude, M. Schein, 
W. P. Jesse, E. O. Wollan—613(L) 
Lifetime and 8-decay, R. Serber—1065(L) 
Saturation of forces, E. Wigner, L. Eisenbud—214(A) 
Slow mesotrons in cosmic radiation, C. G. Montgomery, 
W. E. Ramsey, D. H. Cowie, D. D. Montgomery—635 
Spinor equations for mesotron, A. H. Taub—799, 
853(A); H. Fréhlich, W. Heitler, B. Kahn—961(L) 
Theory; Coulomb law for proton, W. E. Lamb, Jr.— 
384(L) 
Metals (see Crystalline state) 
Meteorology 
Heat transfer in the atmosphere, W. M. Elsasser—855(A) 
Radiation from sky, C. M. Heck—848(A) 
Zodiacal light during solar eclipse, E. O. 
960(L) 
Methods and instruments 
Analysis of Rb in plant sap, L. A. Strait—855(A) 
Arc source for cyclotron, E. McMillan, W. W. Salisbury 
—836(L) 
Bent rocksalt crystal spectrograph, P. H. Abelson—753 
Be targets for nuclear research, M. C. Henderson— 
207(A) 
Bragg spectrometer for low temperature, A. Goetz, A. 
Dember—857 (A) 
Cathode sputtering, C. Starr—216(A) 
Centrifuge; analytical, J. W. Beams, R. E. Fox—848(A) 
Clean-up of H; by Ba, G. W. Johnson, W. A. Hane, P. A. 
Anderson—852(A) 
Cloud chamber expansion ratio and drop count deter- 
mination, D. Bagley—851(A) 
Concentration of isotopes by centrifuge method, J. W. 
Beams, C. Skarstrom—266 
Dental fluoroscope and shield, C. R. Fountain—844(A) 
Differential manometer. W. Hurst—846(A) 
Dispenser cathode; cathode for gaseous discharge tubes, 
A. W. Hull—86 
Double Fizeau toothed wheel, H. E. Ives—215(A) 
Droplets in expansion chamber, growth of, R. M. Langer 
—851(A) 
Electrolytic conductivity at high frequencies, G. Stipe, 
Jr., J. H. Purks, Jr.—843(A) 
Electron counter, R. R. Wilson, D. R. Corson—851 (A) 
Electron microscope, G. A. Morton, E. G. Ramberg— 
705(L); E. F. Burton, J. Hillier, A. Prebus—1171(L) 
Electrostatic generator and ion source, I. A. Getting, 
J. B. Fisk, H. Vogt—213(A), 1098 
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Electrostatic voltmeter, F. H. Nadig, J. L. Bohn— 
216(A) 

Fast neutron counters, S. A. Korff—1241(L) 

Fast neutrons, beam of, P. C. Aebersold—714 

Ferromagnetic impurities; method for detecting, F. W. 
Constant, J. M. Formwalt—373 

Films on interferometer mirrors, W. V. Houston, S. 
Rubin—852(A) 

Geiger counter, modulation of, J. N..Shive—579 

Heat conductivity measurements, C. C. Bidwell—594 

Hertz pumps, modified, K. T. Bainbridge, H. G. Vogt— 
213(A) 

High pressure cloud chamber, W. E. Stephens, W. E. 
Shoupp, R. O. Haxby, W. H. Wells—209(A) 

Illumination of cloud chamber, C. E. Nielsen—855(A) 

Initial performance of 60-inch cyclotron, E. O. Lawrence, 
L. W. Alvarez, W. M. Brobeck, D. Cooksey, D. R. 
Corson, E. M. McMillan, W. W. Salisbury, R. L. 
Thornton—124(L) 

Ions in cyclotron, R. R. Wilson—459 

Low light intensities, observations of, W. N. Latimer, 
H. A. Young—963(L) 

Manometer of small volume, C. Kenty—215(A) 

Melting points by transmitted light, D. Williams, B. 
Hinton—844(A) 

Nomogram for atomic constants, R. A. Beth—208(A) 

Nuclear physics, apparatus for, J. C. Mouzon—846(A) 

Photoelasticity in three dimensions, R. Weller—216(A) 

Photoelectric spectrophotometer, W. C. Bosch, B. B. 
Brown—846(A) 

Raman spectra, excitation of, J. S. Kirby-Smith—844(A) 

Resistance voltmeter for high voltages, W. B. McLean, 
R. A. Becker—859(A) 

Resonance radiation, demonstration of, R. W. Wood— 
1172(L) 

Square-wave voltages at high frequencies, W. H. Fenn— 
850(A) 

Statistical significance of rate differences, F. T. Holmes— 
844(A) 

Temperature from nitrogen bands, O. S. Duffendack, 
Kuang-Tseng Chao—176 

Time distribution of random events, A. Ruark—1165 

Ultracentrifugation of a bacteriophage, T. Davis— 
845(A) 

Use of radioactive tracers, A. Barnett—963(L); H. R. 
Crane—1243(L); L. J. Mullins—1244(L) 

Viscosity measurements, D. R. Inglis—1040 

Work functions, determination of, G. M. Fleming, J. E. 
Henderson—853(A) 

Work functions, measurement: of, W. H. Goss, J. E. 
Henderson—857 (A) 

Young’s modulus; new method, A. King—216(A) 

Mobility of ions (see Ions, mobility) 


Molecular structure and constants (see also Spectra, 


molecular and Raman spectra) 

Coriolis coupling terms in polyatomic molecules, H. A. 
Jahn—680 

Double minima problem of NH; molecule, R. F. Haupt, 
E. Teller—208(A) 


Electronic states of diatomic carbon, R. S. Mulliken— 
778 

Energies of pentatomic molecules, W. H. Shaffer, H. H. 
Nielsen, L. H. Thomas—895 

Force constants of NO: groups, D. Williams, L. Decherd 
—847 (A) 

Forces in molecules; formulas for, R. P. Feynman—340 

Fundamental vibration of NH3, G. B. B. M. Sutherland 
—836(L) 

Vibration of CH;Cl at 7.54, A. H. Nielsen, H. H. Nielsen 
—847(A) 

Vibration-rotation energies of nonlinear triatomic 
molecule, W. H. Shaffer, H. H. Nielsen—188 


Neutrons 


Collimation of fast neutrons, R. F. Bacher, D. C. Swan- 
son—483(L) 

Fast neutrons, beam of, P. C. Aebersold—714 

Production and absorption in U, H. L. Anderson, E. 
Fermi, L. Szilard—284 

Recoil of nucleus in 8-decay, H. R. Crane, J. Halpern— 
232 


Nuclear moments and spin (see also Hyperfine structure) 


Of :3Al*?, S. Millman, P. Kusch—214(A), 303 

Of B isotopes, S. Millman, P. Kusch, I. I. Rabi—165 , 

Of C8, C. H. Townes—850(A); C. H. Townes, W. R. 
Smythe—1210 

Electric and magnetic moments of Li and N™, M. E. 
Rose—1064(L) 

Electric quadrupole moment of In“, D. R. Hamilton 
—30 

Of ;B** and ,B", P. Kusch, S. Millman—213(A) 

Of Ir, L. Sibaiya—768 

Of proton and deuteron, J. M. B. Kellogg, I. I. Rabi, 
N. F. Ramsey, Jr., J. R. Zacharias—213(A), 728 

Of the proton, neutron and deuteron, E. E. Witmer— 
214(A) 

Regularity in moments, E. E. Witmer—203(L) 

Of Rb and Cl isotopes, P. Kusch, S. Millman—527 


Nuclear structure (see also Disintegration of nucleus; 


Energy states of nucleus) 

Binding energies of H* and He®, K. Way—556 

Binding energy of H’, F. W. Brown—1107 

Coupling conditions in light nuclei, E. P. Wigner—519 

Deviation from Coulomb law, W. E. Lamb, Jr.—384(L); 
H. Froéhlich, W. Heitler, B. Kahn—961(L) 

Levels in B®, P. G. Kruger, F. W. Stallmann, W. E. 
Schoupp—297 

Mesotron forces; state of He®, S. M. Dancoff—384(L) 

Nuclear five-body problem, W. A. Tyrrell, Jr.—250 

Pair emission from Ne, J. R. Oppenheimer, J. S. 
Schwinger—1066(L) , 

Phase shifts for error well, H. M. Thaxton, A. M. Mon- 
roe—616(L) 

Phase shifts for square well, A. M. Monroe, H. M. 
Thaxton—616(L) 

Probabilities in double 8-disintegration, W. H. Furry— 
1184 

Spin-dependence of nuclear forces, C. L. Critchfield—540 

Spin-orbit coupling in light nuclei, D. R. Inglis—1175 
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Nuclear structure (continued) 
Of S*, E. Pollard—961(L) 
Theory of nuclear forces, E. P. Wigner, C. L. Critchfield, 
E. Teller—530 
Virtual level of He’, S. M. Dancoff—852(A) 


Optical constants and properties 
Infra-red pleochroism in CH, groups, J. W. Ellis, J. Bath 
—859(A) 
Optical instruments (see Methods and instruments) 


Paschen-Back effect 
In Ne, J. B. Green, J. A. Peoples, Jr.—54 
Theory for intermediate coupling, J. B. Green, J. F. 
Eichelberger—S1 
Photoconductivity 
Of tartaric crystals, J. J. Brady, J. D. Todd—859(A) 
Of willemite at low temperatures, R. Hofstadter, R. C. 
Herman—212(A) 
Photoelectric effect and properties; cells 
Magnetic field effect, N. K. Olson—210(A) 
Of U, W. L. Hole, R. W. Wright—785 





Work functions, measurement of, W. H. Goss, J. E. 
Henderson—857 (A) 
Work functions Ag crystal faces, H. E. Farnsworth, 


R. P. Winch—1067(L) 

Photography 

Graininess of photographic emulsions, W. O. 
A. Goetz—850(A) 
Sensitivity of emulsions for electrons, L. Marton—290(L) 

Photometry 
Photoelectric spectrophotometer, W. C. Bosch, B. B. 

Brown—846(A) 
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Radioactivity (see also Disintegration of nucleus) 
Of Al bombarded by H’, G. Kuerti, S. N. Van Voorhis— 
614(L) 
B-decay of light nuclei, B. O. Grénblom—508 
B- and y-rays from Cu® and Eu, A. W. Tyler—125 
B- and y-rays, coincidences from indium, L. M. Langer, 
A. C. G. Mitchell, P. W. McDaniel—380(L) 
B-radiation of As7*, G. L. Weil, W. H. Barkas—485(L) 
B-ray spectra of Mg??, Cu®; isomers of Rh!, E. C. Crit- 
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B-ray spectra of P, Na and Co, J. L. Lawson—131 
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Kanner, W. T. Harris—207(A) 
Double 8-disintegration, W. H. Furry—1184 
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Emission of neutrons by U, W. H. Zinn, L. Szilard—619 

In Eu, K. Fajans, D. W. Stewart—625 
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Half-life of B®, R. A. Becker, E. R. Gaerttner—854(A) 

Of In, L. D. P. King, W. J. Henderson—1169(L) 
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Internally converted y-rays, G. E. 
Creary—863 

Of I, 8-spectra, G. F. Tape—965 

Isomer of Sr8?, L. A. DuBridge, J. Marshall—706(L) 

Isomerism in Zn, J. W. Kennedy, G. T. Seaborg, E. 
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Positron spectra of N" and Na*, F. Oppenheimer, E. P. 
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Products of disintegration of U by neutrons, P. Abel- 
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287(L) 
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L. A. Delsasso, J. G. Fox, E. C. Creutz—207(A), 512 
Of Si?’ incorrectly assigned, H. A. Bethe, W. J. Hender- 
son—1060(L) 
In Sr and Y, D. W. Stewart—629 
Test for Mg**(a,p) Al, W. J. Henderson, R. L. Doran— 
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Time distribution of random events, A. Ruark—1165 
Use of “‘tracers,”” A. Barnett—963(L); H. R. Crane— 
1243(L); L. J. Mullins—1244(L) 
X-rays associated with radioactive decay, P. H. Abelson 
—753 
Raman spectra 
Of aliphatic molecules in water solution, R. C. William- 
son—848(A) 
In samarium nitrate, L. Sibaiva, H. S. Venkataramiah- 
381(L) 
Relativity 
Electrostatic analog to gravitational red shift, H. T 
Drill—184 
Rotating charged concentric spheres, L. I. Schiff—850(A) 
Wave functions of relativistic hydrogenic atom, L. 
Davis, Jr.—186 
Resonance radiation 
Demonstration of, R. W. Wood—1172(L) 


Scattering of electrons, neutrons, and ions (see also Elec- 

trons, scattering of; Electron diffraction) 

Of a-particles by N, G. Brubaker—207(A), 1181 

Of cathode rays, J. G. Trump, R. J. Van de Graaff, 
R. W. Cloud—215(A) 

Collision cross sections for D-D neutrons, W. H. Zinn, 
S. Seely, V. W. Cohen—260 

Cross sections of C, N, Na, Mg, and Al for neutrons, M. 
MacPhail, J. Giarratana—207 (A) 








19 








ANALYTIC SUBJECT INDEX 1263 


D-D, H-He, D-He scattering, N. P. Heydenburg, R. B. 
Roberts—1092 
Of electrons in Pb and Al; energy losses, M. M. Slawsky, 
H. R. Crane—1203 
Of fast electrons, N. L. Oleson, K. T. Chao, J. Halpern, 
H. R. Crane—482(L), 1171(L) 
Of fast electrons by heavy elements, J. H. Bartlett, R. E. 
Watson—612(L) 
Of fast neutrons, E. R. Gaerttner, L. A. Pardue, J. F. 
Streib—856(A) 
Multiple scattering of charged particles, S. Goudsmit, 
J. L. Saunderson—122(L) 
Multiple scattering of electrons and positrons, C. W. 
Sheppard, W. A. Fowler—849(A) 
Of neutrons from the C+D reaction, E. Amaldi, D. 
Bocciarelli, F. Rasetti, G. C. Trabacchi—881 
Of neutrons by protons, C. Kittel, G. Breit—744 
Phase shift for proton-proton scattering, H. M. Thaxton, 
L. E. Hoisington—1194 
Photoneutrons from deuterium, T. Goloborodko, A. 
Leipunski—891 
Proton-proton and proton-neutron interactions, F. W. 
Brown, M. S. Plesset—841(L) 
Proton-proton scattering, L. E. Hoisington, S. S. Share, 
G. Breit—884 
Proton-proton scattering, E. Creutz—893 
Of protons by protons, N. P. Heydenburg, L. R. Haf- 
stad, M. A. Tuve—1078 
Reflection of electrons by metals, L. A. MacColl—699 
Secondary electrons (see Electrons, secondary) 
Specific heat 
Heat capacity at Curie point, C. C. Stephenson, J. G. 
Hooley—121(L) 
Of monatomic liquid, W. J. Archibald—926 
Of NHs;, R. F. Haupt, E. Teller—208(A) 
Of Ni-Cr alloys, H. Grover, J. Hutzenlaub—212(A) 
Spectra, absorption 
Of AIF, G. D. Rochester—305 
Of CN, J. U. White—853(A) 
Of D.O, P. E. Shearin, E. K. Plyler—845(A) 
Of heavy water vapor molecule, N. Fuson, H. M. 
Randall, D. M. Dennison—982 
Of methyl chloride at 7.4n, A. H. Nielsen, H. H. Nielsen 
—274 
Of NO, A. H. Nielsen, W. Gordy—781 
Of NO, rotation-vibration band, R. H. Gillette, E. 
Eyster—1113 
Of some sugars, near infra-red, E. S. Barr, C. H. Chris- 
man—845(A) 
Water band at 4.7y; effects of alkali halides, T. N. Gau- 
tier, D. Williams—616(L) 
Spectra, atomic 
Calculation of two-electron configurations, A. F. 
Stevenson—586 
D*p configuration of the Ni-like and Pd-like spectra, 
A. T. Goble—209(A) 
Intensities of singlet-singlet and singlet-triplet transi- 
tions, G. W. King, J. H. Van Vieck—464 
Isotope shift in Mg, J. P. Vinti—1120 








Of Ne V, F. W. Paul—1067(L) 
Of Ne, Paschen-Back effect, J. B. Green, J. A. Peoples, 
Jr.—54, 385(L) 
Probabilities of forbidden lines, S. Pasternack—854(A) 
Of Ag, A. G. Shenstone—209(A) 
Of Th III, R. J. Lang—272 
Wave functions of relativistic hydrogenic atom, L. 
Davis, Jr.—186 
Of Y V and Zr VI, F. W. Paul, W. A. Rense—1110 
Spectra, molecular (see also Molecular structure and 
constants) 
Energies of pentatomic molecules, W. H. Shaffer, H. H. 
Nielsen, L. H. Thomas—895, 1051 
Forbidden transitions in Ne, J. Kaplan—858(A) 
Further resolution of NH; bands 10u and 16z, E. F. 
Barker, H. Y. Sheng—854(A) 
Interferometer wave-lengths of lines in spectrum of Ha, 
L. S. Combes, R. M. Frye, N. A. Kent—678 
Normal modes of long chain hydrocarbons, L. H. 
Thomas, S. E. Whitcomb—383(L) 
Rotation temperatures of CO, M. S. McCay, H. P. 
Knauss, B. Archdeacon—847(A) 
Temperature from nitrogen bands, O. S. Duffendack, 
Kuang-Tseng Chao—176 
Transition of C., O. G. Landsverk—769 
Spectroscopy, technique 
Analysis of Rb in plant sap, L. A. Strait—855(A) 
Stark effect 
Broadening of hydrogen lines, L. Spitzer, Jr.—39 
Surface tension 
Temperate dependence, ridges in liquid surface, A. V. 
Hershey—204(L) 


Thermal conductivity 
In gases, I. Brody, F. de K6résy—217(A) 
Heat transfer in atmosphere, W. M. Elsasser—855(A) 
Of Zn and silocel, C. C. Bidwell—594 
Thermionic emission of electrons; Emitting surfaces 
Determination of work functions, G. M. Fleming, J. E. 
Henderson—853(A) 
Dispenser cathode; cathode for gaseous discharge tubes, 
A. W. Hull—86 
Periodic deviation from Schottky line, D. Turnbull, T. E. 
Phipps—663; H. M. Mott-Smith—668; R. L. E. 
Seifert, T. E. Phipps—652 
Of surfaces of tungsten single crystal, S. T. Martin—947 
Of U, W. L. Hole, R. W. Wright—785 
Work function of W, J. G. Potter—216(A) 
Thermodynamics 
Fluctuations and entropy, E. C. Kemble—1013 
Fugacity of CO2, W. E. Deming, L. S. Deming—108 
Statistical mechanics, E. C. Kemble—1146 


Van der Waal forces 

Potential'in He, H. Margenau—1000 
Vapor pressure 

Of N“ and N® ammonia, H. G. Thode—214(A) 
Viscosity 

Of air, J. A. Bearden—1023 
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Viscosity (continued) 
Almost-coaxial cylinders; viscous force, D. R. Inglis— 
1040 
Electric field; effect of, A. A. Dixon—847(A) 
Structural viscosity, F. Breazeale—845(A) 


Work function (see Thermionic emission; Photoelectric 
effect) 


X-rays, absorption 
K absorption edges of Fe to Ge, W. W. Beeman, H. 
Friedman—392 
M, and M; edges of Pb, J. W. McGrath—765, 1068(L) 
M, and M, edges of Au in pure metal and in Au-Cu 
alloy, J. W. McGrath—137 
Resonance lines in K spectrum, L. G. Parratt—295 
X-rays, diffraction, scattering, reflection, refraction and 
polarization 
Diffraction pattern of calcite, F. Miller, Jr —757 
Diffraction by small crystalline particles, A. L. Patterson 
—972 
Diffuse scattering at small angles, G. G. Harvey—247 
Dispersion and atomic electron cloud distortion in Zn, 
E. M. McNatt—406 
Fluorescing coefficient and orientation of crystal, O. J. 
Baltzer, G. E. M. Jauncey—289(L) 
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Reflections from calcite between 300° and 4.2°K, A. 
Dember, A. Goetz—857(A) 

Scherrer formula, particle size determination, A. L. 
Patterson—978 

Theory of diffuse scattering, G. E. M. Jauncey—644 

Total scattering from crystals, G. G. Harvey—242 

X-rays, emission (see also X-rays, spectra and spectro- 

scopy, etc.) 

Fine structure in K lines of Cu, S. T. Stephenson— 
856(A) 

Intensities of lines of W and Pt, W. H. Kliever—387 

K auger transitions, C. J. Burbank—142; W. Veith, 
P. Kirkpatrick—705(L) 

Lines from double ionization of L-shell, R. D. Richt- 
myer—857 (A) 

Spin doublet splittings, R. F. Christy, J. M. Keller— 
856(A) 

Theory of lines LL-LM, R. D. Richtmyer—146 

Weak MoL lines, W. Veith, P. Kirkpatrick—705(L) 

X-rays, tubes, apparatus 
Dental fluoroscope and shield, C. R. Fountain—844(A) 


Zeeman effect 
Hollow cathode source, R. A. Fisher, A. S. Fry—675 
In the hyperfine structure of I, A. S. Fry, R. A. Fisher— 
669 
Of Hg resonance radiation, E. H. Collins—48 
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